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Part  One 

THEORY  AND  EXAMPLES 


I.  INTRODUCTION 


The  scattering  of  electromagnetic  waves  by  a  body  can  be  controlled 
by  impedance  loading  of  the  body,  either  by  lumped  loads  or  by  continuous 
loads.  The  general  problem  has  been  discussed  by  Schindler,  Mack,  and 
Blacksmith  [1].  Methods  are  available  for  the  analysis  of  lump  loaded 
scatterers  [2]  and  of  continuously  loaded  scatterers  [3].  A  synthesis 
procedure  for  controlling  the  electromagnetic  scattering  by  continuously 
loading  the  body  is  also  available  [4].  It  makes  use  of  the  theory  of 
characteristic  modes  of  a  conducting  body  [5,6].  Similar  modes  for  an  N-port 
scatterer  have  been  defined  and  used  for  the  analysis  of  lump  loaded  scat¬ 
terers  [7]. 

This  report  gives  methods  for  synthesizing  loaded  N-port  scatterers 
to  obtain  desired  scattering  patterns  or  to  optimize  the  radar  cross 
section.  The  geometry  of  the  scatterer  is  assumed  known  and  the  loads  are 
chosen  to  control  the  scattering.  The  theory  makes  use  of  the  characteristic 
modes  of  loaded  N-port  scatterers.  By  a  procedure  called  modal  synthesis  we 
can  make  any  set  of  real  port  currents  a  resonant  mode  of  the  scatterer.  If 
this  resonant  mode  is  the  principal  contributor  to  the  scattered  field,  then 
its  field  pattern  will  be  essentially  the  scattering  pattern  of  the  loaded 
scatterer.  To  synthesize  a  desired  scattering  pattern,  we  first  determine 
the  real  port  current  whose  field  pattern  approximate*  the  desired  pattern, 
and  then  we  resonate  it.  The  analogous  procedure  is  used  to  optimize  the 
scattering  cross  section. 
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II.  BASIC  THEORY 


Tha  analyst*  of  loadad  N-port  acattarers  is  dlscuaaad  in  raferanca 
[7].  Tha  thaory  is  expressed  equally  wall  In  terms  of  either  impedance 
parameters  or  admittance  parameters.  A  summary  of  tha  basic  theory  is 
given  here. 

Consider  an  N-port  scatterar  loadad  by  N  lumped  impedances ,  or, 
more  generally,  by  an  N-port  lead  network.  The  total  field  £  is  the  sum 
of  the  impressed  field  E1,  due  to  sources  axtarnel  to  the  acatterer,  plus 
a  scattered  field  E®,  due  to  current  Induced  on  the  scatterer.  In  the 
open-circuit  impedance  formulation,  the  scattered  field  is  given  by  [7] 

£*  -  ir  -  r  t*s  *  v^°c  «> 

Here  [Zg ]  and  [Z^]  are  the  port  impedance  matrices  of  the  scatterer  and 
load,  respectively,  V  is  a  column  matrix  of  the  open  circuit  port  voltages, 

A/QC 

E  is  a  row  matrix  of  the  fields  radiated  by  unit  current  sources  at  the 
oc 

ports,  and  Eg  is  the  field  scattered  when  all  ports  are  open  circuited. 

For  a  modal  analysis,  we  change  the  basis  to  one  which  diagonalizes 

[Z_  +  ZT],  In  particular ,  we  choose  as  a  basis  the  eigenvectors  I  of 
5  L  n 

fX]t  -  A  [R]t  (2) 

n  n  n 

where  [X]  and  [R]  are  the  Hermitlan  parts  of  [Zg  +  Z^],  assumed  to  be 
symmetric.  The  advantages  of  this  choice  are: 

(a)  All  eigenvalues  A  and  eigenvectors  t  are  real.  (More 

n  n 

generally,  the  Ir  can  be  equiphasal.) 

(b)  The  elgencurrents  1  form  an  orthonormal  set  with  weight  [R], 

n 

that  Is 

I  [R]f  -  6 
m  n  mn 


(3) 
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ic»  The  elgancurrents  1  fora  an  orthogonal  aat  with  raapact  to 

n 


(X!,  that  la 


T  [x)l  -  «  x 


(4) 


(d)  For  a  loaa-fraa  loaded  acatterer,  the  elgenflelda  E(I  )  form 

n 

an  orthogonal  set  over  the  aphera  at  Infinity,  that  la 


E*(!  )  •  E(l)ds 
*“  m  n 


6 

on 


(5) 


Here  E(I  )  denotes  the  field  radiated  when  the  eigenvector  1  exists  at 
n  n 

the  scatterer  ports. 

When  the  eigencurrents  1  are  chosen  as  a  basis,  the  solution  (1) 

n 

reduces  to 


A* 


N 


<N/  -*OC 

I  V 
n 


L.  1+jA 

n-1  J  n 


(6) 


Many  examples  of  the  use  of  (6)  for  analysis  are  given  in  reference  [7], 

If  the  scatterer  is  electrically  small  or  intermediate  in  size,  only  a 

few  of  the  eigenvalues  X  are  of  small  magnitude,  even  though  there  may 

n 

be  many  ports.  Hence,  in  a  modal  solution,  only  a  few  of  the  modeB  may 

be  required  for  good  accuracy.  A  mode  having  X  -  0  is  said  to  be  in 

n 

resonance.  In  many  cases  only  that  mode  which  is  at  or  near  resonance 
contributes  significantly  to  the  scattered  field. 

The  dual  short-circuit  admittance  formulation  for  the  scattered 
field  is  [7] 


-  1>S  + 


(7) 


Here  [Yc]  and  [Y  ]  are  the  port  admittance  matrices  of  the  scatterer  and 

"*8C 

load,  respectively,  I  is  a  column  matrix  of  the  short  circuit  port 


currents,  P*  Is  s  row  nstrlx  of  the  fields  rsdistsd  by  unit  voltage  sources 
st  the  ports,  end  P*'  is  the  field  scattered  when  all  ports  are  short  cir¬ 
cuited.  For  a  nodal  analysis,  the  basis  is  changed  to  one  which  diagonalizes 
[Yg  +  Yl ] .  In  particular,  we  choose  as  a  basis  the  eigenvectors  Vr  of 


(81V,  -  u  [G]$ 
n  n  n 


(8) 


where  [B]  and  [G]  are  the  Hermltian  parts  of  [Yg  +  Y^].  The  advantages  of 
this  choice  are  dual  to  those  listed  (a)  to  (d)  above.  Finally,  when  the 
elgenvoltages  are  chosen  as  a  basis,  (7)  reduces  to 


E 


8 


N  Vt,c 

y  -s _ 

ail  l+J"n 


(9) 


Both  (6)  and  (9)  are  exact  formulations  of  the  problem,  and  either  may  be 
used  for  analysis.  However,  our  synthesis  procedures  usually  involve 
approximations,  in  which  case  (6)  and  (9)  may  lead  to  different  results. 


III.  MODAL  SYNTHESIS 


In  this  section  the  general  concept  of  modal  synthesis  is  discussed. 
Any  real  port  current  (or  voltage)  can  be  made  an  eigencurrent  (or  eigen- 
voltage)  of  the  loaded  scatterer  by  choosing  the  proper  load.  The  procedure 
is  called  modal  resonance  if  the  eigenvalue  of  the  synthesized  mode  is  zero. 

In  general,  an  eigencurrent  must  satisfy  (2)  where  [X]  -  [Xg  +  X^J  and 
[R]  ■  [Rg  +  RjJ.  Rearranging  this  equation,  we  have 

[Xj^  -  ARL]t  -  -  [Xg  -  XRg ]t  (10) 

Both  [Xc]  and  [Rc]  are  determined  by  the  scatterer  geometry,  and  hence  are 
fixed.  If  we  specify  a  real  t  and  a  real  1,  then  (10)  becomes  a  set  of 

and  [R^.  In  general 


linear  equations  for  determining  the  elements  of  [X^ 
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there  arc  aore  unknowns  chan  equations,  and  hanca  cha  aolutlon  la  not  unlqua. 
Wa  can  chooaa  cha  load  network  co  ba  loss-free,  in  which  caaa  (10)  raducaa 
co 

(Xjl  -  -  (Xs  -  XRs)t  (11) 

Tha  aoludon  la  aelll  noe  unlqua  In  ganaral.  Ua  can  alvaya  chooaa  cha  load 
natvork  Co  ba  diagonal,  chac  la 


[XJ  -  (dlag 

In  chla  caaa  cha  diagonal  alaaenca  ara  found  from  (11)  aa 


(12) 


X 


1 


-  '"s'1’! 


(13) 


Hara  (  )  danocaa  cha  1-Ch  component  of  cha  c-ncloaed  column  matrix. 
Finally,  If  wa  wlah  co  raaonaca  Che  currenc  1,  wa  chooaa  X  ■  0  and  (13) 
raducaa  Co 


xi  ■  -  f  (IV 


(14) 


A  raal  pore  currenc  may  alao  ba  reaonaced  by  nondlagonal  load  matrices 

lzL]. 

Tha  above  discussion  applies  In  Che  dual  sense  to  pore  modal 
voltages.  Analogous  Co  (10),  the  basic  equation  to  satisfy  la 


[Bl  -  UGJV  -  -  (B$  -  wGsJV  (IS) 

Again  Cha  scatterer  geometry  determines  both  [Gc]  and  [B_].  If  we  specify 
a  raal  V  and  real  u,  Chen  (IS)  becomes  a  set  of  linear  equations  for  de¬ 
termining  Che  elements  of  [G.  ]  and 

L 

than  aquations,  and  the  solution  la  not  unique.  Specializations  of  (IS) 
analogous  to  (11)  through  (14)  can  be  made  according  to  duality.  In  par¬ 
ticular,  for  resonance  by  a  loss-free  diagonal  load  matrix  [Y^]  ■  j[diag  B^j, 
we  have  dual  to  (14) 


[B^].  Again  there  are  more  unknowns 


b 


Bi  -  -  ^  ([BglV^  (16) 

A  real  port  voltage  may  also  be  tesonated  by  nondiagonal  load  matrices 

iV- 


IV.  PATTERN  SYNTHESIS 


In  this  section  we  consider  synthesis  procedures  for  determining  the 
real  port  current  t  (or  port  voltage  V)  whose  radiation  field  pattern  approxi¬ 
mates  a  desired  field  pattern.  The  method  developed  for  general  surfaces  [4] 
can  be  used  for  N-port  scatterers  with  only  minor  changes  in  the  theory.  How¬ 
ever,  this  synthesis  procedure  requires  the  desired  field  pattern  to  be  speci¬ 
fied  in  both  magnitude  and  phase.  For  most  applications  the  phase  of  the 
field  pattern  Is  unimportant,  and  better  approximations  can  be  obtained  by 
specifying  only  the  magnitude  of  the  field  pattern.  This  is  especially  true 
for  our  problem  because  of  the  constraint  that  the  port  current  be  real.  We 
will  develop  the  theory  in  terms  of  N-port  systems,  but  it  also  applies  to 
surfaces,  such  as  considered  in  reference  [4]. 

The  problem  Is  one  of  mixed  antenna  synthesis,  as  defined  by 
Bakhrakh  and  Troytskiy  [8].  The  method  of  solution  is  similar  to  that 
used  by  Choni  [9].  We  desire  to  obtain  the  real  port  current  I  whose 
radiation  field  E  approximates  in  magnitude  a  real  pattern  .F  on  the  radi¬ 
ation  sphere.  Let  Em  denote  a  component  of  E  at  point  m  on  the  radiation 
sphere,  and  F®  the  corresponding  component  of  F  at  the  same  point.  Given 
M  points,  we  define  the  error  of  synthesis  to  be 


e 


(17) 


Next,  let  {1^}  be  a  basis  of  real  port  currents,  and  express  the  port 
current  as 


-*■ 

I 


N 


n-1 


at 

n  n 


(18) 
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where  a  are  real  constants  to  be  determined, 
n 


The  radiation  field  of  I  Is 


N 

E  -  T  a  E 
**  ,  n«n 

n-1 

where  E^  Is  the  radiation  field  produced  by  the  corresponding  1^. 
of  (19),  the  error  (17)  becomes 


(19) 

In  terms 


e 


M 


l 

m*l 


N 

l 

n-1 


„m 
a  ti 
n  n 


(20) 


We  wish  to  determine  the  real  a  which  minimize  e.  It  is  sometimes  con- 

n 

venlent  to  use  fewer  basis  functions  than  ports,  In  which  case  the  N  of 
(18)  to  (20)  Is  the  number  of  basis  functions  instead  of  the  number  of 
ports. 


To  circumvent  the  troublesome  inner  magnitude  operation  In  (20), 
we  first  consider  the  more  general  function 


e(a,^) 


M 

l 

m-1 


N 

l 

n-1 


_m 
u  E 
n  n 


-I* 


ja, 


m 


(21) 


This  Is  the  error  function  used  when  the  desired  field  pattern  is  specified 
both  in  magnitude  F  and  In  phase  8.  Hence,  for  fixed,  the  for  mini¬ 
mum  c  are  given  by  the  analysis  of  reference  [4].  When  the  a  are  fixed 

n 

the  minimum  e  is  obviously  obtained  when  both  quantities  within  the  magni¬ 
tude  signs  of  (21)  are  In  phase  for  each  m,  or 


N 

8  -  angle  of  l  a  E°  (22) 

m  ,  n  n 

n-1 

Because  (21)  is  more  general  than  (20),  Its  minimum  is  less  than  or  equal 
to  that  of  (20).  But  under  condition  (22),  the  e  of  (21)  is  equal  to  that  of 
(20).  Therefore  (20)  and  (21)  have  the  same  minimum. 
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An  iterative  procedure  for  minimizing  (21)  proceeds  as  follows: 

1.  Assume  starting  values  for  * ^2 *  *  *  * *®m* 

2.  Keep  the  8  fixed  and  calculate  the  a  which  minimize 

r  m  n 

e  according  to  the  theory  of  [4]. 

3.  Keep  the  a  fixed  and  calculate  the  8  which  minimize 

n  m 

c  according  to  (22) . 

4.  Go  to  step  2. 

This  procedure  eventually  converges  because  steps  2  and  3  cannot  lncrsase  c. 
While  the  procedure  obtains  absolute  minima  in  a  space  and  in  t  space,  it 
does  not  necessarily  obtain  the  absolute  minimum  in  the  catenated  space 
(at£).  The  procedure  in  general  convergss  to  a  stationary  point,  usually 
a  local  minimum,  which  may  or  may  not  be  the  global  minimum. 

We  now  describe  steps  2  and  3  in  more  detail.  When  ^  is  constant, 
c (a,3)  is  quadratic  in  o  and  thus  has  only  one  stationary  point.  This  is 
the  absolute  minimum  located  at  (4] 


a  -  (Re(?i'*E)]“1tRs(lf*Fe;)S)]  (23) 

where  the  matrices  are 


E 


(24) 


[F“e  "1 


M*1 


(25) 


When  a  is  constant,  we  adjust  the  phase  t  according  to  (22).  An  alterna- 
tlve  way  of  expressing  this  is 


(Ba)m 


(26) 


where  (Ea)  denotes  the  m-th  component  of  the  column  vector  Ea.  Note 
m 

t'let,  since  only  the  exponential  (26)  Is  used  In  the  Iteration  pro¬ 
cedure,  we  need  never  calculate  the  8  . 

m 
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The  synthesis  procedure  In  terms  of  port  voltage  V  Is  the  same  except 
for  a  change  In  definition  of  various  quantities.  Instead  of  (18),  we  ex¬ 
press  the  port  voltage  as 

N 

V  -  l  ct  V  (27) 

r-1  n  n 


where  (V  }  Is  a  real  basis  and  a  are  real  constants  to  be  determined, 
n  n 

Equation  (19)  remains  the  same,  except  E  Is  the  radiation  field  produced 
by  the  corresponding  V  .  The  rest  of  the  theory  applies  unchanged. 


V.  EXAMPLES  -  PATTERN  SYNTHESIS 

The  theory  la  applicable  to  a  large  variety  of  problems,  and  even  when 
the  problem  is  specified  there  still  remains  the  choice  of  the  formulation 
(open-circuit  or  short-circuit)  and  of  the  basis  functions.  For  the  geometry 
wa  choose  the  wire  triangle  with  two  cross  wires  shown  In  Fig.  1.  This  Is 
the  same  object  uaed  In  the  previous  report  on  modal  analysis  of  N-port  sys¬ 
tems  [7],  The  tip  angle  is  30*,  the  parameter  a  Is  one  wavelength  at  a  fre¬ 
quency  fo,  and  the  wire  diameter  is  a/100.  The  four  points  at  which  the  wires 
cross  the  z  axis  are  Input  ports,  labeled  (1),  (2),  (3),  and  (4).  All  com¬ 
putations  are  made  using  38  triangle  functions  in  a  Galerkln  solution.  (The 
computer  input  data  Is  given  in  Part  two,  Section  II  of  reference  [7].) 

The  magnitude  of  the  pattern  chosen  for  synthesis  Is 


coe  6  In  the  x-0  plane, 

(28) 

cos  0|  In  the  y-0  plane. 

(29) 

Twelve  points  are  chosen  for  the  least-squares  pattern  synthesis  as  follows: 

Foi  F  in  the  x-0  plane , 

♦ 

0-0,  30* 
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60*,  90’,  120*,  150*,  180* 


(30) 
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Fig.  1.  Wirt  triangle  with  crott  wlret,  tip  angle  -  30*, 
a  ■  one  vavelength,  wire  dlaaeter  •  a/100. 


11 


and  for  FQ  in  the  y-0  plane, 

6-30°,  60°,  90°,  120°,  150°  (31) 

The  pattern  functions  F.  at  0  and  180°  in  the  y»0  plane  are  the  same  as 
at  0  and  180°  In  the  x=0  plane,  and  hence  are  net  included  in  (31).  The 
phase  angles  of  the  field  were  chosen  to  be  zero  at  all  points  in  the  first 
iteration  of  the  pattern  synthesis  procedure. 

For  the  first  example  we  use  port  currents  as  the  basis  for  pattern 
synthesis.  To  illustrate  convergence,  first  one,  then  two,  then  three  and 
finally  four  port  currents  are  used.  They  are  added  in  the  order  in  which 
the  ports  are  labeled  in  Fig.  1.  This  order  is  arbitrary,  and  any  other 
order  could  be  chosen  if  desired.  Figure  2  illustrates  the  results  of  this 
synthesis  procedure.  The  solid  curves  show  the  magnitude  of  the  synthesized 
E^  and  EQ  in  the  x-0  and  y-0  planes,  respectively,  and  the  crosses  show  the 
desired  field  magnitudes.  The  port  currents  for  the  final  synthesized  pat¬ 
tern  are 


Ix  -  0.3428 

12  -  0.2923  (32) 

13  - -0.1048 

I.  - -0.0013 
4 

Note  that  there  is  little  change  betveen  the  eyntheelzed  patterns  of  Fig.  2c 
and  Fig.  2d.  This  le  reflected  in  the  relative  emallneite  of  1^  above.  Hence, 
for  thle  particular  ayntheela  approach,  one  can  do  almost  ae  veil  by  ualng 
only  three  porta  and  lervlng  the  fourth  one  open  circuited.  It  la  to  be 
emphasized  that  the  ayntheela  procedure  la  an  optimum  seeking  one,  usually 
having  several  local  optima.  Hanes,  starting  from  different  initial  phae<ie 
for  the  field,  or  ualng  e  different  order  of  Iteration,  ve  may  arrive  at  u 
different  final  synthesis. 


(a)  one  port  current 


(b)  two  port  currents 


’ 


(c)  three  port  currents  (d)  four  port  currents 

Fig.  2.  Field  Bsgnltude  pattern  synthesis  using  reel  port  currents 
es  a  basle.  Crossee  denote  the  desired  pattern.  Curves 
labeled  F  .  denote  | E .  |  In  the  x»0  plane.  Curves  labeled 
E  g  denote  j EQ |  In  the  y«0  plane. 
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Once  ch^  desired  reel  pore  cur  rente  ere  obeelned,  they  cen  be  resoneted 
by  the  procedu  e  of  Section  III.  This  hee  been  done  for  eech  of  the  synthe¬ 
sized  petterne  of  Pig.  2.  Since  port  currents  were  used  ee  e  bests,  thie 
mesne  thet  in  Fig.  2s  e  reective  loed  wee  pieced  only  et  port  (1),  in  Fig.  2b 
et  ports  (1)  end  (2),  in  Fig.  2c  et  ports  (1),  (2),  end  (3),  end  in  Fig.  2d 
st  sll  four  ports.  It  is  understood  thet  ports  st  which  no  reective  loed  is 
pieced  sre  open  circuited.  The  resultent  ecetterlng  petterne  ere  shown  in 
Fig.  3.  The  plots  sre  in  terms  of  *£/X,  since  this  is  the  field  msgnltude 
which  wss  used  in  the  synthesis  procedure.  In  eech  '.see  the  incident  weve 
is  sn  x-polsrlzed  plene  weve  propsgstlng  in  the  t-dlrectlon  (incident  on  the 
30*  tip  sngle).  Note  thet  the  scattering  petterne  of  figs.  3s,  3b,  3c,  end 
3d  ere  slmllsr  to  the  synthesized  petterne  of  Figs.  2s,  2b,  2c,  sod  2d,  re¬ 
spectively.  The  degree  to  which  the  ecetterlng  petterne  ere  the  eeae  ss  the 

synthesized  petterne  depends  upon  the  smallness  of  the  "beckground  ecetterlng," 
oc 

thet  is,  of  the  term  in  (6).  In  the  present  problea,  for  the  open-circuited 
scstterer  of  Fig.  1, 

o/X2  -  0.0177  (33) 

In  contrsst  o/X2  is  of  the  order  unity  for  the  loaded  scstterer.  Hence,  our 

OC 

assuaption  that  the  E^  term  is  smell  is  fully  justified  in  this  csss.  The 
rssctlve  losds  used  in  the  final  loaded  scstterer  (Fig.  3d)  are 

Xl  -  -  100.3 
X  -  -  546.7 

1  (34) 

X3  -  -  1448. 

X^  -  47,370. 

Note  that  X^  is  large,  Indicating  that  port  (4)  is  almost  sn  open  circuit. 

It  is  an  open  circuit  in  the  scstterer  of  Fig.  3c,  but  tha  other  losds  are 
then  slightly  different  from  those  of  (34).  This  again  indicates  that  port 
(4)  is  relatively  unimportant  for  this  particular  synthesis  procedure. 


(c)  three  pore  current* 


(d)  four  pore  current* 


3.  Blataelc  scattering  pattern*  (So/\)  for  the  wire  object  of 
Fit*  1  loaded  co  resonate  the  port  currents  which  synthesise 
the  patterns  of  Fig*  2.  Curve*  labeled  x#  denote  t  polarise 
Clou  In  the  x»0  plane.  Curve*  labeled  y6  denote  6  polarisa¬ 
tion  In  the  y»0  plane. 
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In  most  practical  applications  one  seekB  scatterers  having  desirable 
scattering  characteristics  over  some  reasonably  broad  frequency  band.  The 
synthesis  procedure  Is  basically  a  single  frequency  one,  with  no  consider¬ 
ation  given  as  yet  to  variation  with  frequency.  Ae  an  indication  of  fre¬ 
quency  sensitivity,  we  have  calculated  the  backscattering  cross  section 

per  wavelength  squared  (o/X2)  vs.  frequency  over  the  range  f  -  0.8f  to 

o 

1.2fQ,  where  f  is  the  design  frequency.  Figure  4  shows  curves  for  the 
same  loaded  scatterers  as  considered  in  Fig.  3.  In  each  case  the  load  is 
considered  to  be  an  inductance  if  X^  is  positive  at  fQ,  or  a  capacitance 
if  is  negative  at  f^.  Note  that  even  the  singly-loaded  case,  Fig.  4a, 
has  considerable  variation  of  o/X2  vs.  frequency.  Hence,  we  conclude  that 
our  synthesized  loaded  scatterers  are  relatively  frequency  sensitive.  (Any 
non-smoothness  of  the  frequency  plots  are  due  to  relatively  large  incre¬ 
ments  in  frequency  used  over  parts  of  the  frequency  range.) 

The  same  synthesis  problem  is  next  treated  using  the  modal  currents 
of  the  open-circuit  formulation  as  the  basis.  Qualitatively  speaking,  we 
expect  the  modal  currents  corresponding  to  the  smaller  eigenvalues  to  have 
broader  band  characteristics.  Hence,  if  only  currents  with  small  eigen¬ 
values  are  us«d,  structures  with  broader  band  characteristics  should  result. 
To  test  this  hypothesis,  we  added  the  modal  currents  in  the  order  of  In¬ 
creasing  magnitude  of  eigenvalues.  These  mode  eigenvalues  and  elgencurrents 
(normalized  to  maximum  value  unity)  are  tabulated  in  the  following  table. 

Table  1.  Eigenvalue  *  and  slgencurrents  *  for  the 

d  n 

wire  object  o'.  Fig.  1. 


D 

X 

n 

Port  (1) 

Port  (2) 

Port  (3) 

Port  (4) 

n 

-0.1552 

-0.1338 

0.4326 

C.8419 

1.0000 

B 

-10.12 

-0.6078 

1.0000 

0.8054 

-0.6458 

B 

Eg | 

1.0000 

-0.5374 

0.0137 

0.0740 

B 

ESI 

0.5441 

1.0000 

-0.2839 

0.0778 

(c)  three  aode  currents  (d)  four  mode  currents 


Fig.  S.  Field  sugnltude  pattern  synthesis  using  reel  node  currents 

es  e  heels.  Crosses  denote  the  desired  psttern.  Curves 

lebeled  E  denote  | E .  |  In  the  x»0  plene.  Curves  lebeled 

E  .  denote  |E.|  In  the  y»0  plene. 
y9  0 
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(c)  thr**  code  current*  (d)  f<*«ir  aod*  current* 


Fig.  6.  Blatatlc  scattering  pattern*  (*^o/l)  for  tb*  vlr*  object  of 

Fig.  1  loaded  to  r**onat*  th*  port  current*  which  aynthealt* 
the  p*tt*m*  of  Fig.  5.  Curve*  labeled  x+  denote  #  polarise 
tlon  In  th*  x*0  plane.  Curve*  labeled  jr0  denote  6  polarlr*- 
tion  In  th*  y»0  plan*. 
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We  used  first  ons,  than  two,  than  threa,  and  than  all  four  aigancurrenta  in 
tha  synthaaia  procadura.  Tha  results  are  shown  in  Figs.  Sa  to  Sd,  respec¬ 
tively.  Again  tha  solid  curves  are  the  magnitude  of  tha  synthesized  E.  in 

v 

tha  x-0  plana  and  Eg  in  tha  y-0  plana,  while  tha  crosses  era  the  desired 
field  magnitudes.  Tha  synthesized  pattern  using  only  one  baala  function, 

Fig.  Sa,  la  Just  tha  pattern  of  tha  single  mods  currant  1^.  Thla  mode  has 
components  of  currant  at  all  ports,  as  shown  in  tha  first  row  of  Table  1. 

In  fact,  each  synthesized  pattern,  Figs.  Sa,  Sb,  5c,  5d,  now  results  from 
currants  at  all  ports.  Tha  final  synthesized  pattern,  Fig.  Sd,  is  the  same 
as  tha  final  synthesized  pattern  using  port  currents,  Fig.  2d.  However,  the 
final  patterns  in  tha  two  cases  do  not  necessarily  have  to  be  the  same,  since 
tha  synthesis  method  is  Iterative  and  may  arrive  at  a  different  local  opti¬ 
mum  point. 

Again,  once  wa  have  a  desired  real  current  It  can  be  resonated  by 
tbs  procedure  of  Section  III.  Tha  currants  producing  each  of  the  synthe¬ 
sized  patterns  of  Fig.  S  have  bean  resonated  in  this  way.  This  now  involves 
placing  a  reactive  load  at  each  port  in  each  caaa.  The  scattering  patterns 
for  tha  raactlvely  loaded  acattarara  art  shown  in  Fig.  6.  As  before,  these 
plots  are  of  *G/A,  which  la  field  magnitude  as  uaad  in  tha  synthesis  pro¬ 
cadura.  Tha  incident  wavs  in  each  case  is  x-polarlzed  and  z-propagatlng. 
Again,  note  that  tha  scattering  patterns  of  Flga.  6a,  6b,  6c,  and  6d  are 
similar  to  tha  synthesized  patterns  of  Figa.  Sa,  Sb,  Sc,  and  Sd,  respec¬ 
tively.  Tha  reactive  loada  uaad  in  tha  final  case  are  essentially  the  same 
as  (34),  since  era  a  . lvad  at  tha  same  optimum  point.  Note  that  each  scat¬ 
ters  of  Fig.  6  involves  loada  at  all  porta.  This  is  in  contrast  to  the 
preceding  case  where  only  one  load  was  uaad  in  Fig.  3a,  two  in  Fig.  3b, 
and  three  in  Fig.  3c. 

Va  next  wish  to  demonstrate  that  broader  band  behavior  results 
whan  a  partial  modal  basis  la  uaad  for  tha  current.  Figure  7  shows 
graphs  of  backscatterlng  a/X  va.  frequency  for  tha  loaded  scatterers 
of  Fig.  6.  Again  tha  load  is  conaldarad  to  be  an  Inductor  if  X^  la 
positive  at  fQ,  or  a  capacitor  if  X^  la  negative  at  tQ.  Note  that  the 


(a)  one  node  current 


(b)  two  mode  currents 


0.001 


H - t- 


1.0  0.9  1.0  1.1  1.2 

l/So 

(c)  three  mode  currents 


0.0001  + 


H - h- 


0.8  0.9  1.0 


H - 1- 


1.1  1.2 

f/f 


(d)  four  mode  currents 


7.  Backscattering  (o/A2)  vs.  frequency  (f/fQ)  the  same  loaded 
scatterers  as  Fig.  6. 
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single  mode  case,  Fig.  7a,  gives  extremely  broadband  scattering.  The  two 
mode  case,  Fig.  7b,  is  still  relatively  broadband,  the  three  mode  case, 

Fig.  7c,  less  so,  and  the  four  mode  case,  Fig.  7d,  quite  narrow  band.  The 
four  mode  case  is,  in  fact,  identical  to  Fig.  4d,  the  final  design  being 
the  same. 

Other  bases  for  the  current  can  be  used  in  the  synthesis  procedure. 

For  example,  the  mode  currents  for  some  particular  loaded  structure  can 
be  used  [7].  For  the  broadest  band  behavior,  we  can  use  the  Q  mode  cur¬ 
rents,  defined  in  a  manner  analogous  to  that  for  complete  conducting  bodies 
[4].  In  particular,  the  Q  mode  currents  I  are  defined  by 

w[X']t  -  Q[R]t  (35) 

n  n  n 

where  [X']  is  the  angular  frequency  derivative  of  the  port  matrix  [X] . 

For  computation,  we  used  a  finite  difference  approximation  to  [X'].  The 

Q  mode  currents  are  then  obtained  from  (35)  by  the  same  method  used  to 

solve  (2)  for  the  ordinary  mode  currents.  It  turns  out  that,  for  the 

particular  scatterer  of  Fig.  1,  the  Q  mode  currents  did  not  differ  greatly 

from  the  ordinary  mode  currents.  The  eigenvalues  are  different,  but  the 

ordering  according  to  magnitude  of  Q  remained  the  same  as  for  the  X  eigen- 

n  n 

values.  Hence,  computations  using  Q  mode  currents  were  similar  to  those  of 
Figs.  5  to  7.  Plots  of  the  resulting  synthesized  patterns  and  scattering 
are  therefore  not  shown  explicitly.  Suffice  it  to  say  that  the  scatterers 
of  Figs.  6  and  7  are,  for  most  practical  purposes,  the  broadest  band  ob¬ 
tainable. 

A  set  of  curves  similar  to  Figs.  2  to  7  has  been  prepared  using  the 
short-circuit  formulation  and  port  voltages.  The  results  are  qualitatively 
similar,  except  that  the  correspondence  of  the  scattering  patterns  for  loaded 
scatterers  to  the  corresponding  synthesized  patterns  is  not  quite  as  good. 
This  is  because  the  "background  scattering"  for  the  short-circuited  scatterer, 

8C  oq 

that  is,  the  Eg'  term  in  (9),  is  not  as  small  as  the  analogous  Eq  term  in 
(6).  In  fact,  for  the  short-circuited  scatterer  of  Fig.  1, 
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o/A2  -  0.2011  (36) 

which  Is  significantly  larger  chan  (33)  for  the  open-circuited  scatterer. 
For  reference,  the  curves  for  the  short-circuit  formulation  are  given  in 
Appendix  A.  These  correspond  to  Figs.  2  to  7  for  the  open-circuit  formu¬ 
lation.  Finally,  so  that  the  reader  can  compare  the  behavior  of  loaded 
scatterers  to  the  same  scatterer  unloaded,  plots  of  the  open-circuited  and 
short-circuited  backseat  ter  ing  o/A2  vs.  frequency  are  given  in  Fig.  8. 


VI.  MAXIMUM  RADAR  CROSS  SECTION 


Another  problem  of  Interest  is  that  of  maximizing  the  radar  cross 
section  of  a  loaded  scatterer.  In  terms  of  the  open-circuit  modal  analysis, 
the  bistatic  radar  cross  section  is  given  by  [7,  Eqs.  (69)  and  (74)] 


o 


2  2 
UTJJ* 


4ir 


„  r  t 
N  a  a 
oc  y  _n — n 

0  +  L.  1+j  A 
n-1  J  n 


(37) 


Here  F-  is  the  term  due  to  scattering  by  the  open-circuited  scatterer, 

r  t 

and  a  and  a  are  mode  excitation  coefficients  for  the  n-th  mode  due  to 
n  n 

excitation  from  the  receiver  and  transmitter,  respectively.  It  is  desired 
to  reactlvely  load  the  scatterer  to  resonate  a  port  current  (A^  -  0)  such 
that  its  contribution  is  the  only  significant  term  in  (37) .  Then 


o  # 


2 


(38) 


We  here  consider  only  the  case  of  backscattering,  whence  a'  -  a*.  The 
coefficient  a*  is  given  by  [7,  Eq.  (55)]  with  the  port  current  normalized 
according  to  I*[R]I  ■  1.  A  formula  for  |a*|^  which  is  Insensitive  to  the 
amplitude  of  1  is  therefore 


li^l2 

T*[r]I 


p 


(39) 
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where  V  Is  Che  open-circuit  pore  voltage  produced  by  a  unit  plane  wave 
from  the  receiver.  The  problem  is  now  to  find  the  real  port  current  t 
which  maximizes  (39).  Once  the  desired  1  is  found,  it  can  be  resonated 
by  the  method  of  Section  III. 

The  parameter  p  is  also  proportional  to  the  power  gain  g  of  the 
current  f,  the  explicit  relationship  being  [10] 

8  -  ^  P  (40) 


Hence,  maximization  of  radar  cross  section  is  equivalent  to  the  maximiza¬ 
tion  of  gain,  but  with  the  restriction  that  I  is  real.  In  the  more  general 
case  of  bistatic  scattering,  (38)  can  be  written  as 


(41) 


r  t 

where  g  and  g  are  the  power  gains  in  the  directions  of  the  receiver 
and  transmitter,  respectively.  Thus,  maximization  of  bistatlc  radar  cross 
section  is  equivalent  to  maximizing  the  product  of  two  gains,  again  with 
the  restriction  that  I  is  real.  The  relationship  (41)  was  previously  de¬ 
rived  for  one-port  loaded  scatterers  in  reference  [2], 

The  optimization  of  Rayleigh  quotients  of  the  form  (39)  when  I  is 
restricted  to  the  real  field  is  considered  in  reference  [4].  In  this  case 
the  maximum  p  is 


Pr  -  Re(V0C)[R]_1{Re(V°C)  +  c  Im(V°C)}  (42) 

and  the  associated  port  current  is 

Ir  -  [R]_1{Re(V°c)  +  c  Im(V°C) }  (43) 

In  (42)  and  (43),  c  is  chosen  to  give  the  larger  pf  from  the  two  numbers 


c  ■  -  a  +  /a*  +  1 


(44) 
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where 


Re[V°C] fRl~1Re(V°C)  -  Im(V°C) [Rl-1  Im(V°C) 
2  Re(VOC){R]"1  Im($0C) 


(45) 


Even  though  I  must  be  real  (or  equlphase)  for  (39)  to  be  valid,  It  is  of 
interest  to  compare  the  above  result  with  the  maximum  p  that  could  be  ob¬ 
tained  if  I  were  complex.  This  problem  is  treated  in  reference  [10],  the 
maximum  p  being 


„o  -  v‘>c[Rr:l7o<:*  cm) 

and  the  associated  port  current  being 

tQ  "  [R]_1V°C*  (47) 

It  has  been  shown  that  the  maximum  p  for  real  port  currents  is  at  least 
one-half  that  for  complex  port  currents,  and  the  two  cases-  are  equal  if 
V°C  is  equlphasal  [4], 

The  analysis  can  alternatively  be  carried  out  in  terms  of  short- 
circuit  parameters  and  port  voltages.  In  this  case  the  appropriate 
starting  formula  is  [7,  Eqs.  (74)  and  (76)].  It  is  now  desired  to  sus- 
ceptively  load  the  scatterer  to  resonate  (p^  ■  0)  a  port  voltage  such 
that  it  is  the  only  significant  term  in  the  bistatlc  radar  cross  section 
formula.  Analogous  to  (38),  we  then  have 


2, ,2 


a 


4ir 
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(48) 


r  t 
where  8^  and  8^ 

mode.  Again  we 

The  coefficient 

llzed  according 


are  the  mode  excitation  coefficients  for  the  resonant 

T  t 

consider  only  the  case  of  backseat  ter  ing,  whence  8^  -  8^. 

8,  is  given  by  [7,  Eq.  (63)],  with  the  port  voltage  norma- 

1  ^  -*• 

to  V*[G]V  -  1.  Therefore,  analogous  to  (39),  we  have 
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(49) 
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where  I  is  the  short-circuit  port  current  produced  by  a  unit  plane  wave 
from  the  receiver.  The  problem  is  now  to  find  the  real  port  voltage  V  which 
maximizes  (49).  The  maximization  procedure  is  identical  to  that  used  in  the 
preceding  case,  except  for  an  interchance  of  symbols.  Hence,  if  we  replace 
V°c  by  iSC,  I  by  V,  and  [R]  by  [G],  equations  (42)  to  (47)  remain  valid  in 
the  present  case.  Once  the  desired  real  V  is  obtained  it  can  be  resonated 
by  the  methods  of  Section  III.  Because  of  ♦■he  approximations  involved,  the 
solution  in  terms  of  V  will  not  give  exactly  the  same  radar  cross  section 
as  the  solution  in  terms  of  I,  but  the  two  cross  sections  are  usually  very 
close. 


We  can  think  of  the  preceding  solution  as  the  maximization  of  p  when 
I  (or  V)  is  expressed  in  terms  of  the  port  ba*is.  It  ir.  of  interest  to  con¬ 
sider  the  solution  when  an  arbitrary  basis  is  used,  or  even  when  an  incom¬ 
plete  basis  is  used.  To  accomplish  this,  let  {1^}  denote  a  real  basis  and 
let 
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n**l 


n  n 


(50) 


where  the  a  are  real  constants  to  be  determined.  Then 
n 
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1voc  -  I  a  I  V°c  -  oV 
L  n  n 


oc 


n-1 


(51) 


where  a  is  a  row  vector  of  the  a  and  0  is  the  new  voltage  vector 

n 


voc  -  [I  ]*°c 

n 

V  <v 

Here  [InJ  is  the  matrix  with  rows  equal  to  1^.  Similarly, 

N  N 

I*[R]t  -  7  l  a  ?  [R]t  a 

*  '  L  L  m  m  n  1 


(52) 


m«l  n-1 


n  n 


a[R]a 


(53) 
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where  [R]  le  the  new  resletence  matrix  with  elemente 


R  -  I  [R]I 
mn  m  n 


(54) 


In  tenu  of  the  new  basis,  the  p  of  (39)  becomes 


a  [Rjo 


(55) 


which  is  again  the  same  form  as  (39).  Hence,  the  optimization  of  p  pro¬ 
ceeds  as  before.  The  theory  also  remains  valid  if  the  basis  {1^}  is 
incomplete.  In  this  case  the  N  of  the  above  equations  becomes  the  number 
of  basis  functions  instead  of  the  number  of  ports.  A  change  of  basis  for 
the  solution  in  terms  of  port  voltages  is  accomplished  in  a  similar  way. 


VII.  EXAMPLES  -  MAXIMUM  CROSS  SECTION 

Representative  computations  of  maximum  rcdar  cross  section  have  been 
made  for  the  4-port  wire  object  of  Fig.  1.  Again  the  choices  of  the  formu¬ 
lation  and  of  the  basis  are  arbitrary.  For  the  first  example  we  uBe  the 
open-circuit  formulation  and  port  currents  for  the  basis.  Gain  is  optimized 
using  first  one  port,  then  two,  then  three,  and  finally  four  real  port  cur¬ 
rents.  The  ports  are  added  in  the  order  in  which  they  are  numbered  in  Fig.  1, 
which  is  arbitrary.  Once  the  real  port  currents  are  found,  they  are  reso¬ 
nated  by  the  method  of  Section  III.  This  results  in  a  reactively  loaded 
scattcrer  which  gives  approximately  the  maximum  backscattering  radar  cross 
section  in  the  direction  that  gain  was  maximized.  Table  2  summarizes  these 
results.  The  first  column  shows  the  ports  used  in  the  gain  maximization, 
which  are  also  those  ports  loaded  when  the  scatterer  is  resonated.  The 
second  column  gives  the  maximum  x-polarized  gain  in  the  -z  direction  for 
real  currents  at  the  excited  ports.  The  third  column  gives  the  x-polarized 
backscattering  cross  section  for  the  resonated  scatterer  in  the  -z  direction. 
In  this  case,  those  ports  not  used  in  the  optimization  procedure  are  open 
circuited. 
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Fig.  9.  Bistatic  scattering  patterns  (o/l2)  for  the  wire  object  of 
Fig.  1  loaded  for  maximum  gain  in  the  -z  direction.  Curves 
labeled  x$  denctu  $  polarization  in  the  x-0  plane.  Curves 
labeled  y8  denote  0  polarization  in  the  y-0  plane. 
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fable  2.  Maximum  gain  and  backscattering  cross  section 
when  port  currents  are  used  as  a  basis. 


Ports  Used 

Maximum  Gain 

Cross  Section 

1 

1.017 

0.3293  X2 

1,2 

1.185 

0.3322  X2 

1,2,3 

3.972 

5.005  X2 

1,2, 3,4 

4.033 

5.160  X2 

The  bistatic  radar  cross  section  patterns  (o/ A2)  for  the  loaded 
scatterers  of  Table  2  are  shown  in  Fig.  9.  The  curves  labeled  x<f> 
are  for  the  ^-polarized  scattered  field  in  the  x-0  plane,  and  those 
labeled  yO  are  for  the  6-polarized  scattered  field  in  the  y«0  plane. 

When  only  port  1  is  loaded,  Fig.  9a,  the  maximum  gain  pattern  is  the  same 
as  the  field  synthesis  pattern,  Fig.  3a.  (Note,  however,  that  Fig.  3a  is 
a  plot  of  r/a/\,  whereas  Fig.  9a  is  of  o/X2.)  When  ports  1  and  2  are  loaded, 
the  scattering  pattern  is  still  not  very  directive  in  the  desired  -z  direc¬ 
tion.  When  ports  1,2,  and  3  are  loaded  (Fig.  9c),  the  scattering  pattern 
becomes  highly  directive  in  the  -z  direction.  Finally,  there  is  little 
Improvement  over  the  three  port  case  when  all  four  ports  are  loaded,  Fig. 

9d.  The  reactive  loads  for  this  final  case  are  as  follows 

Xx  -  -693.2 

X2  -  -361.7 

(56) 

X3  -  -  933.9 

X,  -  -4949. 

4 

Note  that  is  large  compared  to  the  other  loads,  which  we  would  expect 
since  there  is  little  change  when  it  is  open  circuited.  Plots  of  the  gain 
patterns  for  maximum  gain,  corresponding  to  Figs.  2  and  5  in  the  pattern 


(a)  one  port  current 


ia. 


o/X2 

l.  «" 


0.1 

0.01 

0.001 

O.OQOl 
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0.S  0.9  1.0  t.l  1.2 


,/fo 

(c)  three  port  currents 

Fig.  10.  Backseat ter ing  (c/X2)  vs. 

loaded  scatterers  as  Fig. 


(b)  two  port  currents 


(d)  four  port  currents 

frequency  (f/fQ)  for  the  same 
9. 
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synthesis  problem,  are  not  shown.  They  are  similar  in  form  to  the  correspond¬ 
ing  scattering  patterns  of  Fig.  9. 

The  behavior  of  the  synthesized  scatterers  over  a  frequency  bandwidth 
is  also  important,  just  as  in  the  prttern  synthesis  case.  Figure  10  shows 
the  variation  of  o/A2  vs.  frequency  over  the  range  f  ■  0.8fQ  to  l*2fQ,  where 
f  is  the  design  frequency,  for  the  four  cases  of  Table  2.  Again  a  load  is 
considered  to  be  an  inductor  if  is  positive  at  fQ,  or  a  capacitor  if 
is  negative  at  fQ.  Note  that  each  scatterer  is  relatively  narrow  band,  with 
the  frequency  sensitivity  increasing  as  the  number  of  loads  is  Increased. 

Note  also  that  the  one-port  optimization  result,  Fig.  10a,  is  the  same  as 
the  one-port  synthesis  result,  Fig.  4a.  This  is  because  only  one  basis 
function  is  used  in  each  case. 

The  gain  optimization  problem  is  next  treated  using  the  modal  currents 
of  the  unloaded  body  as  the  basis.  These  modal  currents  are  again  those 
listed  in  Table  1.  The  currents  are  added  to  the  optimization  procedure  in 
the  order  of  Increasing  |*nl»  just  as  in  the  previous  pattern  synthesis 
problem.  Again  the  x-poiarlzed  gain  in  the  -z  direction  is  maximized,  and 
Table  3  summarizes  the  results.  The  first  column  shows  the  number  of  modes 
used,  the  second  column  gives  the  maximum  gain  obtained,  and  the  third 
column  lists  the  back-scattering  cross  section  obtained  when  the  scatterer 
is  resonated.  Note  that  the  fourth  rows  of  Tables  2  and  3  are  identical 
since  a  complete  basis  was  used  in  both  cases.  Also,  it  should  be  empha¬ 
sized  that  for  each  case  of  Table  3  reactive  loads  are  placed  at  all  ports 
in  contrast  to  the  case  of  Table  2  where  loads  are  placed  only  at  those 
ports  used  in  the  basis. 


(c)  three  mode  currents 


(d)  four  mode  currents 


Fig.  11.  Bistatic  scattering  patterns  (o/A2)  for  the  wire  object  of 
Fig.  1  loaded  for  maximum  gain  in  the  -z  direction.  Curves 
labeled  x$  denote  $  polarization  in  the  x«0  plane.  Curves 
labeled  yd  denote  6  polarization  in  the  y-0  plane. 


Table  3.  Maximum  gain  and  backscatttrlng  cross  section 
when  mode  currents  srs  used  as  a  basis 
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Modes  Ussd 

Maximum  Cain 

Cross  Section 

One 

0.609 

0.1934  A2 

Two 

1.044 

0.4757  A2 

Three 

1.554 

0.5659  A2 

Four 

4.033 

5.160  A2 

The  blstatlc  radar  cross  ssctlon  patterns  (0/ A2)  are  shown  in  Fig. 

11  for  the  loaded  scatterera  listed  in  Table  3.  The  curves  are  labeled 
in  the  same  manner  as  those  of  Fig.  9.  Again  the  one-mode  case,  Fig.  11a 
(0/A2  plotted)  Is  the  same  as  the  one-mode  synthesis  case,  Fig.  6a  (/j/A 
plotted).  It  is  also,  of  course,  just  the  bistatic  scattering  pattern  ob¬ 
tained  by  directly  resonating  the  dominant  mode  current.  The  two-mode 
case,  Fig.  lib,  has  a  significantly  higher  cross  section  in  the  desired  minus 
z-dlrectlon  than  does  the  one-mode  case.  The  Improvement  in  backscatter- 
lng  cross  section  in  the  three-mode  case,  Fig.  11c,  is  small  over  the  two¬ 
mode  case.  Finally,  the  four-mode  case,  Fig.  lid,  is  identical  to  the 
four-port  case  of  Fig.  9d,  since  a  complete  current  basis  is  used  in  both 
cases.  Also,  the  final  loading  reactances  must  be  the  same  for  these  two 
cases.  Hence,  the  loads  (56)  are  also  those  used  in  the  final  four-mode 
synthesis  procedure. 

The  broadband  behavior  of  the  loaded  scatterers  listed  in  Table  3 
is  summarized  by  the  curves  of  Fig.  12.  The  one-mode  case  is,  of  course, 
the  same  in  the  pattern  synthesis  problem  as  in  the  gain  optimization 
problem,  hence  Fig.  12a  is  identical  to  Fig.  7a.  The  two-mode  case,  Fig. 

12b,  is  more  narrow  band  than  the  one-mode  case ,  but  the  Improvement  in 
backscattering  cross  section  is  significant.  The  three-mode  case,  Fig. 

12c,  is  still  more  narrow  band,  and  the  improvement  in  backscattering 
cross  section  over  the  two-mode  case  is  only  slight.  Finally,  the  four- 


/ 

I 


same 
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mode  case,  Fig.  12d,  is  very  frequency  sensitive,  and,  in  fact,  is  Identical 
to  the  four-port  case  of  Fig.  lOd. 

We  can  use  still  other  bases  for  gain  optimization  if  desired.  Con¬ 
ceptually,  the  broadest  band  behavior  should  result  when  the  Q-mode  currents, 
defined  by  (35),  are  used  as  a  basis  in  order  of  increasing  lQnl«  However, 
as  noted  earlier,  for  this  problem  there  are  no  significant  differences  be¬ 
tween  the  ordinary  mode  currents  used  above  and  the  Q  mode  currents  defined 
by  (35).  Hence,  results  using  the  Q-mode  currents  for  gain  optimization  were 
very  similar  to  those  of  Table  3  and  Figs.  11  and  12,  and  are  not  shown 
separately. 

Finally,  all  of  the  examples  of  this  section  can  be  run  for  the  dual 
formulation  in  terms  of  short  circuit  parameters  and  port  voltages.  The 
results  are  qualitatively  similar  to  those  obtained  using  the  open-circuit 
formulation.  However,  as  noted  in  the  syr  -heals  problem,  the  "background 
scattering"  for  the  short-circuit  f ormulat:  on  is  larger  than  for  the  open- 
circuit  formulation.  Thus,  neglecting  this  term  is  a  somewhat  poorer 
approximation  for  the  particular  problem  being  treated.  For  reference,  the 
results  for  optimum  gain  and  backseat ter ing  cross  section  using  the  short- 
circuit  formulation  are  summarized  in  Appendix  B. 


VIII.  DISCUSSION 

Given  an  N-port  conducting  body,  any  set  of  real  port  currents  can 
be  resonated  by  reactive  loads  according  to  the  concepts  of  Section  III. 

If  the  resonated  current  provides  the  major  contribution  to  the  scattered 
field,  then  the  reactlvely  loaded  object  will  have  a  scattering  pattern 
approximately  equal  to  the  radiation  pattern  of  the  same  object  excited 
by  the  given  set  of  port  currents.  Various  synthesis  and  optimization 
procedures  can  be  used  to  obtain  real  port  currents  whose  fields  have  de¬ 
sirable  radiation  characteristics.  In  this  report  we  give  two  such  procedures, 
one  which  synthesizes  a  desired  field  pattern,  and  the  other  which  optimizes 
power  gain. 


I 

I 
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The  synthesis  procedure  involves  a  specification  of  the  radiation 
field  in  magnitude  only,  and  then  uses  an  iterative  procedure  to  de  ermine 
the  phase  which  minimizes  a  certain  error  parameter.  The  problem  is  basically 
nonlinear,  and  by  using  different  starting  points  we  sometimes  arrive  at  dif¬ 
ferent  ending  points  (local  minima) .  Some  examples  of  this  for  the  wire 
object  of  Fig.  1  are  given  in  Appendix  C. 

In  the  optimization  procedure,  the  power  gain  is  maximized  subject  to 
the  condition  that  the  port  currents  are  real.  The  solution  in  this  case  is 
unique.  When  resonated  by  reactive  loads,  this  maximizes  the  backscatterlng 
radar  cross  section  under  the  assumption  that  only  the  resonated  mode  con¬ 
tributes  to  the  scattering.  Other  modes  are  usually  negligible  when  the 

scatterer  is  electrically  small  or  of  intermediate  size,  but  not  if  it  is 

oc 

electrically  large.  There  is  also  a  background  term  present,  the  term  ^ 

8C 

in  the  open-circuit  formulation  (6)  or  the  term  in  the  short-circuit 

OC 

formulation  (7).  In  the  first  case,  E.  can  usually  be  made  small  if  enough 

8C 

open-circuits  are  placed  on  the  body.  The  term  E^  is  probably  harder  to 
make  small  in  general,  since  this  involves  placing  short  circuits  on  the  body. 

Sometimes  Eq  may  be  small  and  Eq  large,  in  which  case  the  open-circuit 

formulation  should  be  used.  For  other  bodies  the  reverse  may  be  true,  in 

which  case  the  short-circuit  formulation  is  preferable.  J 

f 

K 


\ 
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APPENDIX  A 

PATTERN  SYNTHESIS,  SHORT-CIRCUIT  FORMULATION 

The  examples  of  this  Appendix  are  dual  to  those  of  Section  V.  The 
wire  object  considered  is  the  same,  Fig.  1,  as  for  the  short-circuit 
formulation.  For  pattern  synthesis,  the  magnitude  pattern  is  again  given 
by  (28)  and  (29).  The  points  chosen  for  the  least-squares  solution  are 
again  given  by  (30)  and  (31) .  The  first  example  uses  port  voltages  as  a 
basis,  with  convergence  illustrated  by  choosing  first  one,  then  two,  then 
three,  and  finally  four  ports.  Figure  13  shows  the  results  presented  in 
the  same  manner  as  in  Fig.  2  for  the  dual  case.  The  port  voltages  for  the 
final  synthesized  pattern  are 


Vx  -  161.4 

V-  -  214.4 

(57) 

V3  —117.3 

V,  -  27.3 
4 

Note  that  the  synthesized  patterns  in  each  case  are  considerably  different 
from  the  correi ponding  patterns  of  Fig.  2.  This  is  because  the  radiation 
fields  resulting  from  voltages  applied  to  ports  with  all  other  ports  short 
circuited  are  different  from  those  resulting  from  currents  applied  to  ports 
with  all  other  ports  open  circuited.  We  should  expect  the  two  formulations 
to  give  different  synthesized  patterns  because  the  space  of  pattern  functions 
for  real  port  voltages  is  different  from  that  for  real  port  currents. 

Next  the  computed  real  port  voltages  are  resonated  by  the  procedure 
of  Section  III.  Figure  14  shows  the  resulting  bistatic  scattering  patterns 
for  when  the  four  cases  of  Fig.  13  are  resonated.  For  Fig.  14a,  we  have 
a  susceptance  at  port  1  and  the  other  three  ports  are  short  circuited.  For 
Fig.  14b,  susceptances  are  at  ports  1  and  2,  and  short  circuits  at  ports  3 


(a)  one  port  voltage 


(b)  two  port  voltages 
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(d)  four  port  voltages 

Fig.  14.  Bistatic  scattering  patterns  (*^o/X)  for  the  wire  object 
of  Fig.  1  loaded  to  resonate  the  port  voltages  which 
synthesize  the  patterns  of  Fig.  13.  Curves  labeled  x$ 
denote  $  polarization  in  the  x»0  plane.  Curves  labeled 
y@  denote  6  polarization  in  the  y*0  plane. 


•f 


(a)  one  port  voltage 


(b)  two  port  voltages 
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and  4.  For  Fig.  14c,  susceptances  are  at  the  first  three  ports,  and  a  short 
circuit  at  the  fourth.  Finally,  for  Fig.  I4d  susceptances  are  at  all  four 
ports.  The  loads  used  in  the  final  case  are 

Bx  -  0.000927 


B2  -  0.002884 
B,,  -  0.001355 


(58) 


3.  -  0.001246 
4 


They  are,  of  course,  different  for  the  intermediate  cases. 

The  frequency  sensitivity  of  the  synthesized  scatterers  is  illustrated 
by  Fig.  15.  In  each  case,  when  a  load  B^  is  positive  at  the  design  frequency 
f  it  is  considered  to  be  a  capacitor,  or  when  B^  is  negative  at  f  it  is 
considered  to  be  an  inductor.  Note  that  the  one-port  case,  Fig.  15a,  happens 
to  be  broadband,  although  we  did  not  design  it  to  be  so.  The  other  three 
cases,  Figs.  15b  to  15d,  are  all  relatively  narrow  band. 

The  second  example  of  this  Appendix  uses  the  modal  voltages  of  the 
body  as  the  basis  for  the  optimization  procedure.  These  mode  eigenvalues 
and  eigenvoltages  (normalized  to  maximum  value  unity)  are  given  in  the  follow¬ 
ing  table. 


Table  4.  Eigenvalues  and  eigenvoltages 
for  the  wire  object  of  Fig.  1. 


n 

yn 

Port  (1) 

Port  (2) 

Port  (3) 

Port  (4) 

1 

0.1552 

-0.0205 

0.1156 

0.6425 

1.0000 

2 

10.12 

0.1904 

0.2524 

1.0000 

-0.9256 

3 

50.54 

1.0000 

-0.2887 

0.7884 

-0.4050 

4 

-816.2 

0.5313 

1.0000 

-0.6749 

0.2066 

Note  that  the  vn  of  Table  4  are  equal  to  -A  of  Table  1,  as  proved  in  [7]. 


(c)  three  mode  voltages 


(d)  four  mode  voltages 


Fig.  16.  Field  magnitude  pattern  synthesis  using  real  mode  voltages  as 

a  basis.  Crosses  denote  the  desired  pattern.  Curves  labeled 

E  ,  denote  I E .  I  in  the  x»0  plane.  Curves  labeled  E  .  denote 
x$  $  y® 

|  E0  |  in  the  y-0  plane . 


(a)  one  mode  voltages 


(b)  two  mode  voltages 


Z 


(c)  three  mode  voltages 


17.  Blstatlc  scattering  patterns  (♦''o/A)  for  the  wire  object  of 
Fig.  1  loaded  to  resonate  the  port  voltages  which  synthesize 
the  patterns  of  Fig.  16.  Curves  labeled  denote  $  polari¬ 
zation  in  the  x*0  plane.  Curves  labeled  y6  denote  6  polari¬ 
zation  in  the  y*0  plane. 
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Again  tha  synthesis  procedure  was  run  using  first  one,  then  two,  then  three, 
and  finally  all  four  eigencurrents  in  the  basis.  The  results  are  shown  in 
Figs.  16a  to  d,  respectively.  The  synthesized  pattern  using  only  one  basis 
function,  Fig.  16a,  is  just  the  field  pattern  of  the  first  mode  alone.  The 
synthesized  pattern  using  all  mode  functions,  Fig.  16d,  is  the  same  as  that 
using  all  port  bases,  Fig.  13d.  This  final  correspondence  is  somewhat  for¬ 
tuitous,  since  the  two  procedures  could  arrive  at  different  points  of  local 
optima. 

The  real  port  voltages  are  next  resonated  by  susceptive  loads  accord¬ 
ing  to  the  procedure  of  Section  III.  Figure  17  shows  the  bistatic  scatter¬ 
ing  patterns  when  each  case  of  Fig.  16  is  resonated.  This  time  there  are 
susceptances  at  each  port  for  all  cases,  regardless  of  the  number  of  basis 
functions.  The  susceptances  used  in  the .  f  inq.l...c,a.*e,»..She  scatterer  for 
Fig.  17d,  are  the  same  as  those  given  by  (58).  This  is  because  we  hav& 
arrived  at  the  same  solution  as  we  did  using  the  port  voltage  basis. 

Finally,  the  frequency  sensitivity  of  the  synthesized  scatterers  is 
illustrated  by  Fig.  18.  Again  a  load  is  considered  to  be  a  capacitor 
if  positive  at  f^,  or  an  Inductor  if  negative  at  fQ.  The  frequency  sensi¬ 
tivity  of  the  loaded  scatterers  is  seen  to  increase  as  the  number  of  modes 
used  in  the  synthesis  procedure  is  increased  (in  the  order  of  Increasing 
|wq|).  It  la  interesting  to  note  that  the  one-port  loaded  scatterer  of 
Fig.  15a  is  considerably  more  broadband  then  the  one-mode  synthesized 
scatterer  of  Fig.  18a.  In  fact,  it  is  not  greatly  different  from  the 
scattering  when  all  ports  are  short  circuited,  Fig.  8b.  In  this  case  the 
"background  scattering"  is  of  the  same  order  of  magnitude  as  that  from 
the  synthesized  scatterer,  violating  our  assumption  that  It  should  be 
small. 


z 


(a)  one  port  voltage 


Z 


(c)  three  port  voltagea  (d)  four  port  voltages 

Fig.  19.  Blatatic  scattering  patterns  (o/X2)  for  the  wire  object  of 
Fig.  1  loaded  for  maximum  gain  in  the  -z  direction.  Curves 
labeled  x$  denote  +  polarization  in  the  x-0  plane.  Curves 
labeled  yd  denote  8  polarization  in  the  y°0  plane. 
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APPENDIX  B 

MAXIMUM  CROSS  SECTION,  SHORT-CIRCUIT  FORMULATION 


The  examples  of  this  Appendix  are  dual  to  those  of  Section  VII.  The 
wire  object,  Fig.  1,  is  the  same,  but  the  short-circuit  formulation  is  used 
for  gain  maximization.  The  theory  is  dual  to  that  summarized  in  Section  VI. 
The  first  example  uses  port  voltages  as  a  basis.  Gain  is  optimized  using 
first  one  port,  then  two,  then  three,  and  finally  all  four  ports,  added  in 
the  order  in  which  they  are  numbered  in  Fig.  1.  Once  the  real  port  voltages 
for  maximum  gain  are  found,  they  are  resonated  by  the  method  of  Section  III. 
The  result  is  a  loaded  scatterer  having  approximately  maximum  (for  the  ports 
used)  backscattering.  Table  5  summarizes  these  results.  The  first  column 
shows  the  ports  used,  the  second  column  gives  the  maximum  gain  for  real  cur¬ 
rents,  and  the  third  column  gives  the  backscattering  cross  section  of  the 
resonated  scatterer. 

Table  5.  Maximum  gain  and  backscattering  cross  sections 
when  port  voltages  are  used  as  a  basis. 


Ports  Used 

Maximum  Gain 

Cross  Section 

1 

2.418 

0.3054  A2  | 

1.2 

2.458 

1.139  A2 

1.2,3 

2.753 

1.546  A2 

1.2, 3,4 

3.730 

4.046  A2 

The  bistatic  radar  cross  section  patterns  (o/A2)  for  the  loaded 

•catterers  of  Table  It  are  shown  in  Fig.  19.  The  curves  labeled  o  ./A2 

x<p 

are  for  the  ^-polarized  scattered  field  in  the  x-0  plane,  and  those 
labeled  o^/A2  are  for  the  6-polarized  scattered  field  in  the  y-0  plane. 
Figure  19a  is  for  the  case  of  only  one  port  loaded,  Fig.  19b  two  ports, 
Fig.  19c  three  ports,  and  Fig.  19d  four  ports.  The  susceptive  loads  for 
the  final  case,  Fig.  19d,  are 
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(a)  one  mode  voltage 


(c)  three  mode  voltages  (d)  four  mode  voltages 

Fig.  21.  Bistatic  scattering  patterns  (o/X2)  for  the  wire  object  of 
Fig.  1  loaded  for  maximum  gain  in  the  -z  direction.  Curves 
labeled  x$  denote  <j>  polarization  in  the  x-0  plane.  Curves 
labeled  y6  denote  6  polarization  in  the  y*0  plane. 
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(c)  three  mode  voltages  (d)  four  mode  voltages 


Fig.  22.  Backscattering  (o/X2)  vs.  frequency  <f /f Q)  for  the  same  loaded 
scatterers  as  Fig.  21. 
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Bx  -  0.001215 

B2  -  0.002870 

B3  -  0.001036 

B.  -  0.000483 
4 


(59) 


The  loads  for  Che  cases  of  Figs.  19a  Co  19c  are,  of  course,  differenc. 

The  frequency  senslclvlcy  of  Che  loaded  scaccerers  of  Fig.  19  is 
llluscraced  by  Fig.  20.  A  load  B^  is  considered  Co  be  a  capacitor  if 
positive  at  fQ,  or  an  Inductor  if  negative  at  fQ.  Note  Chat  Che  one- 
load  optimization  case,  Fig.  20a,  is  Che  same  as  Che  one-load  synthesis 
case,  Fig.  15a,  since  Che  sane  basis  voltage  is  resonated  in  both  cases. 

The  second  example  of  Che  Appendix  uses  the  modal  voltages  as  a 
basis  for  the  optimization  proce^'T; .  These  mode  eigenvalues  and  elgen- 
voltagea  are  those  listed  in  Table  4.  The  gain  is  optimized  using  first 
one  mode,  then  tvo,  then  three,  and  then  four,  added  in  the  order  of  in¬ 
creasing  magnitude  of  eigenvalues.  The  optimum  gain  voltages  are  next 
resonated  by  the  concepts  of  Section  III.  The  results  are  sumr^rlzed 
lr>  Table  6.  It  should  be  noted  that  we  now  have  loads  at  all  ports  of 
the  loaded  scatterer  for  each  case,  in  contrast  to  the  preceding  case, 
Table  5,  where  load*  were  only  placed  at  those  ports  used  in  the  basis, 
other  ports  being  short  circuited. 

Table  6.  Maximum  gain  and  backscatterlng  cross  sections 
when  mode  voltages  are  used  as  a  basis. 


Modes  Used 

Maximum  gain 

Cross  Section 

One 

0.6C9 

0.01990  X2 

Two 

0.616 

0.08658  X2 

Three 

2.003 

0.4682  X2 

Four 

3.730 

4.044  X2 
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The  blstatic  radar  cross  section  patterns  (a/ X2)  for  the  loaded 
scatterers  of  Table  6  are  shown  in  Fig.  21.  These  are  labeled  in  the 
same  manner  as  are  those  of  Fig.  19.  The  susceptive  loads  for  the  final 
case,  Fig.  21d,  are  the  same  as  those  for  the  case  of  Fig.  19d,  given  by 
equations  (59).  This  is  because  the  basis  was  complete  in  both  cases. 

Finally,  the  frequency  sensitivity  of  the  loaded  scatterers  of 
Fig.  21  is  illustrated  by  Fig.  22.  Again  a  load  is  considered  to  be  a 
capacitor  if  positive  at  fQ,  or  an  Inductor  if  negative  at  fQ.  Finally, 
note  that  the  one-mod*  optimization  case,  Fig.  22a,  ia  the  same  as  the  one¬ 
mode  synthesis  case.  Fig.  18a,  sines  the  same  mode  was  resonated  in  both 
cases. 
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APPENDIX  C 

EXAMPLES  OF  DIFFERENT  OPTIMA  IN  PATTERN  SYNTHESIS 


The  field  magnitude  pattern  synthesis  procedure  of  Section  IV  is  a 
nonlinear  one,  and  may  arrive  at  different  optimum  points  if  different 
starting  points  are  used.  Figure  23  gives  some  examples  of  this  phenomena. 
For  Fig.  23a,  the  initial  phase  angles  of  the  field  points  were  chosen  to 
be  0  and  180°,  alternating  between  adjacent  points.  The  corresponding  case 
for  which  all  initial  phase  angles  were  chosen  zero  is  shown  in  Fig.  5c. 

Note  that  the  two  final  patterns  are  different,  although  both  are  reasonable 
approximations  to  the  desired  pattern.  A  second  example  is  shown  in  Fig. 

23b.  In  this  esse  the  initial  phase  angles  of  the  field  points  were  chosen 
to  be  0,  90*,  160*,  270*,  etc.  The  corresponding  case  for  which  all  initial 
angles  were  chosen  zero  is  shown  In  Fig.  5d.  An  example  using  a  voltage 
basis  is  shown  in  Fig.  23c.  In  this  case  the  initial  angles  of  the  field 
points  were  chosen  to  be  0  and  180*,  alternating  between  adjacent  points. 

The  corresponding  case  for  which  all  initial  angles  were  chosen  zero  is  shown 
in  Fig.  16b.  A  final  example  is  shown  in  Fig.  23d.  Once  again  the  initial 
angles  of  the  field  points  were  chosen  to  be  0  and  180*,  alternating  between 
adjecent  point'.  The  corresponding  case  for  which  all  initial  angles  were 
chosen  zerr  is  shown  in  Fig. 


13c. 


(c)  two  mode  voltages,  (d)  three  port  voltages, 

compare  with  Fig.  16b.  compare  with  Fig.  13c. 

Fig.  23.  Different  optimum  points  for  pattern  synthesis  due  to  different 
choices  of  initial  phase  angles. 
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PART  TWO 

COMPUTER  PROGRAMS 


I.  INTRODUCTION 


The  programs  used  to  compute  the  examples  of  this  report  are  described  and 
listed  in  this  part.  Each  program  is  accompanied  by  an  explanation  of  the 
input  data,  a  verbal  flow  chart,  and  sample  input  and  output  data.  In  general, 
the  input  data  for  the  program  of  a  given  section  depends  upon  the  output  of 
programs  of  previous  sections  or  of  reference  [7].  If  each  program  is  run  with 
the  input  data  listed  in  this  report,  the  input  data  for  any  one  of  these  programs 
can  be  verified  in  terms  of  the  output  of  programs  previously  run.  The  Calcomp 
Plotter  is  used  only  in  Section  VIII. 

II.  PATTERN  SYNTHESIS 


The  program  (pattern  synthesis  program)  described  in  this  section  requires 
some  data  computed  by  the  port  parameter  program  on  pages  60-69  of  [7].  However, 
the  impedance  matrix  program  on  pages  45-52  of  [7]  and  the  excitation  vector 
program  on  pages  53-59  of  [7]  must  ^recede  the  port  parameter  program.  The 
punched  card  data  for  the  impedance  matrix  program  is  exactly  the  same  as  in 
[7].  The  punched  card  data  for  the  excitation  vector  program  is  altered  slightly 
to  obtain  the  second  and  third  polarizations  (see  (103)  of  [7])  Instead  of  just 
the  third.  On  the  first  data  card,  NPAT  is  changed  from  1  to  2  and  on  the  third 
data  card,  NPA(l)-3  is  replaced  by  NPA(1)«2  and  NPA(2)»3.  Except  for  NPAT  and 
NPA,  the  resulting  printed  output  is  exactly  the  same  as  in  [7],  but  mors  data 
is  stored  on  record  2  of  data  set  6.  The  port  parameter  program  is  run  next 
with  NT  ■  145  replaced  by  NT  -  290  on  the  first  data  card.  Except  for  NT,  the 
resulting  printed  output  is  exactly  the  same  as  in  [7]. 

In  the  pattern  synthesis  program,  the  activity  on  date  sets  1  (punched 
card  input)  and  6  (direct  access  input  and  output)  is  as  follows. 
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READ(1,4)  NF 
4  FORMAT (201 3) 

DO  62  MF-1,  NF 

READ(1,4)  N,M,NT,N6,N8,K9,I2,I3,NIV 
READ (1, 7)  (RE(I) ,  >1,M) 

7  FORMAT(7E11.4) 

READ (1,4) (NE(I) ,  I-1,M) 

NZ«N*N8 

READ (1, 71) (FI (I) ,  I-1,NZ) 

71  FORMAT(4E14.7) 

READ (1,4) (NST(I) ,  I-1.N8) 

NTN«I2+(NT-1)*I3+N 
REWIND  6 

SKIP  N6  RECORDS  ON  DATA  SET  6 
READ (6) (PP(I)  ,  I-l.NTN) 

DO  25  J-1.N8 
IF(NST(J).EQ.O)  CO  TO  25 
L2-NST(J) 

DO  83  L-l,  L2 

READ (1, 7) (ANG(I) ,  >1,M) 

83  CONTINUE 
25  CONTINUE 
62  CONTINUE 

Virtually  all  of  tha  main  program  la  contained  In  DO  loop  62.  There  are 
N  ports  at  which  reactive  loads  may  be  placed.  The  negative  of  the  normalized 
electric  field  at  the  I1*1  point  on  the  radiation  sphere  from  one  ampere  at  the 
Jth  port  Is  stored  In  PP(I2+(I-1)*I3+J)  for  >1,2..  NT.  However,  the  specified 
pattern  conalats  of  the  magnitude  RE  of  the  electric  field  at  only  the  NE(1), 
NE(2) , . . .NE(M)ck  of  these  NT  points.  The  basis  of  real  port  currents  appear¬ 
ing  In  (18)  resides  In  FI((J-i)*N+l)  through  FI(J*N)  for  J-1,2,...N8.  The  Iter¬ 
ative  procedure  for  minimizing  (21)  terminates  either  as  soon  as  a  of  (23)  falls 
to  decrease  c  or  after  (23)  has  been  evaluated  N9  times,  whichever  occurs  first. 
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Although  the  program  Is  documented  in  terms  of  the  open  circuit  impedance 

formulation  which  is  obtained  when  NIVj*0,  the  short  circuit  admittance 

formulation  can  be  obtained  by  setting  NIV-0  and  by  changing  12  so  that  the 

normalized  electric  field  from  one  volt  at  one  of  the  ports  is  referenced  in 

PP.  The  variable  NIV  is  necessary  because  equations  (60)  and  (52)  of  [7] 

are  used  to  obtain  Esc  •  II  and  -?°c  •  I£  respectively  from  f8c  and  Voc  which 

are  stored  in  PP  and  appear  in  the  table  on  page  62  of  [7].  The  AHG(I)  read 

inside  nested  DO  loops  25  and  83  is  the  starting  value  (degrees)  of  BT  appearing 

th  ^ 

in  (21)  when  only  the  first  second,  and  J  of  the  N8  basis  functions  are 

used  in  which  case  (18)  becomes 

I-  [  a!  (II-l) 

L.  n  n 
n»l 

For  a  fixed  J,  NST(J)  different  starting  values  of  £  are  triad. 

Klnlausi  allocations  are  given  by 

DIMENSION  LR(H8),  C(K6*ti6) 
in  the  subroutine  LINER  and  in  the  main  pro  gran  by 

COMPLEX  PP(I2-MXT-1)M>4K) 

DIMENSION  FI(K*K8),  KE(M),  NST(NB),  RE(M) , 

El (M*N8 ^ ,  E2(M*N8) ,  E3(M*N8),  E4(M*N8), 

B(N8«(N8+l)/2),  BB(N8*N8) ,  E5(M*N8),  E6(M*N8), 

E7(M*M),  E8(M*M),  E9(K*M),  E10(M*M),  CS(M), 

SN(M),  CS2(M),  SN2(M) ,  PAT(M),  ALP(N8),  CUR(N) , 

ERR(N9+1) 

f|C  gQ 

a-  b-  •  So  ’ 

V°c,  and  F°C  •  u  which  are  stored  in  PP  and  appear  in  the  table  on  page  62 
of  [7].  Nested  DO  loops  20  and  21  store  the  real  and  imaginary  parts  of  E 
of  (24)  by  columns  in  El  and  E2.  Nested  DO  loops  17  and  67  store  the  real 

Aj 

and  imaginary  parts  of  EF  (F  is  an  M  by  M  diagonal  matrix)  by  rows  in  E3  and 
E4.  Nested  DO  loops  26  and  27  store  Re(E*E)  according  to  the  symmetric  mode 
of  storage  in  B. 
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The  index  J  of  DO  loop  25  Indicates  that  only  the  first  J  basis  vectors 
f  appearing  in  (18)  will  be  used.  This  means  that  only  the  first  J  columns 
of  E,  only  the  first  T  rows  of  $F,  and  only  the  upper  left  hand  portion  of 
Re(E*E)  will  be  considered.  Nested  DO  loops  32  and  33  store  Re(E*E)  by  columns 
in  BB.  Statement  44  Inverts  Re(E*E).  Nested  DO  loops  34  and  68  store  the  real 

AS  _  1  A/ 

and  Imaginary  parts  of  (Re(E*E)]  [EF]  by  rows  in  E5  and  ,E6. 

The  column  vector  E  of  electric  fields  produced  by  a  of  (23)  can  be  written 

as 


E  -  [ E7 J [cos  6]  +  (E8]aln  B+j[E9][cos  8]  +  j[E10][sin  8] 


(II-2) 


where 


[E7]  -  (Re(E) ) (Re (E*E) ]_1[Re(EF) ) 
(E8)  -  (Re(E) ] (Ra(E*E) J-1 ( Ia(EF) ) 


(11-3) 


(E9]  -  (la(t))(te(2*E)  )'*(>•&)) 

(E10)  -  (Im(r)]|IU(?»E)f 1  (Im(ET)) 

and  where  (cos  8)  is  a  coltstn  vector  whose  aC^  element  is  coe  6^.  [sin  8] 
le  defined  similarly.  With  (I 1-2) ,  (22)  gives  the  new  l  in  terms  of  the 
old  l. 

(E7  cos  8)  ♦  (E8  sin  8) 


cos  8  • 


S3 


(11-4) 


sin  8  ■ 
n 


(E9  cos  8)  +  (E10  sin  8) 

II  I 

S3 


(II-5) 


where 


S3  -  V CCE7  cos  8)  +  (E8  sin  8)  )2  +  ((E9  coe  8)  +  (E10  sin  6)  )2  (II-6) 

mm  mm 

The  notation  (E7  cos  8)  indicates  the  m^  element  of  the  column  vector 
E7  cos  @e  In  (II-4)  and  (II-5)  the  old  o  is  on  the  right  hand  side  whereas 
tbr»  new  £  is  on  the  left  hand  side.  The  index  L  of  DO  loop  83  indicates  the 
Lth  starting  value  of  £.  Equations  (11-4)  and  C II—  5 )  are  evaluated  for  the 


/ 
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Ith  tiu  in  DO  loop  47.  DO  loop  47  stores  c  of  (21)  in  ERR,  the  magnitude  of 

the  approximate  field  in  PAT,  and  cos  6  and  sin  6  of  (II-4)  and  (11*5)  in 

d  ta 

CS2  and  JN2. 

DO  loop  55  puts  the  final  a  of  (23)  in  ALP.  Depending  upon  whether  NIV*0  or 

not,  DO  loop  58  puts  either  the  port  voltages  or  port  currents  in  CUR.  DO 

loop  80  puts  the  real  and  Imaginary  part  a  of  the  approximate  electric  field  in 
CS  and  SN. 

For  the  sample  input  and  output,  the  scatterer  la  the  wire  triangle  of 
Fig.  5  of  (7).  The  magnitude  F  of  the  desired  electric  field  pattern  is  given 
by 

F  •  u^  cos  0|  in  the  x*0  plane  (11*7) 

F  •  u^  cos  6  in  the  y«0  plane  (1 1-8) 

To  obtain  a  continuous  360*  pattern  in  (11-7),  u^  •  u^  while  6*0  on  the  positive 
t  axis,  6-90*  on  the  negative  y  axis  and  6*360*  on  the  positive  x  axis.  Equation 
(11-8)  is  valid  for  6*0  on  the  positive  x  axis,  6*9u*  on  the  positive  x  axis,  and 
6*360*  on  the  positive  x  axle.  Equation  (11-7)  is  specified  at  6*0*,  30*,  60*. 
90*.  120*,  150*,  and  180*  while  (11-8)  Is  specified  at  6-30*.  60*,  90*.  120*, 
and  150*.  Because  of  syesMtry,  It  would  have  been  redundant  to  specify  (11-7) 
et  210*,  240*.  270*,  300*,  and  330*  or  to  specify  (11-8)  at  0*,  180*,  210*, 

240*,  270*,  300*,  and  330*.  The  open  circuit  impedance  formulation  (N1V/0)  is 
used  with  only  three  basis  vectors,  the  dominant  (X  •  -0.1552,  X  *  -  10.12,  and 
X  *  -  50.54)  mode  currents  for  the  unloaded  wire  triangle.  Seven  significant 
figures  were  obtained  for  t.^eee  node  currents  by  changing  format  statement  60 
in  the  elgencurrent  program  mentioned  on  pages  72  and  73  of  ('). 
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2  1  on  13  I  *  1  ,u 

K?* JS*  I 
K  I  *  J  1  ♦  I 
S*C!K?) 

CIK?)*C!K  1) 

CHIOS 
13  CUNT  I  Nllt 

L«t J)*IR<IR J) 

IR(IHJ)*L«.I 
IHJ-LK.II  I  1  <*  *  h  .  1 A 
R  J 1 « J 1 ♦LL 
9  CONTINUE 
HFTUMN 
END 

COMPLEX  PP»3?34),U? 

01  MENS  I  ON  FI ! 141.NF 120)  ,NST<4) ,RE(?0| ,E1! M0|,E?IH0) ,F3(R0),E4(fi0) 
DIMENSION  R| 10) ,RM 14)  .ESIRO) ,E4!R01 ,E7! 144) , ERI 144)  ,E9( 144) 
DIMENSION  E 101 1 44 ) ,CS I  20  I , SN( 20 1 ,CS?  (  ?0 I , SN2 ( 20  I , PAT ( ?0 ) . ALP(4  ) 
DIMENSION  CUR (4)  ,EMR<  100)  ,ANr.l?0) 

P«*3.  141593/lflO. 

R  E AD  I  1 • 4 )  NE 

4  FORMA! ( ?0| 3  ) 

MR  I  1 E 13,41  I  NF 

4 1  FORMA  1  (  *0  NF  •  / 1 X » l 3 ) 

DO  4?  MF- | ,NF 

RFADI 1,4)  N , M ,N T , N4  »NR, N9 , 1 ? ,  I  3,0  '  V 
MR  I TF  (  3 , S )  N,M,NT , N4,NR,N4,I?,I3,NIV 

5  FORMAT!  *0  N  M  NT  N4  MR  N9  12  13  MIW1X, 813,14) 

REAOIl ,7IIRE ( I ) , I ■ 1 ,M) 

7  FORMA  1 ( 7E 1 1 • 4  I 

MRITEI3.RMRE!  I),I*1,M) 

R  FORMAT!  'ORE'/!  1X.7E11.4)  ) 

READ!  1 ,41  !NE(  I  I, I « 1 , M 1 
MRI  TE  I  3.9HNF!  1  ),  t-l.M) 

7  FORMAT! *ONE'/l IX, 2013) I 
02  »N*NR 

READ!  1  ,71  MFI  !  I  )  .  I  «  1  ,NZ  > 

71  FORMA  T ( 4E 14.71 

MRITEI 3,24) IFI ( I  )  , 1-1  ,NZ  I 
24  FORMAT! 'OFI •/< 1X.4E14.7)  I 
READ! 1 ,41  INST! I  I , 1-1 , NR) 

HR|TE(3.m(NST(l  ) ,  I  ■  1 ,  NM  ) 

19  FORMAT! ' ON ST '/(iX, 20131) 

NTN-I ?♦ ( NT- 1 )  ♦  1 3*N 
REMIND  4 
JEINA)  11.11,12 
12  00  41  J«1,N6 
READ! 4) 

41  CONTINUE 

11  READ! 4 ) 1 PP (  I  l,|=|,NTNI 
MS  NN-N-N 

MR  I TE ( 3  «  4? ) ( PP (  I  1,1*1. NN I 

42  FORMAT! 'OVS'/I 1X.SE14.7)  I 
Jl »NN*  1 

J?»2*NN 

MR  I  TE ! 3, R4 I ( PP ( I I.I-J1.J2) 

R4  FORMAT! '02S'/( 1X.SE14.7) ) 

MR  I  T  E  !  3  »  R  T  ) 

H7  FORMAT! 'CISC, FSC.VOC.FCC' I 
J2»I2-N-1 


DO  ft  6  .1=  1  «  M 

JS  =  J?  +  INHJI-1  1*13 

J3=Jb+l 

.)  4  =  J  b  +  l  3 

WRI1E(3*H8)(PP(I  I  ,  1  =.13  »  J4  I 
8ft  EIIRMAK  1 X , bE 14. 7 ( 

H6  CONTINUE 

EE  =  0. 

Dll  4b  J=1  , M 

EE  =  Eb+RE(.l>*REI  J) 

4b  CONTINUE 

WRITE(3,76)  EE 

7 6  FORMA  1  (  * OSUM  RE**2=* .F14.7) 

ERR ( 1  )  =EE 
Dll  20  .1=1  » NH 
J1=(J-1)*M 
.14=  I  J-l  )*N 
Dll  21  l  =  l,M 
J3=12+(MH I >-l  )*I3 
.I2  =  J1  +  I 
02=0. 

Dll  22  K=  1  ,N 
.14= J  3+K 
Jb  =  J6  +  K 

II2=U2+PP<.I4)*FI(.I5> 

22  CUNTINUE 

E1(J2)=RFAL(U2) 

E2( J?)=AIMAG(U2) 

21  CONTINUE 
20  CONTINUE 

IFINIV.E(J.O)  GO  TO  63 
DO  64  .1=1,. 12 
E 1 1 J ) =-E 1 f  J ) 
r-2C  J)=-E?(  J) 

64  CONTINUE 

63  WR  I  TE  I  3 , 36  )  (  E  1  (  I), 1  =  1, 7) 

36  FORMAK  '  OE  1  •  /  (  1 X  ,  7E  1 1  .4  )) 
WRITE(3,37)(E2( 11,1=1,7) 

37  EIIRMATI  '0F2'/(  IX,7E11.4)) 

J4  =  0 

DO  17  J=1 ,Nfl 
DO  67  I  =  1  ,M 
J4= J4+ 1 

E3(  J4)=E1(.J4)*RE(  I  ) 

E6(J4)=E2(.I4)*REI  I  ) 

67  CONTINUE 
17  CONTINUE 
J3  =  0 

DO  26  .)  =  1  , M ft 
J1  =  ( .1-1  )*M 
DO  27  I=l,J 
J3  =  J3+  1 
ft( J3)=0. 

.16=  (I  —  II 
DO  2ft  K  =  1 ,  M 
J4  =  J 1 +  K 
.lb  =  J6  +  K 

R(  J3)  =  BI  J3)+E1(  J4)*E1(.15)+E2(  J4)*E2(  J5» 
2ft  CONTINUE 
27  CONTI  Nil F 


?B  CONTINUE 

WRITE!3,lfl)!R(.)),J=l,J3> 
lft  FORMAT ( '  OR  1 / ( 1 X , 7E 1 1 .4  )  ) 

00  25  J=  1  *  Nft 
I F ( NS  T ( J ) .EQ.G)  GO  TO  25 
Jft  =  0 

00  32  K  =  1  *  J 
J4  = ( K- 1 )*J 
DO  33  L  =  1  *  K 
Jfl= Jft+1 
J5= J4+L 
JB  =  (  L-l  )*J  +  K 
RB!J5)=B( Jfl) 
ftR ( J6 ) =BR( J5  ) 

33  CONTINOE 
32  CONTINUE 

44  CALL  LINFR(J.RB) 

00  34  K=1,M 
J4=K 

00  6B  I=1,J 
JB=< 1-1 )*J 
E5( J4)=0. 

FBI J4)=0. 

.17  =  K 

00  B9  KK= l *  J 
Jfl  =  JB  +  KK 

F5( J4)=E5( J4I+BBI Jfl)*E3(  J7) 

EB(  J4)=E6<  J4H-BH!  JH)*E4<  J7) 

.17  =  J7+M 
B9  CONTINUE 
J4= J4+M 
6fl  CONTINUE 

34  CONTINUE 

DO  70  K=1,M 
J4=(K-1)#M 
DO  75  1  =  1, M 
J5=J4+I 
E7 ( J5 I =0. 

Eft  I J5 ) =0. 

E9( J5)=0. 

E 10 ( J5 ) =0« 

JB  =  K 
J  7  =  I 

DO  72  KK  =  1 ,  J 

E7(j5)=E7(Jb)+El(J7)*E5(JB) 

Eft  I  J5)=Efl( J5I+E1! J7)*E6!J6) 

E9(J5)=E9<  J5I+E2(.I7)*E5(  JB  I 
E10<J5)=E10<J5)+E2(J7)*EB(JBI 
J7=J7+M 
JB  =  J6  +  M 
72  CONTINUE 
75  CONTINUE 
70  CONTINUE 
WRI TE<  3,40) 

40  FORMAT! '0E7, Eft, E9,E10' I 

WRITE(3,14ME7m  ,  I  =  1 ,7  I  ,  (  Eft  1 1  I,  1  =  1,7)  ,  (E9I  I),  1  =  1,71 
WRITE(3,14ME10(  I  )  ,1=1,71 
14  FORMAT! 1X,7E11.4) 

L2=NST! J) 

00  B3  L  =  1 , L 2 


63 


REAOf  1 ,7  )  (  ANGI  !  )  ,  I  =1  ,M  1 
WRI TE( 3* 1 A  » ( ANGI ! ) , 1  =  1 ,M) 

14  FORMAT  (  'OANG' /(  1  X  ,7H  1  .4)  ) 

nn  is  k  =  l  * m 

AN= ANG ( K ) *PK 
C  S  (  K  )  =  COS  (  AN  ) 

SM(K)=S1N(AN) 
lb  CONTINUE 

on  4  7  1=1 , NR 

12=1+1 

ERR ( I P ) =0 . 
on  50  KK=1,M 
S1=0. 

52  =  o. 

J4  =  KK 

00  44  K  =  1 , M 

S1=S1+E7( J4)*CS(K)+E8( J4)*SN(K) 
S2=S2+E4( J4I*CS(K 1+H10I J4)*SN(K) 

J4= J4+M 
44  CONTINUE 

SH=S1-RF(KK)*CS(KK) 

S4=S2-RE(KK)*SN(KK) 

FRRl  !2)=FRR(  I  2 )  +  S8*SH  + S4 *S4 

53  =  SOR  T( SI #S1 +S2*S2 ) 

PAT ( KK  )  =S3 

I F ( S3  1  51,51,52 

51  CS2(KK)=1. 

SN2IKK 1=0. 

GO  TO  50 

52  CS2IKK )=S1/S3 
SN?(KK)=S?/S3 

50  CONTINUE 

0(1  74  KK=  1  ,M 
CS(KK)=CS2(KK| 

SN(KK)=SN?(KK' 

74  CONTINUE 

IF(ERR(l2).GF.FRR«m  GO  TO  73 
47  CONTINUE 
WRITE (3,53  I 

53  FORMAT! 1  TOO  MANY  ITERATIONS') 

73  WRI1EI 3,13) (ERR< I ) ,1=1 ,12) 

13  FORMA  1 ( 'OERR'/I 1X,5E14.7)  ) 

00  55  1  =  1,. I 
J5= I  1-1  )=M 
ALPI 1 )=0. 

00  54  K=1,M 
J4=J5+K 

ALP  I  1  I  =AI_P(  I  )+E5<  .14  )  *CS(  K  I  +  E4(  .14  )"-SN(  K  ) 

54  CONTINUE 

55  CONTINUE 

WRI TE ( 3,54) ( AlP< I ) , 1=1 , J) 

54  FORMAT (  'OALPHA* /(  IX ,4cl4.7  )  ) 

00  57  |  =  1  ,N 
COR  I  I  1*0. 

J4  ■  I 

00  6H  K  *  1  ,  .1 

CUR  m*r.OR(  I  )«ALPIX  I  »F  I  I  J4I 
J4*.|4+N 
CONTINUE 
57  CONTINUE 
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wri  TFn,6oxr.im(  n  ,  1  =  1  ,ni 

60  FORMAT!  'OPDRT  OIIANT  ITI  F  S  •  /  (  1  X  ,  At  14 . 7  )  > 

WR I T6 I  3,77  > ( PAT (  !  > , I  =  1  , M  ) 

77  FORMAT  (  •  OAPPRflX  I  F  I  '  /  (  1  X  ,  7F  1 1 .4  It 
DO  80  K  =  1 ,M 
CS(K)=PAT(K)*CS2(K) 

sn(ki=pat!ki*sn2(k) 

RO  CONTINUE 

WRITE(3,fll)(CS(K) ,  SN ( K  ) , K= 1 ,M ) 

HI  FORMAT ( 'OARPROX  E ' / ( 1 X , 7E  l  1  .4  )  ) 

83  CONTINUE 
25  CONTINUE 
6?  CONTINUE 
STOP 
END 

THAI  A 
l 

4  12290  2  3  40341  10  1 

0. 1000E  +  01  0.R660E+00  0.5000F+00  O.OOOOF+OO  0.5000E+00  0.8660E+00  0. lOOOE+Ol 
0.R660F+00  0.5000E+00  O.OOOOE+OO  0.5000E+00  0.R660E+00 
1  13  25  37  49  61  73158170182194206 
-0. 133R342E+00  0.4326460F+00  0.84186R4F+00  0. 9999999F+00 
-0.6078 1 72E+00  0.9999999F+00  0 .8054470E+00-0. 4457505E +00 
0. 999999 6E +  00-0. 5 37 368 2E +00  0. 1 367921 E-01  0. 739674DF-01 
0  0  1 

O.OOOOF+OO  0. OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOb+OO  O.OOOOE+OO  O.OOOOE+OO 
O.OOOOE+OO  O.OOOOF+OO  0.0000r+00  O.OOOOF+OO  O.OOOOF+OO 
tS  T(1P 
/* 

// 

PRINTED  (IUTP1/T 

NF 

1 

N  M  NT  N6  NR  N9  17  13  NIV 
4  17790  7  3  60341  10  1 


RE 

0.1000E+01  0.8660E+00  0.5000E+00  O.OOOOE+OO  0.5000E+00  0.8660E+00  0.1000E+01 
0.8660F+00  0.5000F+OO  O.OOOOE+OO  0.5000E+00  0.R660E+OO 

NE 

1  13  25  37  49  61  73158170182194206 


F  1 

-0.1338342E+00  0.4326460E+00  0. H4 1 H684F+00  0 . 9999999F +00 
-0,60781 7 2E +  00  0. 9 9999 9 9 F  + 00  0 .8,  5  •  70F +00-0 .64  57505F  +  00 
0.9999996E+or-i.5373(,M?E+00  0. 1367971E-01  0. 7396740F-01 


N  S  T 

0  0  1 


VS 

0.191 72 3 1 1-03  0.3060562F-02-0.39  32  74  7E-03-0.  3 1 965? 3F -07-0. 4  3704  70F-03 
-0. 24  76 1 57 1 -03-0. 6677280E-03  0 . 4  744 35 1 E -03-0. 393? 752F -03-0 . 3 1 96524E-02 
0.10966541-07  0.3973227E-03  0. 19HM4 75E -07  0. 1 4R2670F-02  0. 721 7893F-0? 
-O.50l6743t-03-0.6370674F-03-0.?57615?F-03  0. I9HH475E-02  0. 144766RF-0? 
(1.17866471-02  0. I05447HF-07  0 . 4 3 78 5 6 1 E -0?  0. 7 3969R 3E -03- 0. 66 7 7 ? 7 7 F -03 
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0.42443356-03  0.22128416-02-0.50181736-03  0.437H553E -02  0 . 23970036-03 
O.5248:(866-0?  0.121321  86-02 

7  S 

0.85532236+01-0.18185476+03  0 . 7 385028 F +00  0 . 1 9852776 +03  0 . 1 1 1 2900E +0? 
-0.37882846+03-0.11204006+02  0 . 203 2 7 1 1 E+03  0.73840086+00  0. 1 98527RE+03 
0. 45  3824RE  +  01  0. 72874] 0F  +  03  0 . 1 5280206+02-0 . 24 580R5F +03-0 . 1 7  1423  1 F  +  0  1 
0. 1319428b+03  0. 11 128516+02-0. 3788 257E+03  0. 1 5280 1 36+02-0. 245808 1 F+03 
0. 77 blb03E +  07-0.4788 18HE +  0  3  0 . 444  884 7F  +  0 )  0 . 44b8  748F  +  03-0 . 1  1 7038 3E +0? 
0.20827056+03-0.12191476+01  0 . 1 3 1 44 2 3F+03  0. 4448344F +0 1  0. 4458248E +03 
0.1 1 17707F+03-0.4774771E+03 


i  sr.,Fsc,vnc,Fnr. 

0. l74144bE-01  0. 3044880E-02-0.7802338E-01-0. 1 1 28838E-0 1 -0. 4 54048 3E-0 1 
-0. 237747 bE-01 -0.487bblbE-0 1-0. 3472430E-01  0.8378842E +00  0 . 1 24748 7E-0 1 
-0.2 1428 18F  +  01  0.8 14 748  IE +00-0. 8  14051  3E+00-0.  1 8 1 8085E +  0 1 -0.  788487  IE +00 
0 .58484576+00  0 . 1088378E  + 02  0 . 5553505E+01  0.3882574E-01  0 . 1407508E-0 1 
0. 18881 14E-01 -0.444b 338E-07-0. 3718 383E-01  0. 437544 1 E -02-0. 8 1 7 7 1 74E-0 1 
0.bb204b4E-03-0.88332bbF-01-0.8454431E-07  0.8380722E+00-0.3501515E+00 
-0.580478  lb +  00  0.  1 575 388E +  01-0. 14  78348E  +  00-0.  1 352030F  +  01  0.  1 8  30723E+01 
-0. 1114444E+01  0.U74518F+02  0 . 8744 750E+00  0 . 7 380058F -0 1  0 . 3047889E-0 1 
0. 4447 345E-07-0. 7028884 E-0 1-0. 1287844E-01  0. 448 7018E-0 1 -0. 7794789E-01 
0.8088 7 58F-01-0.7157417E-01  0.88 5404 7F-01 -0.7877 14 5E +00-0 . 89 37 255E+00 
0. 1 5 304 55F+0 1-0. 8422108b +00-0. 31* 1 404E +00-0. 7 1 1 97R8F+0 1 -0. 1 5 728 3 1 E+0 1 
-0. 4 4 84454 b +  01  0 . 3430484F +  0 1 -0 . 8 558  1  78F  +  0 1  0 . ] 4 3844  5E +00-0 . 4424 1 74E-0 1 
-0. 1444873E-01-0.4051385E-02  0 . 47 30245E-0 1  0. 1 024248E -0 l  0. 4 1 1 4557E-0 1 
0. 1T87571E-01  0.1071B17E+00-0. I 8855B 5E-07-0 . 4 1 2001 0E +00  0 . 58 24828E +00 
-0. 1448554E+01-0. 108 7 748 E +00-0. 2807508E+0 1  0. 142 5888E +00-0. 1501971 E +02 
0. 7355479 E +00-0. 84408 57E +01  0 . 1 8 382 1 OE+O l -0 . 1 322 157E +00-0 . 1 59809 3E+00 
-0.8873806E-03  0 .435285 5F-02-0 . 27 1 4045E -02-0. 2BR9325E-0 1 -0. 9 380578E-02 
-0.7143170E-01-0. 1954884b -01 -0 .9028485E-01  0 . 7 148822E +00  0.88830586+00 
0.15358206+01  0. 74284026+00  0 . 34 7 348 1 E-02  0.2147792E+01  0. 1452803E+00 
0.98878856+01  0 . 5825748b +0 1  0 . 88 14459E +0 1 -0 . 50 1 03 14E -0 1  0 . 1802 144F+00 
-0.18818876-02-0.38334426-02-0. 1998757E-01  0. 4851 8 93E-02-0. 504 1 1 746-01 
0.4234578 6-02-0.85 18740E -01-0. 7 188884F-03  0 . 850444 OE +00-0 .4 349 3 1 2E+00 
-0.3 145  39  lb  +  OO-O.  18 34  94 1 E +0 1-0 . 30524 73E-02  0.  10351 18E  +  01  0. 241 1313E+01 
—  0 . 8  894  44 1 E  +  00  0.9515347E+01-0.2054931E+01  0.  18 54 1 54E +00  0. 1052294E-0 1 
-0.5152740E-02-0.38R2023E-02-0. 1282415E-01  0. 1421350E-01-0. 3345028E-01 
0.2893584E-01-0.4874B45F-01  0 . 285340 l E -01  0 . 308872 7E +00-0 . 8 145 549E+00 
-0. 7022444E +0 1-0.1 295 329E +01-0. 5801 583E +00  0. 143O530E +01-0. 101 7804E+01 
-0.2748825E+01  0 . 809 3R03E +0 1 -0 . 5 198 1 58E +0 1  0 . 1 742824E +00-0 . 1 08502 3E +00 
0. 1431801E-01-0. 1 85845 1 F -02-0. 2404482E-01  0. 3H20483E -02-0. 44344 74F-01 
0.5770248E-02-0.5593947E-01-0.5390318F-02  0.4450209E +00-0 . 38 1 2340F+00 
-0.245430RE+01  0. 2 150214E +01-0. 87352026+00-0.931 2 1536+00-0. 1570337E+01 
0.281 1038E+01  0.1154047E+02-0.31 17700F+01  0.21 1484 7E -02-0. 8987 198F-01 
0.51 2838HE-02  0. 1770929E-03  0. 38859 38E-02  0. 7904314E-02  0. 8R39853E-0? 
0.2830038F-01  0.283331 7F -03  0 .434 1742F-01 -0.47 24085E +00-0 . 28 3088 3F+00 
0.521 534 1F+00  0. 3884882F+0 1  0 . 9 544 1 72F +00-0. 9H74884E-0 1  0. 4000002F+01 
0.457  8  14  21  +01-0.55449  1  IF +  01-0. 1034  1 70F +  02-0 . 2048484F +00  0 .9 7 7802 3E -02 
0.  13955196-01  0.2344323F-02-0. 121 2004F -01-0. RR4 384 3E -02-0.444 519  IE-02 
-0.17458886-01  0.1498838E-01-0. 201487 OE-Ol  0. 1 2 l 7503E +00  0 .45 1 3527E+00 
I). 3534557E+00  0. 7494 1 04E-01  0. 1988842E+00-0.2759507E+01  0.44130846+00 
-0. 44841296+01-0. 1480445E+01  0 .9 3 5 24 54F +0 1  0 . R9 1 4B84E -0 1  0.20513816+00 
0. 88434 256 -02-0. 8 708 7486 -02-0. 2882H57E-02  0. 20885076-01-0. 2947438E -02 
0.2442843E-01-0.1 1747806-01  0 . 222904 1 F-0 l  0. 1 727545E -01 -0. 38388356+00 
-0. 44587 876 +00  0.30784B2E+01-0. 1230785E +01-0.411 8199E +00-0. 5848402t +01 
(1.31789496+01  0.80158H2F+01 -0.74408826+01  0 . 88 30072t -0 1 -0 . 1 708338F+00 
0.  17(170996-02  0. 44087586-03  0 . 1 8 8 2 7 74F -0 1 -0.  1  39484 36-02  0.40078846-01 
-0.32089486-02  0.44738176-01  0 . 5 1 4 1  1 8 56-02-0 . 5224 34  RE *00  0 . 2 5088 3 8 F +00 
0. 1954428F +01  0.248B438F+01  0 . 584 58 l 1 E +00-0. 84 33392F +00  0.38150986+00 


37 

0.43H2909F+01-t).HtH107HF+()l-0.7  33 7049 F ♦ 0 1 -0 . 1  943 339F +00-0 . 7 24735R F -0 1 
SUM  RF**?=  0.5499H?) F+01 


F  1 

-0.  1027F  +02-0.  OPR^  +  OP-O.  17  661-  +01  0.  20 1  1 F+02-0.  5344F  +01-0 . 1  1  39E  +02-0. 7233E+0 1 


F  2 

-0.5133F+01  0.I105E+01  0. 1777E+02-0.2334E+01-0. 13R0E+07  0.2110E+01  0.371RE+01 


H 

0. IR33E+04-0.2344E+03  0.8790E+03  0.  1 003E+02-0. 2425E  +  00  0.7154E+02 


E7,EH,F9,E10 


0.110RF+00  0.953HE-01  O.R992E-02-0.2233F-01  0.2493E-01  0.7490E-01  0.H747E-01 
0.7542E-01  0. 7 104E-01  0. 89 1 1 E-0 1-/< 3R5RE-03  0.9512F-01  0.570RF-01  0.9009E-02 
0.7542E-01-0.  1R75E-01-/<5737E-01  0.133RE-01  0.3R92F-01  0. 1 195E-0 1  -0 . 1 3 1  3E -01 
0.1079E+00  0.H4H1F-01-0. 1R9RE-01  0.4103E-02  0, 72RHE-07-0 . 9791 E -01 -0 . 1 073E+00 


ANT, 

O.OOOOE+OO  0 . OOOOF+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 

o.ooooe+oo  o.ooooe+oo  o.ooooe+oo  o.ooooe+oo  o.ooooe+on 


ERR 

0  .  3499  8  2  1  E  +  0  1 
0.941 003RE  +  00 
0.9124485F+00 
0.9032344F+00 
0.9044H94E+00 
0.90390 1  BE  +00 
0.903ARHRE+00 
O. 903307  3E  +00 
0.9035755E+00 
0.9035323E+00 
0. 903337 3E+00 
U.9035333E+00 

ALPHA 

-0 . 3H  7  1 034E -01 


0.435037 IE  +  01 
0.9309457E+00 
0.9  104964E  +  00 
0.903  7  273F  +  00 
0.9043177E+00 
0.90384 1 2E+00 
0 .9033361 E  +00 
0.90339R7E+00 
0.90357 1 RE+OO 
0.9035314F+00 
0.9035533E+00 
0.9035530E+00 


0.3347580E-01 


0 . 1 01 037 1 E+0 1 
0.923HR21 E+00 
0.9089837E+00 
0.9053047E+00 
0 .904 1793 E+00 
0.9037908E+00 
0.9033439E+00 
0.90359 1 OE+OO 
0.90353R9E+00 
0.9035593E+00 
0.9O35333E+00 
0.9035543E+00 


0.1049RR5E+00 


0. 9739403E+00 
0.9 1  R8403E  +  00 
0.907R334E+00 
0. 904964RE+00 
0.9040370E+00 
0. 9037502E  +  00 
0.90333! 2E+00 
0. 9035R5 1 E+00 
0.9035333E+00 
0. 9035591 E+00 
0.9035559E+00 
0. 9035543E+00 


0.9355495E+00 
0.915 1 R57E+00 
0.90395 1 1E+00 
0. 904  7029E  +  00 
0.9039757E+00 
0. 9037 135E  +00 
0 .9036 1 R 1 E+00 
0. 9035  797E  +00 
0.9035344E+00 
0.90355H IE +00 
0.9035555E+00 
0. 903554RE  +00 


PORI  OUAM  T 1  TIES 

0.89R2205E-01-0.3938939E-0 1-0. 4 1901 2  IE-02-0. 5253191 E-01 


APPROX  | F  | 

0 . 7 53  3E  +  00  0.3R98E  +  00  0.5140E  +  00  0.3970E  +  00  0.4h34E  +  00  0.5433E  +  00  0.5917E+00 
0. H731F+00  0.9036F+00  0.3372E+00  0.7437E+00  0.70R9F+00 


APPROX  E 

0.7414E+00  0.1490E  +  00  0.3393F  +  00-0.1337F  +  00  0. 23R5F -01 -0. 5 1 34E +00-0 . 33 30E +00 
0.1043F+00  0.15R.F+00  0.4173E+00  0.53R4E+00  0.9454E-01  0. 5R03E +00-0. 1 1 l 7E+00 
0.7HR4E+00-0.3H71F+00-0.7R3RF+00-O.R579E+00-0.9377E-01  0.3542E+00  0.32R7E+00 
-O.37C5F+00-0.3147F+00-0.3531F+00 

ONLY  1HREE  MOIJF  CURRENTS  APPEAR  IN  THE  ABOVE  PRINTFO  OUTPUT.  FOR 
REFERENCE,  ALL  THE  MUL)FS  OF  THE  WIRE  TR I  ANCLE  ARE  LISTED  BELOW. 

PORT  MODE  VOLTAGES 

EIGENVALUE  PORTIll  P0RTI2)  PORTO)  P0RTI4) 

0. 1551H43E  +  00-0.2053380F-01  0 . 1  1 57R2 3E  +  00  0.3424H25E+00  0.9999995E+00 
0.  101 1907E+0?  0. 1904054F  +  00  0. 2323323E +  00  0. 9999999E  +  00-0. 9255374E  +  00 
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0.t5  0b3486E+02  0.9449996E+00-0.2887278E +00  0 . 7RR39 1 6F +00-0 .4044 70RE +  00 
-0.R162234E+03  0.  t>31 3363F  +  00  0.4994999E+00-0.674H93RE+00  0. 2066367E +00 

PORT  MODE  CURRFNT S 

EIGENVALUE  PURT(l)  Pri8T(2)  P()RT(3)  P0RT(4) 

-o. lbbiHoqE+oo-o. 133R342E+00  0.4324440F+00  o. R4 1 R4R4E+00  o. qqqqqqqE +00 
-0.101 1435E+02-0.607R172E+00  0. 9999999 F +00  0.8054470E+00-0.64b7505E+00 
-0.8044034E  +  02  0.9499996E+00-0. 53736R2E+00  0. 1347921E-01  0.  7396740F-01 
0.81641 28E+03  0.E441 235F+00  0.qqqqqqqE+00-0.?R3q342F+00  0.77R1506E-01 


69 


III.  STORAGE  OF  SPECIFIED  AND  SYNTHESIZED  PATTERNS 

The  program  of  this  section  stores  the  specified  and  synthesized  patterns 
on  direct  access  data  set  6  so  that  they  can  be  plotted  by  the  program  on  pages 
104-110  of  [7].  The  activity  on  data  sets  1  (punched  card  input)  and  6  (direct 
access  input  and  output)  is  as  follows: 

READ(1,10)  N,  NT,  NPAT,  NF,  N6,  N7,  NPP 
10  FORMAT  (613,  14) 

REWIND  6 

SKIP  N6  RECORDS  ON  DATA  SET  6 
READ(6) (PP(I) ,  1-1,  NPP) 

DO  17  L-l,  NF 

READ(1,10)  M,  12,  13,  N8,  NIV 
READ (1,45) (NE(I) ,  I-1,M) 

45  FORMAT(20I3) 

READ(1,19) (RE(I) ,  I-1,M) 

19  FORMAT(7E11.4) 

READ(1,19) (AE(I) ,  1-1,  M) 

DO  22  K-l,  N8 

READ (1,23) (CUR(I) ,  1-1,  N) 

23  FORMAT (5E14. 7) 

22  CONTINUE 

17  CONTINUE 

SKIP  N7  RECORDS  ON  DATA  SET  6 
WRITE(6)(G(I),  1-1,  JG) 

There  are  N  ports.  There  are  NPAT  polarizations  and  NT  angles  as  described 
on  page  53  of  [7].  The  same  port  parameters  used  by  the  pattern  synthesis 
program  are  read  in  through  PP.  The  real  and  imaginary  parts  of  the  speci¬ 
fied  desirable  electric  field  pattern  are  read  in  through  RE  and  AE.  In  this 
instance,  the  program  is  more  general  than  necessary  because  the  desirable 
electric  field  pattern  is  specified  in  magnitude  only.  This  magnitude  has 
been  read  in  through  RE  and  AE  has  been  set  equal  to  zero.  For  RE (I)  and 
AE(I),  NE(I)  ■  (K-1)*NT+J  indicates  the  KC^  polarization  to  be  stored  in  PP 
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and  Che  J1*1  angle.  12,  13,  and  NIV  are  Che  same  as  in  Che  paccern  synthesis 
program.  An  approximate  pattern  is  defined  by  its  port  currents  CUR  which 
appear  in  the  printed  output  of  the  pattern  synthesis  program.  For  each  value 
of  L,  first  the  M  points  on  the  specified  pattern,  next  the  NT*NPAT  points  on 
the  approximate  pattern  for  K-1,2,...N8  are  stored  in  G. 

Minimum  allocations  are  given  by 


COMPLEX  PP(NPP) 

DIMENSION  ANG(NT) ,  RE(M),  AE(M) ,  G(JG), 
G2(M),  CUR(N),  NE(Mfl) 

where 


NF 

JG  -  l  (M  +  N8*NPAT*NT) 

L-l 

DO  loop  28  stores  the  specified  pattern  in  both  G  and  G2.  DO  loop  22  stores 
in  G  the  synthesized  pattern  defined  by  CUR.  DO  loop  25  is  necessary  to 
account  for  the  minus  sign  in  (52)  of  [7].  DO  loop  43  obtains  the  NPAT  pola¬ 
rizations  and  DO  loop  26  the  NT  angles.  DO  loop  27  is  able  to  accumulate 
the  normalized  electric  field  in  E  because  PP(J2)  is  the  negative  of  the 
normalized  electric  field  from  a  unit  current  at  the  Ith  port. 


For  the  sample  input  data,  CUR  is  the  set  of  port  currents  which 
radiate  the  pattern  synthesized  by  the  pattern  synthesis  program.  The  syn¬ 
thesized  and  desirable  patterns  are  printed  at  the  12  points  at  which  the 
desirable  pattern  was  specified  in  the  pattern  synthesis  program  plus  at 
the  12  redundant  points  enumerated  in  the  previous  section  as  well  as  at 
360*  for  both  polarizations  for  a  total  of  26  points.  The  synthesized  pat¬ 
tern  is  stored  at  every  2.5  degrees  in  G  for  a  total  of  290  points.  The 
contents  of  G  are  written  on  record  4  of  direct  access  data  set  6. 


LISIING  |)F  PROGRAM  TO  STORK  SPECIFIED  AND  SYNTHESIZED  PATTERNS 
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//  ( 0034, EE, 70S,  1  )  ,  'MAUTZ , JOE ' , R EG !UN=700K 

//  EXEC  WATFIV 

/  /GO .  E  TOEiEOO  1  DO  l)SNAME  =  EE0034.REVl,OTSP  =  0LD,ONIT  =  2314,  X 

//  VOLUME=SER=SU0004,DC«=  IRECFM=VS,HLKSIZE=2596,LRECL=2592,X 

//  RUFNO= 1 ) 

//GO.SYSIN  DD  # 

tJOtt  MAUTZ»TIMF=1«PAGES=30 

COMPLEX  PP ( 323ft) »E 

DIMENSION  ANGI  1 45 ) , RE ( 26 ) , AE ( 26 )  , G ( 4692 ) , G2 < 26 ) , CUR < 4 ) , NE ( 77 ) 

READ! 1,10)  N,NT,NPAT,NF,N6,N7,NPP 

10  FORMA  1(613, 14) 

WRITE! 3,11)  N,NT  ,NPAT ,NF»N6,N7*NPP 

11  FORMAT ( ' 0  N  NT  NPAT  NF  N6  N7  NPP • / l X , 2  I  3 , 1  5 , 3 1  3, 14 ) 

REWIND  6 

IF(N6)  12,12,13 

13  DO  14  J= 1 *N6 
READ! 6) 

14  CONTINUE 

17  READ ( 6 ) ( PP ( I ) , 1=1, NPP) 

WRI T  E I  3 , 16  )  (PP(l ) » I  =  1  *  2  ) 

16  FORMAT ( 'OPP* /( 1 X , 4E  1 1  .4)  ) 

0FL  =  360./ (NT-1  I 
DO  31  J=  1  ,NT 
ANG! J)=( J-l )*DEL 
31  CONTINUE 
JG  =  0 

DO  17  L  =  1 ,  NF 

READ (1,10)  M,12,I3,N8,N1V 
WR I TE ( 3,  18  )  M,I?,I3,N8,NIV 

18  FORMAT (  ' 0  M  12  13  NH  N I  V ' / 1 X , 4  I  3  ,  I  4 > 

HEAD! 1 ,45)  <NE( I  ),I*1.M1 

45  FORK  A  T ( 20 1  3 ) 

WR I  TE( 3.46)  ( NE  <  I  )  ,  1  =  1 ,M) 

46  FORMAT!  'ONE'/!  IX, 20(3)  ) 

NE ( M+ 1 ) =0 

REA(>(1,19)!RE!!),I=1,M) 

19  FORMAT  I 7E 1 1 .4 1 

WRI TE( 3,20)  (RE(I ) , I  *  t ,M ) 

20  FORMAT! 1  ORE  •  /  ( 1X.7E11.4) ) 

READ( 1,19)  (AE( I ) , I  *  1  »M ) 

W'H  I  TE  (  3, 7  1  )  (  AE  (  I  1,1*1, M) 

21  FORMAT! ' OAE • / ( IX ,7E1 1  . 4)  ) 

DO  28  J  *  1 ,  M 

JG=JGe1 

G( JG  )  *  SORT  (RE(J)*RE(J)*AE(J)*AE«J) ) 

G? ( J I *G ( JG ) 

28  CONTINUE 

DO  22  K»1  ,NR 

READ! 1 , 23 ) ( CUR ( I ) » J  *  1 ,  N ) 

73  FORMAT ( 5E 14, 7 ) 

WRI TE ( 3 , 24  )  ( CUR ( I l,!*l,N) 

24  FORMAT! •0CUR'/(1X,5E14. 7)  ) 

IF(NIV)  40,39,40 

39  WRITEI3.37)  K 

37  FORMAT! 'O' ,3X, 'PATTERN  0F',I3,'TH  SET  OF  PORT  VOLTAGES ', 8X ,' SPEC  I F 
1IE0  DESIRABLE  PATTERN* ) 

GO  TO  41 

40  WRI IE! 3,38)  K 
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3«  FORMAT!  '0* ,3X, 'PATTERN  0F',I3,'TH  SET  (IF  PORT  CURRENT S ', flX SPFCI  F 
1IFU  DFS1RAHLF  PATTERN' ) 

1)1)  25  I  =  1  *  N 
CURT  I )=-CUR<  1  ) 

25  CONTINUE 

41  J4=  1 
Jb  =  0 

WRITE! 3*42) 

42  FORMAT!'  ANGLE  KF All E ) ' , 5X , • I  MAG < E  )  •  ,  7X  ,  «  |  E  I  '  , 7X , ' REAL  (  E 1 ' ,  5X , • 

1  I  MAG ( E  )  •  ,  7X, ' IF  I  '  ) 

on  43  LL=  1  ,NPAT 
I.LK=  (  LL-1  )*13*NT+I2 
00  2<S  J=  1 , N T 
F  =0  . 

011  27  I  =  1  » N 
J  2  =  L  LK+ I 

E=F+PP( J2)#CUR<  I  ) 

27  CONTINUE 
LLK  =  LLK+I  3 
JG=JG+ 1 
G  (  JG ) =  C  ARS ( E ) 

.I5  =  J5+1 

IFINE! J4J.NE.J5)  GO  TO  24 

WRITE  (  3,24  )  ANG!  .11  ,E,G(JG)  ,RE(J4),AF  (J4)  ,  G?  ( .14  ) 

.14  =  J4+1 

29  FORMAT!  IX  ,F4. 1 ,4E12.4) 

24  CONTINUE 

43  CONTINUE 
22  CONTINUE 
17  CONTINUE 

IFIN7I  33,33,34 
34  DO  34  .1=1  ,N7 
READ! 4 ) 

34  CONTINUE 

33  WRl TE(4) (G( I  )  ,  1=1 , JGI 
WRITE(3,47)(G< 11,1=1,7) 

47  FORMAT! ' OG ' / 1 X , 7E 1 1 .4 ) 

STOP 

END 

JDA1A 

414b  2  1  2  03234 
24341  10  l  1 

1  13  2b  37  49  41  73  H5  97 109 1 2 1 1 33 1 451 44 1 5H 1 70 1 P2 1 94 2042 1  A 
2 3024 2 2 542 4 42 7 R 290 

0.1000F+01  0.R440F+00  0.5000E+00  O.OOOOE+OO  0.5000E+00  0.R440E+00  0.1000E+01 

D.R440E+00  0.5000E+00  O.OOOOE+OO  0.5000E+00  0.H440E+00  0.1000E+01  0.1000E+01 

(I.8440E+00  0.5000E+00  O.OOOOE+OO  0.5000E+00  0.R440F+00  0.1000E+01  O.R440F+00 

0.5000E+00  O.OOOOF+OO  D.5000E+00  O.R440E+00  0.1000E+01 

O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOF+OO  O.OOOOE+OO  O.OOOOE+lJ  O.OOOOE+OO 

O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOF+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 

O.OOOOF+OO  O.OOOOF+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO 

O.OOOOE+OO  O.OOOOF+OO  O.OOOOE+OO  0 . OOOOE +00 . 0. OOOOE +00 

O.R9R2205E-O1 -  0.39 489 49E-01-0.4 1901 21 E-02-0.  5254191 E-01 
1ST0P 
/  + 

// 

PRINTFD  OUTPUT 


N  NT  NPAl  NF  N4  N7  NPP 


4  14S 


2  1  2  OH?  if. 


73 


MM 

0.  HI  7F-03  0.30.Me-0?-0.3933E-03-0.3|97F-0? 

M  I 2  13  N8  N I  V 
2634  1  10  1  1 


NF 

1  13  25  37  49  (SI  73  85  9  7  1  09  1  2  1  1  33 1  45  14f  I  Sfl  1  701 8?  1 94206?  1  8 
230242254266278290 


HP 

0. 1000F+01 
o.H(Sf,oe  +  oo 
o.H66oe+oo 
o. sooof  +00 

AF 

o.ooooe+oo 

o.oooof+oo 

o.ooooe+oo 

o.ooooe+oo 


o.866oe+oo 

0.5000E+00 

o.soooe+oo 

o.ooooe+oo 


n. ooooe+oo 
o  ooooi+oo 

o. oooue+oo 
o.ooooe+oo 


o.soooe+oo 

O.OOOOE+OO 

O.OOOOE+OO 

O.SOOOE+OO 


O.OOOOE+OO 

O.OOOOE+OO 

o.ooooe+oo 

O.f'OOOE  +  OO 


0.0000F  +00 
0.5000E+00 
0.5000F  +00 
0.  R660F  +  00 


o.ooooe+oo 

O.OOOOE+OO 
0.0000E  +00 
O.OOOOF+OO 


0 .  5<>00t-  ♦  00 
0.  8660E  +00 

o.ronof+oo 
o. iooof+oi 


O.OOOOE  +00 
O.OOOOE+OO 

o.ooooe+oo 

o.ooooe+oo 


0.  R460E  +00 
0. 1000F+01 
0.  1000E+01 


O.OOOOF+OO 

O.OOOOF+OO 

O.OOOOF+OO 


0.  1000F  +01 
0. 1000F+01 
0.H66OE  +00 


O.OOOOF  +00 
O.OOOOF+OO 
O.OOOOF +00 


r.UH 

0.8982206E-01-0.396H959E-01-0.4190121F-02-0.5236I91F-OI 


MA11FRN  OF  1  1  H 

set  lie  port 

CURRENTS 

SPECIFIED  DESIRAHLF 

PATTERN 

ANKLE 

REAL ( E  1 

IMAf.lEI 

IF  1 

BFALtEI 

IMAf.lEI 

IE  1 

O.C 

0.  /414E+00 

0.  14B9E+00 

0.7542F+00 

0. 1000E+01 

O.OOOOE+OO 

0. 1000F  +  01 

30.0 

0.6693F+00 

-0. 1668E+00 

0. 48988+00 

0.8660E+00 

O.OOOOE+OO 

O.P660F+00 

60.0 

0.237HE-01 

-0.5134E+00 

0.6140E+00 

0.5000F+00 

O.OOOOE+OO 

0.5000F+00 

9Q.0 

-0.3830F+00 

0. 1046F+00 

0. 3970E+00 

O.OOOOE+OO 

O.OOOOE+OO 

o.ooooe+oo 

120.0 

0. 1585E+((0 

0.4173F+00 

0.4464F+00 

0.5000E  +00 

O.OOOOE+OO 

0.5000F+00 

150.0 

0 . 5  3  84  E  +00 

0.9461E-01 

0.5466E+00 

0.8660F+00 

O.OOOOE+OO 

0.R660E  +00 

180.0 

0.5H06E+00 

-0. 1 1 16E+00 

0.5912E+00 

0. 1000E+01 

O.OOOOE+OO 

0. 1000E  +  01 

210.0 

0.5384E+00 

0.9461E-01 

0 . 5466E  +00 

0.8660E+00 

O.OOOOE+OO 

0.8660F+00 

240.0 

0.  15R5E+00 

0.41 7  3E  +00 

0.4464F  +00 

0.5000E+00 

O.OOOOE+OO 

0.5000E+00 

270.0 

-0.3830F+00 

0. 1046E+00 

0. 3970E+00 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOF  +00 

300.0 

0.237PE-01 

-0.5134E+00 

0.5140F+00 

0.5000E+00 

o.ooooe+oo 

0.5000F+00 

130.0 

0.6693E+00 

-0. 166HE+00 

0.6898E+00 

0.8660E+00 

O.OOOOE+OO 

0.8660F+00 

360.0 

0 , 74 1 4E  +00 

0. 1489F+00 

0.7562F+00 

0. lOGOE+Ol 

O.OOOOE+OC 

0. 1000E+0) 

0.0 

0,74 1 4F  +  00 

0. 1489F+00 

0.7562E+00 

0. 1 000E+01 

O.OOOOF+OO 

0. 1C00F+01 

30.0 

0.78R4E+00 

-0. 3P22E+00 

0.8762F+00 

0.8660E+00 

O.OOOOF+OO 

0 • 8  660F ♦ 00 

60.0 

-0.2H37E+00 

-0.H580E +00 

0.9036F+00 

0.5000E+00 

O.OOOOE+OO 

0.5000F  +00 

60.0 

”0 .  9  676E  — C 1 

0. 3541F+00 

0.3671 F  +  00 

O.OOOOF+OO 

O.OOOOE +00 

O.OOOOE+OO 

120.0 

0.32H6F+00 

-0 . 6705E+00 

0.7467E+00 

0.5000E+00 

O.OOOOE+OO 

0. SOOOF  +00 

150.0 

-0.6147F  +00 

-0.3531F+00 

0. 7089F  +  00 

O.R660F+00 

O.OOOOF+OO 

0.8660F+00 

180.0 

-0.5806E+00 

0.  1  1 16F+00 

0.5912E+00 

0. 1000F+01 

o.ooooe+oo 

0. 1000E  +01 

210.0 

-0.6147E+00 

-0.3531F+00 

0.7090F+00 

0.8660F  +00 

o.ooooe+oo 

0.8440F+00 

240.0 

0.32H6F+00 

-0.6705F+00 

0.  7467E  +00 

0. 5000E+00 

O.OOOOE  +00 

0 . SOOOF +00 

270.0 

-0.9675E-01 

0. 354 1 E  +00 

0 . 367 1 F  +00 

O.OOOOF+OO 

O.OOOOE  +00 

o.ooooe+oo 

300.0 

-0.2837E+00 

-0.HS/9E+00 

0.9036F  +  00 

0.5000F+00 

O.OOOOE+OO 

0. SOOOF  +00 

330.0 

0. 7H84F  +00 

-0. 3822E+00 

0.8762E+00 

0.8660E+00 

O.OOOOE+OO 

0.8660E+00 

360.0 

0.74 14e+00 

C. 1489F+00 

0. 7562E+00 

0.  1000E  +  01 

O.OOOOE+OO 

0. 1000F+01 

r. 

0.1000E  +  01  0 .  R660F  +  00  O.SOOOE+OO  O.OOOOF+OO  O.SOOOE+OO  0.8660F+00  0.1000F+01 
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IV.  LOADS  FOR  MODAL  RESONANCE 

The  prograa  (nodal  raaonanca  program)  of  ehla  aacelon  calculataa  tha  aat 
of  raactlva  loada  (14)  which  aakaa  a  glvan  raal  port  currant  an  algancurrant 
vhoaa  alganvalua  la  zero.  Tha  activity  on  data  aaca  1  (pun chad  card  input)  and 
6  (dlract  accaaa  Input  and  output)  la  aa  followa. 

READ(1, 10)  NF 
10  FORHAT(13) 

DO  12  L-l,  NF 
READ(1,13)  N,  N6,  N8,  12 
13  FORMAT(3I3,  14) 

N2  -  N*N 
NPP  -  12  ♦  NZ 
REWIND  6 

SKIP  N6  RECORDS  ON  DATA  SET  6 
READ (6) (PP(I) ,  1-1,  NPP) 

DO  20  J-l,  N8 

R£AD(1 ,21) (CUR(I) ,  1-1,  H) 

21  FCRHAT(5E14. 7) 

20  CONTINUE 

12  CONTINUE 

Tha  raactanca  matrix  Xg  In  (14)  la  tha  Imaginary  part  of  tha  lmpadanca  matrix  found 
In  PP(I2+1)  through  PP(I2+N*N).  Tha  Jtl‘  aat  cf  raal  port  curranta  1  of  (14) 
la  raad  In  through  CUR. 

Minimum  allocatlona  ara  glvan  by 

COMPLEX  PP(NPP) 

DIMNSION  X(N*N),  CUR(N),  XL(N) 

DO  loop  19  atoraa  IXg J  of  (14)  In  X.  DO  loop  26  puta  the  largaat  |ltl  of 
(14)  In  CU.  DO  loop  23  atoraa  Xj  of  (14)  In  XL(X).  Statamant  27  anauraa  that 
1^  In  tha  danoodnator  of  (14)  la  at  laaat  CU*l.E-8. 

Tha  aampla  output  XL  la  tha  aat  of  raactlva  loada  for  modal  reaonanca  of  tha 
aat  of  raal  port  curranta  which  radlata  tha  pattarn  aynthaaizad  by  tha  pattarn 
aynthaala  prograa. 


/ 


I  HI  INC.  Ilf-  Ml  II 1  At  R  p  M  IN ANC  f-  PHI 1C  H  AM 
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//  (  OOH.  Ef-  .  I  O'..  1  )  .  •  M AU I /  ,  JIIF  1  ,HPr.|t)N»700K 

n  me  mifiv 

//l.l'.HONPOOl  Dll  (ISNAMP  .PP0034.HEV1  ,01  5P*OLI),UNI  T  »  2314  ,  X 

II  VUlUMP-SPR»SU0004,liCh»<PFCFM-VS,rtLKSI/F»2*94,LRFCL*2‘>92.x 

//  Ml  IP  M I  *  I  1 

//r.n.sxsiN  no  • 

*  Jl  IH  M  All  I  7  t  T  I  pfc  ■  I  «PAf.ES*20 

r i imp t P x  pp(wno) 

DIMENSION  XI 900), CURT  30  I, XL  I  301 
P  P  All  (  1,10  1  NP 
10  PIIHMA1  1131 

MM |  IP  I  3, 1  l  I  NP 
1  1  PlIHMAll  'ONP.*  ,131 
Dll  l  2  I  «  I  .  NP 
RP  A(l(  1,13!  N.N4.NH, |  2 

13  PIIMMAl  I  313.14  1 

HR  lit  I  3.  14  I  N.N4.NH ,  | ? 

14  PORMAll'O  N  N4  NH  I  ?  '  /  1 X , 3  I  3 . I  4  I 
N7  «N*N 

NPP» I ?«NZ 
KpHlNi)  4 
IKN4I  1S.1S.1P. 

1  4  Dll  17  Jm  1  ,N4 
RPAOI 4) 

17  CONTINUE 

15  RFAOIMIPPI  I  I  .  I«1  ,NPP» 

HR  IIP!  3,  1«  MPPI  I  1  ,  l«l,3l 

|R  KIRMATI  'OPP'/I  1X.10E11.4I1 
iki  19  i « i . ni 
J7*l 2*1 

XI  I  I  «  A  I  MAC,  I  PP  (  J7  I  1 

19  CONTINUE 

no  20  J«1.NH 

READ! 1 ,21  I (CUR  I  I  1 , I  * l ,N 1 
21  FORMA  T  I  t>F  14 . 7  I 

HR|  TE(3,2?HCUR(  II  ,  T*1  ,NI 
27  FORMAT  I  •  OCllR  •/!  lX»bE14.7)  I 
fU'O. 

1)0  24  I  ■  1 ,  N 
S1*AHS(C'»R(  I  )  I 
IFISl.r.T.CU)  ClI-SI 
24  CONTINUE 

CUR«C0*1 ,F-R 
DO  23  I '1  ,N 
Jl«l 1-1 1 -N 
Sl-O. 

1)0  24  K  *  1  ,  N 
J  2* J 1 *K 

S1*S1*XI  J2I-CURTKI 

24  CI1NTINUF 
S2«CUR( I  I 

27  IKARSIS2I  .LT.CUfl)  S2-CU8 
XL  I  I  I—S1/S2 
23  CONTINUE 

NRI  TEI3.25  )  imill  .  1  =  1  ,N1 

25  FORM  A 1 1  •OXL'/I  IX, *>614.71  1 

20  CONTINUfc 
17  CONTINUE 
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STOP 

END 

SOA1A 

1 

•4  2  11^ 

0.8982205E-01-0. 39684 69E-01-0. 4 19012  IE-02-0. 5256 19  IE-01 
$S  TOP 
/* 

// 

PRINTED  OUTPUT 

NE=  1 

N  N6  N8  12 
4  2  1  16 


PP 

0.1917E-03  0.306  IE -02-0.  39  33E-03-0.  319 7E-02-0. 637 IE-03-0. 2i> .  oE-03 
COR 

0. 898  2205E-01 -0.39684  69E-01 -0. 419012  IE-02-0. 5256 19  IE-01 
XL 

0. 37279 10E  +  03  0.6710123E  +  02-0. 1086468E+05  0.6431934E+03 
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V.  RADAR  CROSS  SECTION 

The  program  (o/X2  versus  angle  program)  of  this  section  calculates  the 
radar  cross  section  per  wavelength  squared  o/X2  of  the  N-port  loaded  scatterer. 
For  a  given  Incident  plau.3  wave,  this  program  obtains  patterns  of  o/X2  In  one 
or  more  of  the  coordinate  planes.  The  activity  on  data  sets  1  (punched  card 
Input)  and  6  (direct  access  Input  and  output)  Is  as  follows. 

READ(1,10)  NF,  N6W 

10  FORMAT (213) 

DO  12  JF-1,  NF 

READ(1,13)N,  N6,  NQ,  NS,  NL,  NT,  NA,  NPAT,  NIV,  NZ,  12,  13,  BK 
13  FORMAT (9 13,  314,  E14.7) 

REWIND  6 

SKIP  N6  RECORDS  ON  DATA  SET  6 
NPP  -  12+13*  (NT*NPAT-1)  +  N+l 
READ(6) (PP(I) ,  1-1,  NPP) 

DO  22  JL  -  1,  NL 
READ(1,24) (XL(J) ,  J-l,  N) 

24  FORMAT (5E14. 7) 

Z2  CONTINUE 

12  CONTINUE 

SKIP  N6W  RECORDS  ON  DATA  SET  6 
WRITE(6)  (SIG(I)  ,  >1,  J7) 

Virtually  all  of  the  main  program  Is  Inside  DO  loop  12.  Referring  to 
the  table  on  page  62  of  [7],  the  given  Incident  plane  wave  Is  specified  by 
Voc  and  F^c  •  u*  which  reside  In  PP(I2+(NA-1)*I3+1)  through  PP(I2+(NA-1)*I3+N+1) . 
For  each  of  NPAT  patterns,  a/X2  is  evaluated  at  NT  points  but  is  printed  at 
only  the  first,  the  (NS+l)t^1,  the  (2*NS+l)t^1, . . .  of  these  points.  The  N-port 
impedance  matrix  Z  at  propagation  constant  BK  resides  in  PP(NZ+1)  through 

->  HQ 

PP(NZ+N*N)  while  V°c  and  F  •  u  are  in  PP(I2+((J-1)*NT+I-1)*I3+1)  through 

'^o  /*r  .  , 

PP(I2+((J-1)*NT+I-1)*I3+N+1)  for  the  I  value  o/X2  on  the  J  pattern.  The 
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JLth  set  of  reactive  loads  is  read  Into  XL  inside  DO  loop  22.  N1V  t*  0  obtains 
the  open  circuit  impedance  formulation  while  N1V  ■  0  obtains  the  dual  short 
circuit  admittance  formulation.  If  NQ  is  neither  1  nor  2,  both  •/ajx7  and 
o/X2  are  stored  in  SIG.  If  NQ  -  1,  only  /o/x^  is  stored  in  SIG.  If  NQ  ■  2, 
only  o/X2  is  stored  in  SIG.  All  of  the  data  (either  /o/X7  or  o/X2)  for  a  given 
set  of  loads  is  stored  in  a  block  in  SIG.  For  the  first  set  of  loads,  /o/x2  at 
the  I**1  point  on  the  pattern  is  put  in  SIG((J-1)*NT+I)  and  the  corresponding 
o/X2  is  put  in  SIG(NT*NPAT+(J-1)*NT+I)  provided  NQ  is  neither  1  nor  2. 

Minimum  allocations  are  given  by 

COMPLEX  C(N*N) 

DIMENSION  LR(N) 

in  the  subroutine  L1NEQ  and  in  the  main  program  by 

COMPLEX  PP(I2+I3*(NT*NPAT-1)+N+1),  ZS(N*N),  Z(N*N),  CUR(N) 

DIMENSION  ANG (NT),  XL(N),  ND(N) ,  D(N) ,  El (NT*NPAT ) , 

E2(NT*NPAT)  ,  SIG(J7) 

DO  loop  21  puts  Zg  of  (1)  in  ZS.  The  JLth  set  of  loads  is  considered  in 

DO  loop  22.  DO  loop  23  puts  ZS  in  Z.  DO  loop  26  adds  the  loads  to  the  diagonal 

til 

elements  of  Z.  If  the  ratio  of  | XL ( J ) j  to  the  magnitude  of  the  J  diagonal  ele- 

til  til 

ment  of  ZS  is  greater  than  10,  then  DO  loop  29  divides  the  J  row  and  the  J 
column  of  Z  by  the  square  root  of  this  ratio.  Statement  28  inverts  Z.  DO  loop  32 

is  similar  to  DO  loop  29.  DO  loops  29  and.  37  scale  [11]  the  matrix  Z  to  avoid 

excessive  round  off  error  in  the  subroutine  LINEQ.  Do  loops  29  and  32  have  no  net 

effect  in  the  absence  of  roundoff  error.  DO  loop  34  puts  [Z  +  Z  ]  *Voc  of  (1)  in 

CUR.  DO  loop  50  is  necessary  because  (1)  is  combined  with  (52)  of  [7]  while,  for 
the  dual  short  circuit  admittance  formulation,  (7)  is  combined  with  (60)  of  [7]. 
The  J  o/X2  on  the  LLth  pattern  is  obtained  in  nested  DO  loops  38  and  39.  In 
DO  loop  39,  PP(J5)  represents  Eoc  of  (1).  DO  loop  39  stores  o/X2  in  E2  and 

_  MO 

/o/x2  in  El. 

The  loads  XL  for  the  sample  input  data  are  taken  from  the  printed  output 
of  the  modal  resonance  program.  The  incident  electric  field  is  an  x  polarized 
plane  wave  traveling  in  the  z  direction.  According  to  the  figure  on  page  42  of 
[7],  this  plane  wave  is  incident  on  the  tip  of  the  wire  triangle.  The  two 
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patterns  appearing  In  the  printed  output  are  the  second  and  third  polarizations 
of  (103)  of  [7].  For  proper  comparison  with  the  |e|  patterns  of  the  pattern 
symthesis  program,  only  /o/Az  is  stored  on  record  5  of  data  set  6.  The  SoJ A7 
patterns  on  record  5  of  data  set  6  can  be  plotted  by  the  program  on  pages 
104-110  of  [7] . 


80 

LISlINf.  OF  Sir, HA  OVER  LAMBDA  SOI  J  A  R  tf  I J  VF  R  SI  IS  ANGLE  PROGRAM 

//  (0034  ,FF  .70S,  1  )  ,  'MAIIIZ  ,.IOF  '  •  PEG  |  (IN  -  200K 

//  EXEC  WA1EIV 

//GO.E 108F001  1)1)  I)SNAME  =  EE0034.REV1  , I)  I  SP=OLI) ,IIN I  T  =  23  14 ,  X 

II  VOLUME  =  SER  =  SU0004,IICR  =  (  RECEM=  VS,  RL  KS  I  7  E  =  2c>9h ,  LR  FCL  =  2^42 ,  X 

//  HI  IF  NO  =  l  ) 

//GO.SVSIN  1)1)  * 

tJOH  MAUTZ  «T1MF=1, PAGES =40 

SIIHRIHMINF  L  I  NEO  (  LL  ,C  ) 

COMPLEX  mom  ,STOR,STO,ST,S 
DIMENSION  LR I  40  ) 

DO  20  1  =  1  .LL 
LR (  I  )  =  I 
20  CON  1 1  NOE 
M  l  =  0 

DO  1R  M= 1  ,LL 
K  =  M 

DO  2  I =M,LL 
K 1 =M 1 + I 
K?  =  M  l  +  K 

IFICAHSICIKl 1 I-CAHSICIK2)  ) )  2,2,8 
4  K  =  I 

2  CONTINUE 
LS=LR(M) 

LR (MI=LR(K) 

LR(K)=LS 
K2  =  Ml +K 
S 10R  =  C ( K2  ) 

J1  =  0 

DD  7  .1=1, LL 
K 1  =  J 1 +  K 
K2=J 1+M 
ST0=C(K1  ) 

CIK1 )=C(K2) 

C(K2)=Sin/ST0R 
J1=J1+LL 
7  CONTINUE 
K  1  =M  1 +M 
C(K1)=1./ST0R 
DO  11  1  =  1, LL 
I F (  I -M  )  12,11,12 
12  K 1 =M 1+ I 
ST  =  C(K1  ) 

C ( K 1 ) =0 . 

Jl=0 

DO  10  J=  1  ,LL 
K  1  =  J 1  +  I 
K  2  =  J 1 +M 

C(Kl)=C(Kl)-C<K2>*ST 
J 1 = J 1+LL 

10  CONTINUE 

11  CONTINUE 
M 1 =M1 +LL 

18  CONTINUE 
J1=0 

DO  4  J  =  1 ,  LL 
IFIJ-LRIJI)  14,8,14 
14  LRJ=LR(.I) 

J2=ILRJ-l)UL 
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21  (Hi  13  1=1  .LL 
K2  = J2+I 
K  1  =  J  1  ♦  1 
S=C(K2I 
c  <K2i=r.  iKi) 

CtKl  i  =  s 
13  f.llNl  I  NUE 

LR ( J)=LR!LRJ) 

LR  I  LR  J  )  =LR.I 
IHJ-LRIJII  lA.fi, 14 
fi  J1=J1+LL 
A  C<  >N  T 1  NUE 
RElURN 
END 

COMPLEX  U.PP 13234)  ,7  SHOO)  ,Z  (  100 ), ClIR  110  I  , E  .CON  JO 
01  MENS  10M  AND  (  1 4b  ) , XL ( 10 ) «N0(  10)  , f) <  101 ,  E  1  (  290  1  ,  E 2  (  290  )  ,  S  I G  1  4440 1 
E  1  A=  37 4 . 730 
PI=3.l41b43 

C  1  =  .  2S»ET  A/SORT  (  P!  *>PI*PI) 

0=10., I, ) 

R  FAI)(  1  »  10  1  NF.N4W 

10  FORMA  T ( 2  I  3  I 
WR1TE13.11)  NF.N4W 

11  FORMAT ( '0  NF  N4W • / 1 X , I  3 , 1 4 ) 

.17  =  0 

(HI  12  JF  =  1 , NF 

READ! 1 ,13)  N.Mfi.NO.NS.NL.NT .NA.NPAT ,N I V ,NZ , !?, 13.BK 

13  FORMAT |RI3,3!4,E14.7) 

WRITE! 3, 14 )  N,N4 ,NQ,NS,NL,NT .NA.NPAT ,NIV»N7 , 12, I3.RK 

14  FORMAT ( *0  N  NO  MO  NS  NL  NT  NA  NPAT  NIV  N7  12  I  3 ' , AX , • BK • / 1 X , 7 1 
13.I5.4I4.E14.7) 

NP 1 =N+ 1 

JP=I 2+INA-l 1*13 
NTP=Nl*NPAT 
OEL  =  340 . / ( NT- 1 ) 

UO  lb  J=1.NT 
ANGl  J)  =  (  J-l  )*I)FL 
lb  CONTINUE 
C2  =  RK’>RK*C1 
REWIND  4 
I F ( N4 )  14,14,17 
17  00  18  J= 1 , N4 
READ! 4) 
lfi  CONTINUE 

1 4  NPP=I2+I 3*<NT*NPAT-1  l+NPl 
READ! 4) (PP ( I ) , I = 1 , NPP ) 

WR! TE(3, IP) <PP< I )  ,  1  =  1 ,2) 

14  FORMAT! 'OPP'/< 1X.4E11.4) ) 

N1P=NT#NPAT 
NN=N*N 
00  21  J=1 .NN 
J 1 =NZ  + J 
7  S ( J  )  =PP ( J 1 ) 

21  CONTINUE 

00  22  JL= l ,NL 
00  23  .1=1, NN 
7  (  J )  =  7  S  I  J  ) 

23  CONTINUE 

READ! 1,24) (XL! J) ,J=1,N) 

24  FORMAT (bE 14.7  ) 
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WRI  1E<3,25)  (  XL  ( .1 )  ,  J=1,N) 

25  FORMAT!  '()XL'/(  1X.5F14.7)  ) 

J  1  =  1 

.12  =  0 

(HI  26  J=  1  ,N 
Z  (  .11  )=Z  (  J1  )+U*XL< .)) 

S1=ABS( XL ( J)  » /CAMS  <  7  S I J1 )  I 
J1=J1+NP1 

IEIS1-10.)  26,26.27 
27  J2  =  .I2  +  1 
NIX  J  2  )  =  J 

IX  J2)=1./S0RT(S1 ) 

26  CONTINUE 

I F ( J2 . EQ.O )  on  TO  7H 
nil  24  J  =  1 , .12 
.11  =NIX  .)  I 
J  3  = ( J 1  —  1 )*N 
no  30  I  =  1  ,N 
J  3  =  J  3  +  1 

7  ( .13  )  =Z  (  J3  I  *l)(  .1 1 
Z  (  J1  )=Z  (  J1  )  *1)1.1) 

.1 1  =  J 1  +N 

30  CONTINUE 
29  CONTINUE 

2R  CALL  L I NEO ( N , Z ) 

IFIJ2.E0.0)  r,n  TO  31 
DO  32  J=  1 ,  .12 
J 1  =NIX  J  ) 

.1 3  =  (  Jl-1  )*N 
00  33  I  =  1 ,  N 
J3=J3+ 1 

7! J3)=Z( J3)*D(J> 

Z  (  J1  )=Z  (  J1  )*IH.J) 

J 1  =  J 1  +N 

33  CONTINUE 
32  CONTINUE 

31  00  34  .1  =  1, N 
J3=< J-l  )*N 
CUR  I J ) =0. 

00  35  I  =  1  ,  N 
.12  =  1  +JP 
J4=J3+I 

CUR  (  J  I  =CUR  (  J  )  +7  <  J4  )  »PP  ( .12  ) 

35  CONTINUE 

34  CONTINUE 
IF(NIV.EO.O)  GO  TO  49 
00  50  J=1 ,N 

CUR ( J ) =-CUR ( J ) 

50  CONTINUE 
49  WRITE  (3, 36)  ,)L 

36  FORMAT ( 'OSCAT TEH ING  PATTERN  FUR  THE',I3,'TH  SET  OF  LOADS') 

IHR  I  TE  (  3 , 37  ) 

37  FORMAT (  1  ANGLE  RE AL ( E ) • , 5X , • I  MAG ( E ) ' , 7X , • | E I • , 6X , ' S I G/ ( L AM ) **2 ' 
1  ) 

J  1  =  I  2 
J6  =  0 

1)0  3R  LL=  1  ,NPAT 
00  39  J=!,NT 
.15  =  J  1  +  NP  1 
E  =  PP  ( .15  ) 


i 

.1 

i 

vs 


83 


DM  4  n  I  =  1  ,  N 
I?=J1+I 

E  =  I — CUR  (  I  )*PP(  .1?  | 

40  CONTINUE 
h  =  fc  **  C  2 

J  (S  =  J4+I 

t-?(  jM=»-»ciiN.ir,(e  i 
H  (  JA  >  =  sow  T  (  E  2  (  J4)  ) 

J  3  =  J  - 1 
J  1  =  J 1  + 1  3 

IE(J3/NS*NS.NE.J3)  f,G  TO  39 

WRI  It  13.41  )  ANG  (  J  (  ,E.E1  (  JM  ,E2<.)4) 

41  FORMAT!  IX, Eft. 1  ,4E 1 2 .4  I 
34  CONTINUE 

38  CONTINUE 

1EIN0.E0.?)  CO  TO  42 
1)0  43  J  =  1  «NTP 
J7= J7+1 
SI0(J7)=F1!J) 

43  CONTINUE 

42  I E (NO.FO. 1 )  CO  TO  2? 

DO  4b  .1  =  1, NTP 

.1 7  =  J  7+  l 
SIG(.I7)=E2(J) 

4b  CONTINUE 
22  CONTINUE 
17  CONTINUE 

IF!N4W)  44,44,47 

47  DO  43  .1  =  1, N4W 
RE  ADI  4  I 

48  CONTINUE 

44  WR  I  TE ( 4 ) ( SI G(  I  ) , I  =  1 , J7 ) 

WRITE  I  3.51 1  (SI  GUI  ,1=1,4) 

b 1  FORMAT! • OS  I G • / ( IX, ftE 11,4)1 
STOP 
END 

IDA  1  A 
1  1 

4  7  1  17  114b  73  7  1  14  34  1  10  0 . 1 94349bE  +  00 

0.37774 10E+03  0. 47 1 0 1 2 3E +02-0. 1 08 4448E+05  0. 443 1 934E+03 
4. S10P 
/* 

// 

PRINTFD  OUTPUT 

NF  N4W 
1  1 

N  N4  NO  NS  NL  NT  NA  NP AT  NIV  NZ  17  13  BK 

4  2  1  12  1 1'  5  73  7  1  14  341  10  3 . 1443495E  +  00 


PP 

0.  19176-03  0.304 IE -02-0.  39  33E-03-0. 31 9  76-0? 

XL 

0. 37 279 10E  +  03  0. 67101236  +  02-0. 1084448E+0b  0.44319346  +  03 

SCATTERING  PATTFRN  FOR  THE  ITH  SFT  OF  LOADS 
ANGLE  RtAL(E)  I  MAG ( E )  IE!  SIG/(LAM)**2 
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0.0  0.651  71  +00 

30.o  o .  5  39of  +oo 
60.0  -0.1  ?<S0h -0  l 
40.0  -0.38646+00 
170.0  0 .  11046  +  00 

180.0  0. 54616  +  00 

180.0  0.56646+00 

710.0  0.54616+00 

740.0  0.11046+00 

770.0  -0.38646+00 
300.0  -0.12606-01 
330.0  0.53906+00 

37,0.0  0.651  7 F  +  00 

0.0  0,6517 F  +  00 

30.0  0. 57896+00 

7,0.0  -0 . 4 7  1 3 F  +00 
40.0  0. 42286-01 

170.0  0.18336+00 

150.0  -0. 7,78  2F  +00 
180.0  -0.56696+00 
710.0  -0 . 678 2E  +00 
740.0  0.18326+00 

770.0  0.4229E-01 

300.0  -0.47 1 3F+00 
330.0  0.  57846+00 

360.0  0.65176+00 


-n.'.os  1F-01 
-0.27906+00 
-0.48301  *00 
0.376/6-01 
0.4777F+00 
0.54476-01 
-0. 18606+00 
0.54476-01 
0.47771+00 
0.3747F-01 
-0.48306+00 
-0  •  774  OF  +  00 
-0.50516-01 
-0.5051 F  — 0 1 
-0.  57476+00 
-0 . 68  37F  +  00 
0.35406+00 
-0.75706+00 
-0 .76628+00 
0. 18506+00 
-0 .  25576  +00 
-0.75706+00 
0. 3 64 01- +  00 
— 0 • 68  376+00 
-0.5?426+00 
-0 . 506 1 E-0 1 


0.66376+00 
0.7,0648  +  00 
0.4  8  3 1 F  +00 
0. 38H7E+00 
0.4  366F+00 
0. 544  3E  +  00 
0.5463E+00 
0. 6443F+00 
0.4366E+00 
0. 3HH2E+00 
0.4P31E+00 
0 • 6064F  +00 
0.66376+00 
0.6537F+00 
0.7843E+00 
0.8300F+00 
0. 35666+00 
0.77406+00 
0.7246F+00 
0.6463F+00 
0 • 7246F+00 
0.7740F+00 
0.3566E+00 
0.83006+00 
0.78436+00 
0.65376+00 


0.4273F+00 
0. 3684F+00 
0.7334E+00 
0.15076+00 
0. 14066+00 
0. 30176+00 
0.36566+00 
0.30176+00 
0. 14066  +  00 
0. 15076+00 
0. 7334F  +  00 
0.36846+00 
0.42736  +00 
0.42736+00 
0.61516+00 
0. 68H9F+00 
0.17716+00 
0.5441E+00 
0.5251  E  +00 
0.35566+00 
0.525  IE +  00 
0. 59416+00 
0.  1271F  +  00 
0. 68H9E+00 
0.6151E+00 
0.4273E+00 


sir, 

0.65376+00  0.65 3 36+ 00  0.6574E+00  0.65086+00  0.64856+00  0.64566+00 


i 


& 

A 
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VI.  OPTIMUM  GAIN 

The  program  (optimum  gain  program)  of  thla  section  obtains  the  set  of 
real  port  currents  (50)  which  maximizes  (55).  The  activity  on  data  sets  1 
(punched  card  input)  and  6  (direct  access  input  and  output)  is  as  follows. 

READ(1,7)  NF 
7  FORMAT (2013) 

DO  9  JF  -  1,  NF 

READ (1,10)  N,  N6,  N8,  NPP,  NZ,  NV,  BK 
10  FORMAT (313,  314,  E14.7) 

READ(1,7) (NST(I) ,  I  -  1,  N8) 

REWIND  6 

SKIP  N6  RECORDS  ON  DATA  SET  6 
READ(6) (PP(I) ,  I  -  1,  NPP) 

NN8  -  N*N8 

READ (1,17) (FI (I),  I  -  1,  NN8) 

17  FORMAT (4E14. 7) 

9  CONTINUE 

Virtually  all  of  the  main  program  is  inside  DO  loop  9.  For  propagation 
constant  BK,  the  N-port  impedance  matrix  Zg  of  (1)  is  in  PP(NZ+1)  through 
PP(NZ+N*N)  while  Voc  of  (36)  is  in  PP(NV+1)  through  PP(NV+N).  The  nth  basis 
vector  of  port  currents  resides  in  FI((n-l)*N+l)  through  FI(n*N).  If 
NST(J)  -  0,  then  the  optimization  of  (55)  with  J  basis  vectors  (II-l)  is 
omitted. 

Minimum  allocations  are  given  by 

DIMENSION  LR(N8),  C(N8*N8) 

in  the  subroutine  LINER  and  in  the  main  program  by 

COMPLEX  PP(NPP) ,  VC(N8) 

DIMENSION  NST (N8) ,  FI(N*N8),  R(N*N) , 

RI(N*N8) ,  RCS (N8* (N8+1 ) / 2 ) ,  VI (N8), 

V2(N8),  RC(N8*N8) ,  A1(N8),  A2(N8), 

ALP(N8),  CR1(N) ,  CR2 (N) ,  CUR(N) 
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DO  loop  19  puts  [R]  of  (34)  in  R.  DO  loop  20  puts  (R]ln  of  (54)  in 
RI((n-l)*N+l)  through  Rl(n*N).  DO  loop  23  puts  [R]  of  (55)  in  RCS  using  the 
symmetric  tnMe  of  storage  [12],  DO  loop  27  puts  the  real  and  imaginary  parts 
of  V°C  of  (55)  in  VI  and  V2. 

DO  loop  30  maximizes  (55)  using  J  basis  vectors  (11-1).  DO  loop  31  puts 

R  in  RC  by  columns.  Statement  52  Inverts  R.  DO  loop  33  puts  the  real  and 

—  1 A  oc 

imaginary  parts  of  [R j  V  in  A1  and  A2.  DO  loop  35  and  statement  53  put  the 
maximum  gain  for  complex  a  in  GC.  The  constant  appearing  in  (40)  resides 
in  Cl.  The  maximum  gain  for  real  a  will  be  put  in  GR.  If  all  of  the  elements 
I  of  either  Re(V  )  or  Im(V  )  are  zero,  then  GR  ■  GC,  C  ■  1,  and  the  evaluation 

of  (45)  is  avoided.  It  was  permissible  to  set  C  ■  1  because  (42)  does  not 
depend  on  C  when  Voc  is  constant  phase.  The  logic  between  statements  36  and 
38  puts  (44)  in  C  and,  with  the  help  of  (42),  puts  the  maximum  gain  for  real 
a  in  GR.  Upon  exit  from  DO  loop  40,  the  real  a  for  maximum  gain  is  in  ALP 
while  the  real  and  Imaginary  parts  of  the  complex  a  for  maximum  gain  are  In 
A1  and  A2.  DO  loop  48  puts  the  real  port  current  I  for  maximum  gain  in  CUR 
after  storing  the  real  and  imaginary  parts  of  the  complex  port  current  I  for 
maximum  gain  in  CR1  and  CR2. 

-►oc 

For  the  sample  input  data,  V  is  the  set  of  open  circuit  port  voltages 
for  an  x  polarized  unit  plane  wave  traveling  in  the  positive  z  direction.  The 
basis  vectors  (II-l)  are  the  three  dominant  (X  -  -  0.1552,  X  ■  -  10.12,  and 
X  ■  -  50.54)  mode  currents  for  the  unloaded  wire  triangle. 


IIMlNf  Ilf-  IIP  1  |  Ml  IM  I. A|N  PROGRAM 
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/  /  i  0034  . M-  ,20s.  1 1 .  •  h aii t / ,  ji ip  1  ,«pr. iiiN.?orn' 

/  /  p  xFf  HA  in  v 

//f.ll.P  10FF001  I'll  IlSNAMP  «PP0034.RP  VI  .1)1  SP»!IU'.MNt  I  =2  i|4,  » 

/  /  VI  III  IMP  .  Mi  M.  <,1)0004  ,  I  )C  H  ■  IPPCFM.  vs.  Ml  *S  I  7  P  >  2  S<)4  ,  L  “  P  C  L  ■?•>■»?.  * 

//  Hi)FNO  =  1  I 

//cii.sysin  do  * 

S.IOK  MAUTZ  ,  T  I  MF  U  1  .PACE  s»  30 

c  this  program  talus  thp  .'atrix  inversion  surriiih  inf  linpr  listed 

r  W I  1 H  1HP  PATTfcPN  SYNlHFSIS  PRtlORAM  HP  SECTION  f| 

COMPLEX  PP|  323F)  .  V  f.  (  4  ) 

DIMENSION  NSr(4),F|(lA),R(liS),RI(lN),Rf. SI10), Vl(4),V?  141, RCIIM 
f)l  MENS  I  ON  A1I4  |  ,  A?(4)  ,ALP(4  )  ,CR|  (4  )  ,CR?I4  )  ,CUR  14) 

P|=3.141SR3 
E1A=37F.730 
C?  =  ET A/(4.*PI  ) 

R  P  A  I)  (  1  , 7  )  NE 
7  FORMA  1(20131 
WHITPI 3,P )  NF 
H  FORMA  1 (  '0NH=* ,13) 

DO  0  jr-  =  1  , NF 

RI  AD!  1,10)  N, NF, NR, NPP, NZ ,NV,HK 

10  FORMAT! 313, 314, E14. 7) 

HR  I  TP  I  3, 1 1)  N,NF,NH,NPP.NZ  ,  MV  .  HK 

11  FORMAT (  •  0  N  NF  MR  NPP  NZ  N'  • , NX , • RK 1 / 1 < , 3 1 3 , 3 1 4 , 6 14 . 7 ) 

RFAIX  1 , 7  )  (  M  ST  (  I  1,1=1, MA) 

HR  I T  E  I  3 » 1 2  ) ( NST (  I  )  ,1  =  1, NR) 
l?  FORMAT! • ONST • / ( IX, 2013) I 
RFwINf)  F 
IF(NF)  13,13,14 

14  1)0  IS  J  =  1  , N F 
R  F  A  0  (  F  ) 

15  COMTINOF 

13  P  F  AO ( F ) ( PP  (  I  )  ,  I = 1  ,  M P P ) 

WR11E(3,1F)(PP(1  )  ,1  =  1,3) 

IF  FORMAT!  'OPP* /( 1 X.FFll. 4)  ) 

NNH=N*NH 

REAO! 1 , 17) ( Fid  ) , 1=1, NNH) 

17  FORMAT (4F14. 7 ) 

HR  I  Tfc  (  3, 1R ) (FI  (  I  )  .  1  =  1, NNH ) 

1R  FORMAT! • OF  I •/( 1X.4F14.7)  ) 
mn=n*n 
00  IN  J=1,NN 
Jl= J+NZ 

R< J)=RFAL(PP( J 1  )  ) 

10  CONTINUE 

00  20  .1=1, OR 
Jl  =  ( J-l  )*N 
00  21  1  =  1, N 

.13=  (  I  -1  )*N 
.12  =  J  1  +  I 
RI I J2)=0. 

00  22  K  =  1 , N 
J4= J3+K 
JS=J1+K 

RI  (  J2)=RI(  J2)+R(.)4)fF!(.JS) 

22  COMTINOF 
21  COMTINOF 
20  CONTINUE 


M 

•  I  •  0 

nn  ?  3  i *  |  ,n** 

J7» l j- I  l «n 
nn  ?<.  I  « l  ,  i 
J !  •  J  I  •  I 
WCS( j| I «0. 

J<.»  (  I  -  I  1  »N 
Oil  7b  k  »  |  ,n 
.H*.I7<« 

J  b  *  J  ♦  K 

RCSI  Jl  »=R0<>I  .111  *F  I  l  .tb»*R|  (  J3) 
7b  CUN  II  MIE 
2 <*  CONTINUE 
23  CONllNilt 

WRI  TH  3t2M  IRCSI  I  )  »  I  =  I »  J 1  ) 

?b  FORMAT)  'DKCS'/I  1  X  ,  7  e  1  1  .4  11 
Dll  27  I  =  1  ,NH 
VC  C I  ) =0. 
ji  =  <  i-n*N 
1)0  28  K  =  1  ,N 
-12  =  J 1  ♦K 
J3=NV+K 

VCU  )«VC(  !  »+F  I  (  J7  )  *>PP  ( .13  ) 

78  CONTINUE 

VI  (  I  )=REAL< VCI  I  )  ) 

V?  (  I  ) =  A  I M  AG ( VC (  I  1  ) 

27  CONTINUE 

WRITE<3,29)  (VC(I  )  »  I  =  1  » NP ) 

7B  FORMAT ( ,CVC'/(1X,7E11.4)  I 
C1=C2*BK*BK 
00  30  J=1,N8 

IFINSTU)  .  E  0 . 0  1  GIT  TO  30 
J2  =  0 

DO  31  K  =  1  » -I 
J3= ( K- 1 )*J 
Jb  =  K 

DO  32  I  =  1  *  K 
.12  =  J  2  + 1 
J4= J3+ I 

RC( J4)=RCS( J  2 ) 

RC< J5)=RC< J4) 

Jb=J5+,l 

32  CONTINUE 
31  CONTINUE 

b 2  CALL  LINER  ( -l»RC  ) 

51  =  0. 

52  =  0. 

00  33  K  =  1  ♦ .! 

A  1  (  K  )  =  0  . 

A  7  (  K  )  =  0 . 

.13=  (K-l  )  *J 
0(1  34  1  =  1  ,  J 
Jl= J3+ I 

A1(K)=A1(K)+RC( J1 >*V1<  I  ) 
A?(K)=A2(K)i-RC(  J1)»V2<  I  ) 

34  CONTINUE 

51  =  S1  +  ABSI A1IK1  ) 

52  =  S2  +  ABS( A2(K  )  I 

33  CONTINUE 
GC=0. 
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ihi  35  k  =  i .  j 

GC=GC+Vl  I  K  I*  A)  IK  I+V2I  K  )+A2!K  ( 

36  CUN  T I  I'll l(- 

s  -i  gc=gc*c  l 

IF  <  SI *S2.Nb .0. )  Oil  ril  36 

gr=gc 

r.=  l . 

011  1(1  38 
36  S1=0. 

52  =  0. 

53  =  0. 

on  34  k = i » j 

S1=S1+V1IK|»A1IKI 
S?  =  S2  +  V7(K)»'A2<K) 

S3  =  S3  +  Vl(K)’*A2(K) 

34  CONlINUb 

A=(S1-S2)/(2.*S3) 

SA=SORT< A*A+1.  I 
C=-A+SIGN< SA, S3  ) 

GR=C1*( S1+C*S3  ) 

38  (ill  40  K  =  !,.l 

ALP(K)=A1(K)+C’!‘A2IK) 

A2 ( K  )  =-A2 1 K ) 

40  CONI  INUF 

WR I Tb ( 3,45 )  GR.GC 

46  FORMA  1 (  »0*,5X, «GR»  ,9X,*GC'/1X,2E11.4) 

WR I Tb  <  3,46) <  ALP  (K  )  ,K=1  ,  J) 

46  FORMA  1  (  1  ORFAl  ALPHA  FOR  MAXIMUM  GA I N • / < 1 X , 5E 1 4 . 7 ) ) 

WR I Tb ( 3,47 ) I A1 <K) , A2IK) ,K=1,  J) 

47  FORMAT (  •OOOMPLEX  ALPHA  FUR  MAXIMUM  G A  IN ' / ( 1 X , 4F 14 . 7  I  I 
III!  48  K  =  1  , N 
C  R  1  (  K  )  =  0  . 

CR2(K  )=0. 

J  1  =  K 

DO  44  1=1, J 

CR 1 ( K ) =CR 1 1 K  )  +F I (J1I1A1III 
CR2(K  )=CR2(K)+FI(J1 )*A2<  I  ) 

J 1 = J 1+N 
44  CONTINUE 

CIIR<K)=CR1(K)-C*CR2(K) 

48  CONTI NUF 

WR  I TF( 3,60) (COR(K)  , K  =  1 , N ) 

50  FORMAT (  'ORFAL  PORT  OUANTJTIFS  FOR  MAXIMUM  GA IN • / I  1 X ,4F 14. 7 )  I 
WRITE  <3,51  HCR1(K),CR2(K),K=1,N) 

51  FORMAT!  'OCOMPLFX  PORT  QUANTITIES  FOR  MAXIMUM  G A  I N ' / ( l X , 4F 14 . 7  I  I 
30  CONTINUE 

9  CONTINUE 
STOP 
FNO 

tllATA 

l 

4  2  31066  161061  0. 1463495E+00 

0  0  1 

-0.133R342E+00  0 .4326460E+00  0 . 84 186H4F +00  0.9949999F+00 
-0.6078172E+00  0.9999999E+00  0. 80544 70E+00-0. 645750SE+00 
0, 999999 6E +00-0. 5 37 368 2F +00  0 . 1 3679 2 l F-0 1  0 . 7396740E -01 
J.STOP 
/«■ 

II 


m 


:$ 
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PRINTED  OU1PU1 
NE  =  1 

N  N6  NH  NPP  N 7  NV  RK  t 

4  2  31066  161061  0. 1963495E+00 

MSI 

0  0  1 
PP 

0.1917E-03  0.306 IE -02  0 . 39 33E-03-0 . 3197E-02-0 . 637 1 E -03-0 . 2576E -03 
El 

-0. 133R342E+00  0.4326460E+00  0. R41H6fl4E+00  0. 9999999E +00 
-0.607ft 172E  +  00  0.9999999E  +  00  0 . ft  0544 70E  +  00-0 . 64 5 7505E +  00 
0.9999996E+00-0.53736A2E+00  0. 1 36792 1 E-0 1  0 . / 396740F-0 1 

PCS 

0. 1H57E  +  03-0.2441E-03  0 . 1 047E  +  03-0 . 5« 1 7E-04-0 . 2R23E -03  0.6  226E  +01 
VC 

0. 7 265 F +01 -0. 67  IRE  +  01-0. 537 7 E  +  01  0. 3362E  +  01-0. 1 13  7E  +  01-0. 24R6E+01 

OR  GC 

0. 1554E+01  0 . 2440E+01 

PEAL  ALPHA  EOR  MAXIMUM  GAIN 

-0. 109049 3E +01  0.9509 159E+00-0. 1 26505 7E+02 

COMPLEX  ALPHA  FOP  MAXIMUM  GAIN 
0.39 12077 E-0 1  0.361759HE-01-0.5133666E-01-0.3209720E-01 
-0. 1H25690E+00  0.3992ftfl2E+00 

REAL  PORI  QUANTITIES  FUR  MAXIMUM  GAIN 

-0. 130R260E+0?  0 . 7 27 7 1  29E+01 -0 . 325 1 R9HF+00-0 . 264027RE +0 1 

i 

COMPLEX  PORT  QUANTITIES  FOR  MAXIMUM  GAIN 

-0.  1566013E  +  00  0.41 R9557E+00  0. 6369549E-01 -0. 23 1 0 1 05E  +  00 
-0.1091181E-01  0.1 006477E-0 1  0 . 5P 76725E-0 1  0 . ft 643705F -0 1 


i 

i 
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Next,  Che  eec  of  reel  pore  currente  listed  under  Che  heeding  "Reel 
pore  quantities  for  maximum  gain"  In  the  printed  output  of  tho  optimum  gain 
program  was  fed  Into  the  modal  reaonance  program  of  section  IV.  The  loads 
appearing  In  the  printed  output  of  the  modal  resonance  program  were  then  fed 
Into  the  o/X2  versus  angle  program  of  aactlon  V  with  the  variables  NQ  and  N6W 
changed  so  as  to  score  o/X2  on  record  6  of  data  set  6.  The  following  printed 
output  resulted. 
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P8IN1PD  OUlPlIl  IIP  IMP  MlllAI  MPSIINANCP  PPIIGHA8 

NF  ■  | 

N  NA  NH  12 
A  2  I  1  A 


PM 

0.  | >4 1  /(- -03  0.  3<:MF  -02-0.  39 3 3F -03-0.  3|9  7P-02-0.  A  3 7  I F -0  3-0. 25 7AF -03 
CUP 

-0.1 3082A0F«02  0.7?771 29E ♦ 0 1 -0 . 325 1 89«F *00-0 . 2a A0278f »0 1 


*1. 

0. 25M30A2P *03  0. 1 AlA 3AAF*03  O.A5223a7F«OA-<).  30O593mF*03 

PPIN1M)  OUlPItl  IIP  I H  P  SICMA/L  AMHIIA  SOIIA8F0  VERSUS  ANCLE  PHor.PAH 
NP  NAW 
I  2 

N  NA  NO  NS  NL  Nf  NA  NPAI  N1V  UJ  J7  |3  HK 

A  2  2  1?  I  IAS  73  2  1  lA  3A1  10  0.  1*PA3A9DP*00 


0. |9| 7F-03  0.  30AIF-02-0.  3933F-03-0. 3197F-02 


0.25830A?F*03  0. 1 A1A3AAP «G3  O.A5223A7F *0A-0. 300593«F »03 
SCAUFBINC.  PAlTPPN  HIP  IMF  11H  Stl  IlF  LOADS 


angle 

AFALII 1 

1  NAG  IF  1 

IF  I 

SIG/ILA8)**2 

0.0 

0.H379F«00  -0.312AF-0I 

0.8335F  *00 

0. A9A  7F  *00 

30.0 

0. 72 |9F  *00  -0 . A  78  1 F  «00 

0. AA59F  *00 

0, 7A97F  *00 

AO  .  0 

-0.2990F*00  -0.  70 (OF *00 

0 . 7A2 1 F  *00 

0.580«E»00 

00.0 

-0.A32AF-01 

0. 3059P*00 

0. 3090F  *00 

0.95A7E-OI 

170.0 

0.23H2F*00  -0.A|97E*00 

0.AA3AF  *00 

O.AAOlE«IO 

150.0 

-O.a800F*00  -0.3753F*00 

0. 77A7E  *00 

0,A033F«00 

1  AO.  0 

-0. 75 10F  *00 

0.AA09E-01 

0. 7523E  *00 

0.5A5«E*00 

210.0 

“0. A AOOt ♦ 00  -0. 3753F  *00 

0. 77A7C*00 

O.A033F»00 

2A0.O 

0. 23A2F  *00  -0. A 197F  *00 

0,AA3AE»00 

0.A501E  »00 

270.0 

-0.A32AF-01 

0 . 3059E  *00 

0.3090E*Or 

0.95A7E-01 

300.0 

-0. 2990F  »00  -0.7010F*00 

0. 782  IF  *00 

0.5808F  *00 

330.0 

0.72 1 9F  *00  -0.A7H1F*00 

0. 8 A59t  *00 

0. 7A97E  *00 

3A0.0 

0.8379F+00  -0.312AF-01 

0.833SP*00 

0,  A9A  7E  *00 

0.0 

0.8  379F  *00  -0.312AE-01 

0.8  335E  *00 

O.A9A7E»OII 

30.0 

0 . A2A 7F *00  -0.A57AE-01 

0.A3A5E  *00 

0.A025E  *00 

AO.O 

O.A39AF*00 

0 . 2297E  *00 

0. A9A0F *00 

0. 2AA0F  *00 

90.0 

O.RASAF  *00 

0. 10AAF*00 

0 .9  5 1 AF  *00 

0.905  IE ‘00 

1  20.0 

0. A27 l F  *00  -0. 32AAF  *00 

0. 70AI E  *00 

O.A9HAE«00 

150.0 

0 . 58  7 AF  *00  -0.A2A9F-0I 

0. 5907E  *00 

0.3A89E»00 

1  AO.O 

0. 75 1  OF  *00  -0.AA09F-01 

0. 7573F*00 

0. 5A59E  *00 

210.0 

0 . 58 7 3F *00  -O.A251F-OI 

0.5R07F  »00 

0.3A89E*00 

2A0.0 

0. A27 1 F  *00  -0.32AAF*00 

0. 70A1 F  *00 

0.A9MAE*00 

270.0 

0.RA5AF  *00 

0. 10AAE*00 

0.951AF  *00 

0.905  IE  *00 

300.0 

0. A  39AF  *00 

0.2292E*00 

0.A959£*00 

0. 2A59F  *00 

330.0 

0.A?B7E«00  -0.R533F-0I 

O.A3a5F»00 

0.AO25E  «00 

3A0.0 

0.8329F*00  -0.3128F-01 

O.A335E  *00 

0.A9A7E  ♦  0:> 

SIC. 

0.A9A7F*00  0.A952F  *00  0.  A9A7F  *00  0.A993E+00  0.7027F*00  0.7071E*0() 
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VII.  BACKSCATTERINC  VERSUS  FREQUENCY 

The  program  (o/X2  versus  frequency  program)  of  this  section  requires 
the  port  parameters  Zc,  Voc,  and  ?°*  •  u  at  the  frequencies  of  lntereet. 

As  a  preliminary  step,  the  Impedance  matrix  program,  the  excitation  vector 
program  and  the  port  parameter  programs  on  pages  45-69  of  [7]  were  run  after 
Insertion  of  the  statements 

DO  200  I  -  1,  KV1 
V(I)  -  C0NJG(V(I)) 

200  CONTINUE 

just  after  statement  37  in  the  main  program  of  the  excitation  vector  program. 
These  3  additional  statements  change  the  Incident  field  from  a  plane  wave 
traveling  in  the  minus  x  direction  to  a  plane  wave  traveling  in  the  positive 
x  direction.  The  impedance  matrix  program  and  the  port  parameter  program  were 
run  with  the  G  level  compiler  but  the  excitation  vector  program  wae  run  with 
WATFIV.  It  has  been  observed  that  the  particular  propagation  constant  read 
in  as  0.1963495  is  printed  correctly  by  the  WATFIV  compiler  but  is  printed  as 
0.1963494  by  the  G  level  compiler. 
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MINIMI  •  KJ 1  Ml  •  1  (It-  IMPEDANCE  M  A  T  H  I  *  PRIIGRAM 
NF  N6  NP  NM  HAD 

ir  h  44  6  o. 16oooooe*oo 


P  X 


2.0706 

1.H117 

1  .bb24 

1.2441 

1.0353 

0. 7765 

0.5176 

0.2588 

0.0000 

-0.2588 

-  0  ,  b  1  7  6 

-0.776b 

-1 .03b3 

-1.2441 

-1.5629 

-1.8117 

-2.0706 

-1.5529 

-1.8117 

-2.0706 

-I. 03b 3 

0.0000 

1.03b3 

2.0706 

1.8117 

1.5529 

-4. 6b«7 

-4.3499 

-4.1411 

-3.1058 

-2.0706 

-1 .03b3 

0.0000 

1  .0363 

2.0706 

3.1058 

4,1411 

4.3904 

4.6587 

-4.14H 

-4.3444 

-4 . 6bR  7 

-4.4176 

-b. 1764 

-5.4352 

-5.6940 

-5.4528 

-6.2117 

-5.1764 

-4.  1411 

-3. 1058 

-2.0706 

- 1 . 03b3 

0.0000 

1  .0353 

2.0706 

3.  "08 

4.1411 

5.1764 

6.2117 

5.9628 

b . 6440 

6.4  332 

b , 1 7  64 

4.4176 

4.6587 

4.3949 

4.141) 

2.0706 

1. 1038 

4.1411 

3.RR23 

3.623b 

3.3646 

3.1058 

2.8470 

2.68R2 

2.3244 

2.0706 

1.0353 

0.0000 

0.0000 

-1 .03b3 

-2.0706 

-2.3244 

-2.688? 

-2.8470 

-3.1058 

-3.3646 

-3.6?33 

-3.HM23 

-4.1411 

-3.1 ObH 

-2.0706 

PY 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

o.oooo 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0,0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

o.oooo 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

o.oooo 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

o.oooo 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

o.oooo 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

p? 

7.7274 

6.761b 

S.74b6 

4. 8246 

3.8637 

2.8978 

1.9319 

0.9659 

0.0000 

0.4659 

1.4314 

2.R47R 

3.H637 

4.H296 

5.7956 

6.7615 

7.7274 

5.7936 

6.7615 

7.7274 

7.7274 

7.7274 

7.7274 

7.7274 

6.7615 

5.7956 

17.3867 

16.4207 

15.4548 

15.4548 

lb.4b4R 

lb.4b4R 

1 b , 4b4H 

1 b ,454H 

15.4548 

15.4548 

1  6 ,4548 

16.4207 

17.3867 

15.4548 

16.4207 

17.3R67 

18.3626 

14.3185 

20.2844 

21.2504 

22.2163 

23.18?? 

23.1822 

23.182? 

23.1R22 

23.1R22 

2*.  1R22 

23.1822 

23.1822 

23.1822 

23.182? 

23.18?? 

23.1822 

23.182? 

22.2163 

21 .2b04 

20.2H44 

14.3185 

18.3526 

17, 3867 

16.4207 

16.4548 

15.4548 

15.4548 

lS.4b4R 

14.4RR4 

1 3.b230 

12.5570 

11.5911 

10.6252 

9.6503 

8.6433 

7.7274 

7.7274 

7.7274 

7.7274 

7.7274 

7.7274 

8. 64 3 3 

9.6593 

10.625? 

11.5411 

12.5570 

13.5230 

14.4PR4 

lb.4b4H 

lb.4b4R 

16.4548 

ll 

1  1H  2  7  40  64  H? 

HKK 

0.1  ,707966*.j0  0 •  16684 716*00  0 . 1  767  1451: +00  0.  1  8260506  +00  0. 1  865320F*00 
0. .9046906*00  0.19242246*00  0.19438596*00  0.19536786*00  0.  1  9634946*00 
r. 19733126*00  0.19831306*00  0 ,200?766t +  00  0 . 2022344b *00  0.20616706*00 
0.2 1009 396*00  0.21 598456*00  0. 225801 RE +00  0.  ? 356 1 v?E *00 

IMPFOANCF  MATRIX  OF  ORIIFR  3R 

0.19626*01-0.66566*03  0.19436*01  0.28976*03  0.18866*01  0.54136+02  0.12R6F+00 
0.14626*01-0. 16066+01-0.1383E+0? 

|  MR6 OANC  6  MATRIX  (16  (IKI)hR  3H 

0.22146*01-0. 61 3B 6* 03  0.21R46+01  0.27446*03  0.21176+01  0.51486*0?  0.14476*00 
0. 18806*01-0.  1 7446  *0 1 -0, 1 3326*02 

PLUS  17  MORE  APPEARANCES  OF  IMPEDANCE  MATRIX  OF  IIRD6R  3R 

PRINTED  OUTPUT  OF  EXCITATION  VECTOR  PROGRAM 

Nb  N6  NP  NW  NT  NPAT 
14  2b  44  6  1  1 


UK 

(t.  |  K  70744K  401) 
o.  1404440K  +on 
O.  197331  7M0O 
o. > i oooiomoo 


o.  i #»#*«**?  1  h  *ito 
o.  14?4??4K  ♦oo 
0.19831  oik  +no 
0.21448441  400 


n.  r/47i44Mnn  n 

o.  14438401- #on  o 
0.2002744K ♦OO  I) 
o.224hoi hi  4(to  o 


1  8  2404(1*  •III' 
I*<si4  /hi  »nrt 
20 2219*<K  *  /mi 

?  144  l'l  it  <in» 


0.  i  mi,*,  i?ni  4ft« 
0. 1*11  14441  <00 
II.  70',  |  *>  II  I  t  (Ml 


l'  l  A 


7 


M  X 


7.0704 

1.8117 

1 .4479 

1.7941 

1.0141 

(1.  7744 

0.41 74 

(I.248M 

0.0000 

-O.  74*<m 

—  0.4 1 74 

-0.7  7  44 

-1  .034  1 

-1.7441 

-1.4429 

-1.8117 

-2 . 0  704 

-1.4429 

-1.8117 

-2.0704 

-1 .0383 

0.0000 

1  .014  1 

7.0704 

1.8117 

1  .442-1 

-4 .448  / 

-4 . 3444 

-4.1411 

-  1.  |04- 

-7.0/04 

-1 .0343 

0.0000 

1.0343 

2.0704 

3. 1048 

4.1411 

4.  1  <44 

4 , 4*1  H  7 

—  4.1411 

-*.  .  3944 

-4 .448  7 

-4.4174 

-4.1744 

-4.434? 

-4  .4040 

-4.44?x 

-4.211/ 

-4.1 744 

-4.14M 

-1.1  OSH 

-7.0704 

-1  .034  1 

0.0000 

1.0343 

2.0704 

3.1048 

4.1411 

4.  1  //.** 

'*.211/ 

S  .  4  S  2  M 

4.4440 

4.4  14? 

4. 1 744 

4.41 74 

4.4487 

4 . 14  49 

4.I4H 

7.0/04 

1.  |04l 

4.14)1 

3.8871 

3.4734 

3.  1444 

3.1048 

2.8470 

2.488? 

2. 3744 

?.o  704 

1  ,n  14  « 

0.0000 

0.0000 

-1.0143 

-7.0704 

-2.3244 

-2.4882 

-2.8470 

-3.  10411 

-3,3444 

-  1.  +  /  14 

-3.887  1 

-4.1411 

-3.1048 

-7.0704 

MY 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.01100 

O.OO/lO 

0.0000 

0.0000 

0.0001' 

0.0000 

0.0000 

0.0000 

(1. 000(1 

o.oono 

o.ooon 

o.ooon 

0.0000 

0.0000 

O.i  00(1 

0.0000 

0.0000 

0.0000 

0.0000 

n.oooo 

0.0011(1 

o.nnn< 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0. ooon 

0.0000 

0.0000 

7/.  00(11* 

0.0000 

0.0000 

0. 0000 

0.0000 

0.0000 

n.oooo 

n.oooo 

0.0000 

O.OiK'O 

OIK* 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

n.oooo 

0.0000 

- .OOOn 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0,0000 

(1,000' 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0. 0000 

0.0000 

0.000(1 

0.0000 

".noon 

0. 0(100 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

n.COOO 

0.0  100 

o.onoo 

0.0000 

0.0000 

0.0000 

0.0000 

M  7 

7.7274 

4.7414 

4.7444 

4.8744 

3.8437 

2.8978 

1 .4114 

0.944*1 

o.ooon 

0.944 u 

1.4314 

7.8478 

3.8437 

4.8744 

4.7944 

4.7414 

7.7274 

4.7414 

4  .  7  /»  |  4 

7.  7?  74 

7.7774 

7.7774 

7.7774 

7.7774 

4.7414 

4.7944 

17.  1847 

14.4707 

14.4448 

14.4448 

1S.4S4H 

14.4448 

14.4448 

14.4448 

14.4448 

1 4.444H 

14.4  4**  8 

14,4707 

1  7.3X4/ 

14.4  84*. 

14.4707 

17.3847 

18.3474 

14.3184 

20.2844 

71 .2404 

72.2143 

23.1*7? 

23. 1 “7? 

?  4 .  1 8?? 

23. 1H?7 

71.187? 

73. 1 H?7 

73.1872 

2 3.1822 

73.18?? 

21.187? 

7  1.18?? 

71. lx?? 

2i.l“22 

77.21M 

21 .7404 

70. 7844 

19.1184 

18.3424 

17.3X4? 

14.4207 

1S.444H 

14.4448 

1  4.444*- 

1S.4S48 

14 .4884 

13.4730 

12.4470 

11.44  1  1 

10.424? 

4.4491 

8.4431 

7.  7?  7<* 

7.  '2  74 

7.7774 

7.7774 

7.7774 

7.7274 

8.44  33 

4.4441 

10.424? 

1 1 .491  1 

12.4470 

13.4/  10 

14.4884 

14.4448 

14.4448 

14.4448 

LI 

1  |H  27  AO  49  8? 

0.  104448  3(  +  01-0.41?4414t-*03 
-n.  ?i  Moa7  h  400-0. <.MbH?nF4no 
0.1017711K -03  0. 3040442K-02 
0.  lOOOOOOf-401  0.8444449M00 
0.4)  1  3738K*00  0.41  1  7  7H‘>K400 
0,140  2279 fc  +00  0.  1494A10K+00 
0.  |94??44(:+0l— 0.4444  74PK*03 
0.??l3A40fc*fl1-0.4n747  3K*n3 
(I ,  ?4  7448  4t  +  n  1  -  0  •  4744  2  34F  +  0  3 
0. 2444884 6 +01-0. 4448408F +03 
0.774  11  ?4I- +01-0.44  28  34  4F  9  03 
0.?H7  77O9K+Ol-0.S30.i?84X  +  03 
0.24  348  7  71: +01 -0.47424441: +03 
(>.2494436t+01-(1.4  1M?97  IK  403 
0. 30247. 34t  401 -0.^1*>3404t  403 
0 . 304498  3t  +  0 1 -0. 4 l ?44 1 4F  +  03 
0.  108 7398F+01-0.4094 /03F  +  (H 
0.31  1  79041: +  01-0. 4067 143K +  03 
0.31  794  34b  401  -0.  *0108471-403 
0.324  1  421K+01-0.4944404F+03 
0.  33474  9  4X  +  01 -0. 484*1  (147K +03 
(1.349  477RK  +  0 1 -0.4  7442801: +03 
0.349248  31  +01-0. 4494244F  +  03 
0.40372071-401 -0. 434 3047K  403 
(1.4  384?  141- 4.',  1  -0.4  148  844K4  03 


war 
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-o.  31  aii  .oo-n.  aovsi  *no-n.A?i  S7  *00-0.70777  *oo-n.  amqkf  *no-o.  1S3A7  ♦oo-o.s  I  3?7  400 

-O.MA?*  -  0|  -0.  AMOH7  4  00-0.  |S3A|  »oo 
m'llln  OIIIPlII  ID-  P(|M  |  M  AM  AMI-  1 1  M  PM00MAM 

St-  0)  H  U  N av  NA/  Nap  svt 

I V*  |  4M  A  ?•,  -70  I  I 

tin 

0, | S  7 P70A7  4  00  0. |AA10717*00  0. I 7A7)aS7*00  0.  I  M  ?  A0S07  •  00  0.  I  1AS3707  400 
0.  lOOASOOf  *00  0. 10?A??A7 *0(  0.  |<4A  31S07  400  0.  |  QS^A  7Hh  »00  0  .  I  <4  A  3AQAt  400 

0.  |9733l?7  400  0. |0H3| 307*00  0.  700?  7 AAf  4  00  0.  ?0??3‘497*00  0.  70A  | A  707  *00 
0.71000307*00  o.?isomam*oo  0.77MOIH7.0O  o.?  3S<  1077  400 

NA 

A  MM??  I  7  3  S  A  7  H  10  n  l?  |A  IS  |A  17  |«  l't 

711  7 1  73  ?A  ?■>  ?A  ?7  ?M  ?0  30  31  17  33  3A  3S  3A  37  31 

nv  N/  miv  h i  i  n|sc  Msr.  *7SC  <*ior.  Avuf  one 

l  i  o  o  o  i  i  o  i  i 
O.10S77 .oj-o.si7m  «ni  n. 30117*01  0.71077403 

-0. 7  1  a  77  >  00-0.  A  A3  A7  *00-0. 370A7 *00-0. 3S I  37*00 
O.I9177-C3  0.30A1* -07-0.10337-03-0.11077-0? 

0.10007*01  0.«AA07‘O0  0. SP007  *00  0.00007*00 
0 . A | | A7  »00  O . A  1  IMF  *  no  0 .A | 307  ♦ 00  0.A1S07«0C 
O. 1 A077  *00  0.  1  AO A7  *00  0,  17  OAF  *00  0.|7??7*00 
0.  10A77  *01-0.  ASSA7«03  M.lOA37*f)l  0.71977*03 
0.??|A7*01-0.A13M7*0»  0.71107*01  0.77AA7*03 
O.7aaP7*01-O.S7aa»  *03  O.7AAO7*01  0.7A0<»f»O3 
0,  ?AA  7  401-0.SSS97.03  O.  ?A  1  ?f  *01  0.7l3At*03 
n.?7Al7«OI-O.WR7*03  0.77737401  0.?ahO7*03 
n.?«7M»-*oi-0.  *>3037*03  0.  ?h  377  *01  0.?aaS7*03 
0.79377*01 -O.S?A37 *03  0.7HOA7401  0,?A?At-*03 
0. 79977*01-0. S1H37»03  0.?9S?t:*01  0.7AO3E403 

n.  30777. 01-0. S1SA7*03  0.?ot»17*01  0.73037*03 
0.10A77*01-n.*>l?S7*03  0.  301  17*01  0.731?£*03 
O. 30177*01-0. *»OOAt*03  0 , 30A07  *01  0.73777*03 
0.31 1*7. 01-0. SOa77*03  0.30707*01  0.73A?7*03 
0.3170t-. 01-0. 1)01  17*03  0.31707*01  0,?3A37*03 
0.37A77  *01-0.A0r,AF*03  0.31107*01  0.73?A7*03 
O.33a71»01-O.aia17*33  0.33117*01  O.??RA7*0! 

O.3AOA7*01-0.a7aa7«03  0.3a3S7*01  0.7?SH*03 

o.  )AO37*01  -O.A«>Oi7*03  0.3A?S7*«H  n. 77007*03 
O.A03?7*01-O.A3A37*03  O.3os?t*01  0.71717*03 
O.A31A7*01  -O.A  (A07  *03  0.A70  |7*01  0.70S07*03 


-0. 3 IA17  *00  O.A0S37*00-0.A?lbF*00  0. ?o?77 *00-0. A AOHt *00  0.1S3A7*00 
} 

0. 10A?7*01  -0.ASSA7 *03  0 . 10A  37  *01  0.71077*03  0.)mia7*01  0.SAI37.0? 

7 

0.77JA7  *01-0. A1  317  *03  0.71107*01  0.?7AAt*03  0.71177*01  O.S1A17*0? 

PUIS  17  fOiF  APP74P1NC7S  07  /.  Tm7  7INH  7  UN7S  117  PP1N770  (HlTPIM  A«7 
PP 

0. 10 AA 7 -03  0. 1  3737-07-0. 1 7  717-0  3-0.  33 377 -07-0.7AO77 -03-0. A0A07 -03 
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According  Co  Che  preceding  prlnctd  output,  Impedance  matrices  are  or 
records  7  Co  25  of  data  sac  6.  The  excitation  vectors  and  port  parameters 
era  on  records  26  and  27  respectively.* 

In  the  o/A2  vert us  frequency  program,  the  activity  on  data  sets  1 
(punched  card  Input)  and  6  (direct  access  Input  and  output)  Is  as  follows. 

REA1> (1,10)  NC.N6K 
10  F0RMAT(20I3) 

DO  12  JC  "  .1,  NC 

READ(1,13)  N,  NF,  N6,  NL,  NIV,  NZ,  NV,  13,  B 
13  F0RMAT(5I3,  314,  E14.7) 

READ (1,15) (BK(I) ,  I  -  1,  NF) 

15  FORMAT (5E14. 7) 

REWIND  6 

SKIP  N6  RECORDS  ON  DATA  SET  6 
NPP  -  NV  +  (NF-1)*I3  +  N+l 
READ (6)  (PP(I) ,  1-1,  NPP) 

DO  22  JL  -  1,  NL 
READOL,  23) (XL(J) ,  J-l,  N) 

23  FORMAT (5E14. 7) 

22  CONTINUE 

12  CONTINUE 

SKIP  N6U  RECORDS  ON  DATA  SET  6 
WRITE<6) (SIG(I) ,  I  -  1,  J7) 

Virtually  all  of  the  main  program  Is  Inside  DO  loop  12.  Of  all  the  propagation 
constants  BK,  B  la  the  particular  propagation  constant  at  which  the  load  react¬ 
ances  XL  are  evaluated.  The  N-port  parameter  Zg  at  the  .T  n  frequency  resides 

In  PP(NZ  +  (J-1)M3+1)  through  PP(NZ+(J-1)*I3+N*N)  while  Voc  and  F0C  •  u  reside 

1*0 

in  PP(NV+( J-l)  *13+1)  through  PP(NV+(J-1)*I3+N+1)  .  For  the  open  circuit  impedance 
formulation,  NIV  +  0  whereas  NIV  -  0  for  the  dual  short  circuit  admittance 
formulation.  The  JL**1  set  of  reactive  loads  is  read  into  XL  inside  DO  loop  22. 

Minimum  allocations  are  given  by 
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COMPLEX  C(N*N) 
DIMENSION  LR(N) 


In  Che  subrouclne  LINEQ  end  in  the  mein  progr.  m  by 


COMPLEX  PP(NV+(NF-1)*I3+N+1) ,  Z(N*N),  V(N) 
DIMENSION  BK(NF) ,  C2(NF),  XL(N),  ND(N) ,  D(N) , 
NC 

SIG(  l  NL*NF) 

JC-1 


The  index  JF  of  DO  loop  25  denotes  the  JFS  frequency.  DO  loop  26  puts 

Z  in  Z.  DO  loop  27  adds  the  reactive  loads  to  the  diagonal  elements  of  Z. 

If  the  ratio  of  the  magnitude  of  the  ND(J)  reactive  load  to  the  magnitude  of 

the  ND(J)t^1  diagonal  element  of  Z  is  greater  than  10,  then  DO  loop  33  divides 
th  th  ^ 

the  ND(J)  row  and  the  ND(J)  column  of  Z  by  the  square  root  of  this  ratio. 
Statement  35  Inverts  Z.  DO  loop  37  is  similar  to  DO  loop  33.  DO  loops  33  and 
37  scale  [11]  the  matrix  Z  to  avoid  excessive  round  off  error  in  the  subroutine 
LINEQ.  DO  loops  33  and  37  have  no  net  effect  in  the  absence  of  round  off  error. 


DO  loop  39  puts  Voc  in  V. 


DO  loop  40  accumulates  V°C[Z  +Z  ]  *Voc  of  (68) 

of  [7]  in  EL.  In  statement  43,  PP(J3)  is  F  •  u  of  (68)  of  [7],  The  constant 

k2n  **°  '*'r 

C2(JF)  -  — rW  is  necessary  to  obtain  o/Azwhich  is  subsequently  stored  in 

4tt  /z 

SIG((JL-1)*NF+JF)  for  the  JLth  set  of  loads  and  the  JFth  frequency. 


There  are  two  sets  of  reactive  loads  in  the  sample  input  data.  The  first 
set  of  reactive  loads  resonates  the  set  of  port  currents  which  radiates  the 
pattern  synthesized  by  the  pattern  synthesis  program  of  section  II.  The  second 
set  of  reactive  loads  resonates  the  set  of  port  currents  generated  by  the 
optimum  gain  program  of  section  VI.  The  38  values  of  o/A2  appearing  in  the 
printed  output  are  put  on  record  28  of  data  set  6.  The  center  frequency  men¬ 
tioned  in  the  printed  output  is  the  frequency  at  which  the  propagation  constant 
is  B.  Note  that  o/A2  at  the  center  frequency  is  the  same  as  o/A2  at  180°  on 
the  pattern  printed  by  the  o/A2  versus  angle  program  of  section  V. 


I 
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LI  SI  INI.  Ilf-  |HL  SIGMA  OVER  LAMBDA  SODAREI)  VERSOS  EREOOFNf.V  PROGRAM 

//  i  003*.  pe  ,  20s.  n .  'maot; ,  jne  •  .meg idn.2ook 

//  P*P(.  Ko’HV 

//GlI.P  10AF001  nil  PSNAMP  .EE003A.REV1  .0!  SP'IILO.ONI  I  "2  31  A  ,  1 

/  /  VI IL I  If  I  -',FR*S0000A,l)CH.(RFfFW«VS.BL,<SI/E  ■  251  A  ,  LPF  C  L  ■  2512  ,  < 

//  NnFNO"  t  I 

//GO.SVSIN  DO  • 

f.M'H  mAh!/ .  1 1  hh i  ,PAr,fs«<.n 

(  1h|S  PROGRAM  CALLS  IMF  NATH  I  y  |  NVE  H  S  1 1  If  SUBROUTINE  LINEO  LISlPI) 

C  Hi  1H  Imp  SIGMA  OVfR  LAMBDA  SOWAR  E  0  VERSUS  ANCLE  PROGRAM  flE  SECTION  v 

CWMPLf  *  II.PP  (  7VB  I  .7  II  OO  I  .VI  A  I  .EL  .CWW  ,E  .Cl  IN, If. 

WIDENS  I  (IN  HK I  I V  )  « C 2 (  IV  I . XL  I  A  I  ,  M)l  A  I.IIIaI.SIGIaORI 
IA»37A.730 
PI «3. 1AI5V3 

Cl».25*FTA/SQRTIP|*P|»P|  I 
ll>  10.  .  I  .  I 
R  )  AO  (  I  .1(11  NT. NAM 
TO  FORMA  1 ( 201 31 

MR I  1 P  (  3 «  1  1  I  NC.NAm 
11  FORMA  1 (  ' 0  NC  HAH' / 1 > • I  3.  I  A  I 
.17.0 

Oil  |2  JOI.NC 

READ!  t  . 1 3  I  N.NF .NA.NL .Nl V.N7  .NV  ,  I  3.M 
13  EIIHNA1(5|3.3|A,E1A.7> 

NR I1E 1 3. 1A  I  N.NE,  S.NL.NIV.N7  ,  NV  .  I  3 , P 
1A  F  ( IR  M  A  1  (  *  0  N  NF  Ni  NL  NIV  N?  MV  |  3  •  ,  A  *  ,  •  H  •  /  I  1 ,  A  I  3 ,  A  I  A  .  F  1 A  .  7  I 
UEAD(1«I5I(HK (11.1*1. NF) 
l5  FORMA  1  (  5F  I  A  ,  7  I 

MR  ME  (  3,  \A)  (MKT  |  )  ,  |-  1  ,NF  I 
1A  FORMAK  '0HK.'/<lX,5EtA.7ll 
REMIND  A 
IFINAI  17.17. 1H 
1«  1)0  IV  J ■  1  « ElA 
READ! A  I 
IV  CONTI NDF 
17  NP 1 *N* l 

NPP=NV*  (  NF  —  1  )*|3*NP1 
RFAOI Al |PP( I  I  ,  I  >1  ,NPP) 

WR I TE ( 3 « 20 1  I PP  ( I 1,1*1,21 

20  FORMAT!  'OPP'/dX.AFll. All 
DO  21  J=1 .NF 

C2(  J  I  *C  1  *RK  ( .1 1  *  RK  ( .1 1 
RK!Jl*hK< Jl/B 

21  CONTINUE 
NN=N*N 

1)0  22  .IL  =  1  »NL 
RE  AD  I  1.23)  (XL(.I)  ».l*  1  »N  I 
23  FORMAT ( 5E 1 A. 7  I 

WRITEI3,2AMXL<J)..I*1.N> 

2A  FORMAT!  '  OXL  '  /  I  1  X  .  5E  1 A  . 7)1 
MR  ME !  3 . 5  1  I  .IL 

51  FORMAT! • OBACKSC AT  TER  I NG  VERSOS  RATIO  OF  FRFOWENC V  TO'/'  CENTER  FR 
1F00ENCY  FOR  THF',I3.'TH  SET  OF  LOADS') 

WR I TE  (  3 . 5  2  I 

52  FORMAT! •  FREQ*  ,  3X  ,  'REAL ( E 1 ' ,5X . ' |M AG ! F ) ' , 7X , ' I E I  • » AX , • S 1 0/  (  L AM  I  * 
1*2'  ) 

K7=N7 
KV  =  NV 

00  25  JF  = 1 »NF 
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Ilf  I  ?b  ,l»|  iNN 
.17  *  J  *K  7 
/  (  .1 1  »PM  .17  I 
?4  CMN 1 1 Nlit- 
.17*0 
.11*1 

1)11  77  .1*1. N 
I  *■  (  XI  (  .1 II  7H  ,7H  ,  7n 

?h  *• *1 ( J I/kk I JF  ) 

Oil  1(1  10 

74  *«ltU  J  I  •  MK  (  Jh  | 

10  M*AHS(  X  I  /f AHSI7  (JIM 

/ 1  jii>;  (  jii  *•»•  * 

JI*JI *NP| 

IMSI-IO.I  77.77.JI 

n  j7* j?* i 

Mil  J7  I*  J 

III  J7I°  I  ./S<»R  I  t  S|  I 
77  CON  1 1  Mil 

lMJ7.fo.oi  r.d  in  is 

Ml  33  J ■ I . J 7 
J I ■NIX Jl 
J 3* I J 1  - 1  l»N 

111.  34  I  *1  .N 
J1*J3«  1 

7  (  J3  I  */  (  Jl  I  •()(  Jl 
/(  Jl)*/(JI  I  •**(  J  I 
J  I  •  J  I  *N 
14  CONllNUf 

33  CONI  IN' It 

3S  CALL  MNF0IN.7  I 

IFU7.E0.0I  f.ll  111  34 

Ml  37  J-1.J7 

Jl  *NI)I  J  I 

J3«( Jl-I  l*N 

01)  30  I  ■  I  ,N 

J3* J3* 1 

7 ( J3 I *Z  <  J3 1 *01 J I 
2 ( J 1 1*7  I Jl I »0( Jl 
Jl-Jl+N 
3«  CON  T I  Ni IF 
37  Cl  IN  3  I  Nile 

34  IK)  34  J  *  | ,  N 
J  7*  J ♦K  V 

VI J 1  * PP ( J?  I 
34  CONlINtlfc 

el*o. 

!))•  40  J*1  ,N 
J  3  * ( J  —  1  )*N 
CUR*0. 

IK)  41  1*1, N 
J4«J3*I 

C'iR«CllR*7(  J4»*V(  I  1 
41  C0N1INIIF 

fl*fl ♦C»IR*V (  J  I 
40  C0N1IN0F 
J3«KV+NP1 
IF(NIV)  43.47,43 
47  F*(PP(  J3I-EL  l*C?l  .IF  I 
Oil  10  S3 
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<.3  f » (  ppi  .H  i  *ki  *r  ?<  je 1 

*>3  .1  7  «  J  7  *  I 

SI  f.  I  17  I  *7  •  Cl  IN  Jf.  (  E  I 
M- Sum  I  Mi;  I.i7  )  i 

WR I  IE  I  J.SA I  HKI  .IK  I  .  E  .n  ,  SIM  .17  I 
SA  ElIRMA  1  (  I  X  ,1-6.3  .AM  ?.A  | 

K/«K/«| 3 
KV*KV« I  3 
?•)  CONTINUE 
7?  Cl  IN  1  I  NUt' 

I  7  CONTINUE 

I  E  INAW I  A  7 .A  7  , AM 
AM  nil  AW  J  •  |  ,NAW 

«nuiM 
AW  CUNTINIIE 

A7  mm|  1E<  All  SIC.I  II.  |i|  ,J7I 
R«l  7K  3.  SOM  Sir,  l  |  I.  |  ■  J  ,M 
SO  EOHMA  II  «OS  |  r.  •  /  I  I  *  ,  At  11  .A  I  I 
S  ini’ 

F  MI* 

mill 
i  o 

I  |w  7k  7  I  |*  37  A 7  O.IVASAVSE*O0 

0.  |S707waE  .00  n.  1AARV7  1E«00  0 . 1  7*  7  1 A  Sfc  *00  0.  1  h?A050E  »00  0.  |HAS3?0E  ♦00 
(I.  |W0ASW0E«00  O.  lW?AP?Sf  *00  0.1WA3MA0f»00  0. I WS3A7ME.OO  0.  |WA3aWsE»00 
0. |w  7  3  3 1 ?E  *00  0.1W«31 3OE.O0  O.?O077» AE *00  0. 707?  Jwwf  *00  0.70MA70E  *00 
().?I00W3WE«00  O.JISWMASfc.OO  0.?>3HO|ME«00  0.?  3SMW3fc.no 
0.37?7w|OE«03  O.A71O|?3E«O?-O.10MaaamE»0S  0 . Aa 3 1 w 3AE • 03 
O.^SM JOA?E«03  0.1MA3AAE.03  O.AS?73A7F*OA-0.300SW3ME*03 
AS  1 1 J  M 
/• 

// 

P«  INlEI)  UUlPlll 

NC  NMI 
l  0 

H  NE  NA  NL  Nl V  N7  NV  13  8 

a  |i  71  7  1  1A  37  A?  0. lWA3AWSE.no 

MR 

0.  |S707WAE«00  ll.IAAHW71E.00  0  .  I  7  A  7  I A  SE  ♦  00  0. 1 RPAOSOE »00  0.  I  HAS3?0E  .00 
0. 1W0ASW0E.00  0.1W?A??SE*00  0. IWA3HAOE.OO  0. 1WS3A7HE *00  0.  IWA3AQSE *00 
0. |W73  312E*00  0.1WH3I 30E.00  0. 70C? 7 AAE »00  0. 707? 3WWE ♦ 00  0. 70A | A70E .00 
0 ,  ?  1 00W  39  E  *00  0.?|SWHASE*00  0.  ??SMO I  ME »00  0. ? 3SA l W 3E .00 

PP 

0. 10AAE-03  0. 137 3E -07-0. 1 77ME-03-0. 3337E-0? 

U 

0.37?7V|OF*03  O.A7  1Ol?3F*n?-0.10EAAA«E.0S  0 . aa 3  I w 3aE *03 

HAf  KSC  A 1 1EM  I  NO  VERSUS  RATIO  IlF  FREQUENCY  in 
CENlER  FREQUENCY  FI1H  THE  1  IH  SET  HE  LIIAOS 
ERFO  REAL  I  E  1  IHAMFI  I  E  I  S  1 0/  1 L  *  7 

II.HOO  -0.A733E-07  0 . W?0?E -0 1  0.«??7E-0|  O.MSlAC-O? 

0.850  0. 10WAE .00  O. 1 A 1 ?E .00  0.17R7F.00  0.31WAE-01 

O.WOO  -0.?0'.7E»00  -0.7AWHE.00  0.2<?3WE*00  0.10AWE.00 

O.W30  -0.13AOE.00  0.838AE-0I  O.ISWRE.OO  O.?SS?E-01 


102 


0, 

.950 

-0, 

>  5H49F 

-01 

0 

.970 

0, 

.75R7F 

-01 

0. 

.  4  RO 

0, 

.  22 1 2F 

+  00 

0 

.94  0 

0, 

.4  T 12E 

+  00 

0, 

.494 

0, 

,579ftF 

+  00 

1 

.000 

0, 

,5ftft4F 

+  00 

1 , 

,004 

0. 

.4275F 

+  00 

1 

.010 

0. 

,2ft4ftF 

+  00 

1 . 

.020 

0. 

,  5254F 

-01 

1 

.030 

-0. 

,  17Mf 

-01 

1 . 

,050 

-0. 

,  ft25  7F 

-01 

1 , 

.070 

-0. 

,5923F 

-01 

1 . 

.100 

-0. 

,  ft7ft5F 

-01 

1 , 

.150 

-0. 

,  6774F 

-01 

l. 

.200 

-0. 

,  77ftOF 

-01 

0.1  T32E+00 
o.  2ft58e+oo 

0. 3001F  +  00 
0.212RE+00 
0.4335E-01 

-o. irsof+oo 

-0. 34R0E+00 
-0.3471E+00 
-0. 3340F+00 
-0.2544F+00 
-0.1  54ftE+00 
-0. 1088F+00 
-0 • 7074E-0 1 
-0.440ftF-01 
-0 . 2744F-0 1 


0.  182HF+00 
0.2754T  +00 
0.372RF+00 
0.51711+00 
0.58121:+  0 
0.54ft3F+00 
0 . 55 l 2E  +00 
0 .4 77?F  +00 
0. 3447F+00 
0 . 2ftOOF  +00 
0.  1714E  +  00 
0.1 2R4F+00 
0.47  8  RE-0 1 
0.80R2F-01 
0.H231E-01 


0. 3342F-01 
0.7ft39E-01 
0. 1390E+00 
0.2ft74E+00 
0.337«P+00 
0.355ftF+00 
0. 303RF+00 
0.2277F+00 
0. 1 1RRF+00 
0.5750E-01 
0.2939E-01 
0. Iftft2fc -0 1 
0.9580E-02 
D.ft532E-02 
0. 5774E-02 


XL 


0 . ?5fl  30ft 2E  +03  0. lftl*34ftF+03  0.ft5223ft7E+04-0.3005938F+03 


HACK  SC  A  T 1 FR I  NO  VFRSUS  RAT  10  OF  FRFOIJENCY  TO 
CENTER  FREQUENCY  FOR  THE  2TH  SET  OF  LOADS 


FRFO 

RFAL(F) 

I  KAO ( F  1 

1  £  1 

SIR/ (LAM  1**2 

O.ROO 

0. 2503E-02 

0. 84 14F-01 

0.8417F-01 

0. 70P5F-02 

0.H50 

0.57H3F-01 

O.H8ME-01 

0. 1058F+00 

0.1120P-01 

0.900 

0. 1205F+00 

0.4872E-01 

0. 1300E+00 

0. 1590E-0 1 

0.430 

0.151  3F  +  O0 

-0.7158E-02 

0. 1515E+00 

0. 2245E-01 

0.950 

0.1814E+00 

-0.585 l F-0 1 

O.1754F+00 

0. 3075E-0 1 

0.970 

0. 1433F+00 

-0. 1755E  +  00 

0. 22ftftE+00 

0.51 34E-01 

0.980 

0.9220E-01 

-0.2734E+00 

0.28R5F+00 

0.8324E-01 

0.990 

-0.1124E+00 

-0.4315E+00 

0.4459E+90 

0. 1989E+00 

0.995 

-0.4309F+00 

-0.4385F+00 

0.8148F+00 

0.37ROE+00 

1  .000 

-0.751 OF+OO 

0.441 0F-01 

0. 7522F+00 

0.5ft59i+( 

1.005 

-0.3540E+00 

0.4775E+00 

0.59441  +00 

0.3533E+00 

1.910 

-0.3072F-01 

'.4151F+00 

0.4 1  ft 2F  +  00 

0.  1732E  +  00 

1  .020 

0.151 7F+00 

0.21 ftOE+OO 

0 . 2ft 34F  +00 

0. ft9ftftE-0 1 

1.030 

0. 1H4RF  +  00 

0. 1034R+00 

0.21 ft4E  +00 

0.4ft81F-01 

1.050 

0.1952E+00 

-0.2047E-01 

0.14ft4F+00 

0.3R5ftE-01 

1.070 

0.1746E+00 

-0. 1 003E+00 

0. 2014E+00 

0.4057E-01 

1.100 

0. 1 1H9F+00 

-0. 1854F+00 

0.2202F+00 

0.4R50E-01 

1.150 

-0.2448E-01 

-0.2612F+00 

0. 2523E+00 

O.ftRROE-Ol 

1.200 

-0. 2025F+00 

-0. 2384E+00 

0.3128F+00 

0.4787E-01 

SIC, 

0.H514F-02  0.3194E-01  0.1049F+00  0.2552E-01  0.3342F-01  0.7ft39F-01 
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VIII.  PLOTS 

The  synthesized  pattern,  the  scattering  pattern  of  the  wire  triangle 
loaded  by  the  set  of  reactive  loads  that  resonate  the  set  of  port  currents 
which  radiate  the  synthesized  pattern,  and  the  scattering  pattern  of  the 
wire  triangle  loaded  by  the  set  of  reactive  loads  that  resonate  the  set  of 
port  currents  which  radiate  the  optimum  gain  pattern  are  plotted  from 
records  4,  5,  and  6  of  data  set  6  by  the  plot  program  on  pages  104-110 
of  [7]. 


PRINTED  Din  PUT  OF  PLOT  PROGRAM  ON  PAGES  104-110  OF  REFERENCE  7 

NF 

3 

NT  NE  NEP  NS  NA 
14b  4  3  12  3 

N 1 

0  0  11 
N2 

1  3  4 

N3 

OOl 

NSRL 

4  0  0 

SCAl 

0.2000E+01  0. 2000E+0 1  0.2000E+01 
0.3057E+01-0.5125E+03  0.30UE+01  0.23R2E  +  03 
-0.2142E+00-0.4636E+00-0.370AE+00-0.3513E+00 
0.1917E-03  0.3061E-02-0.3933E-03-0.3197E-02 

SIG 

0.1000E+01  O.flAAOE+OO  0.5000E+00  0.0000E+00 

NT  NF  NEP  NS  NA 
145  2  2  12  4 
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N  1 

1  1 

M2 

1  2 
M3 

0  1 


NSRL 
4  4 

SCAL 

0 «  2000E  +  01  0 . 2000E  +0 1 

0.3057E+01-0.5125F+03  0.3011E+01  0.2382E+03 
-0.2142E  +  00-0.4ft36E+00-0.370<SE  +  00-0.3513E+00 
0.1917E-03  0.30ME-02-0.39  33E-03-0.319  7E-02 
0 . 1000E+01  0.R660E+00  0.5000E+00  0.0000F+00 

SIG 

0 . 65  37  E+00  0.6533E+00  0.6524E+00  0.650RE+00 

NT  NE  NEP  NS  Nfi 
145  2  2  12  5 

N 1 

1  1 
N2 

1  2 
N3 

0  1 

NSRL 
4  4 

SCAL 

0«  2000F+01  0.2000F+01 

0.3057E+01-0.5125E+03  0.3011E+01  0.23R2E+03 
-0.2142E+00-0.4636E+00-0.3706E+00-0.3513F+00 
0.1917E-03  0.3061E-02-0.3933E-03-0.3197E-02 
0. 1000E+01  0.8660E+00  0.5000E+00  0.0000E+00 
0. 6537E+00  0.6533F+00  0.6524E+00  0.6508E+00 

SIG 

0.6947E+00  0.6952E+00  0.6967E+00  0.6993E+00 
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A  new  plot  program  (o/A2  varaua  frequancy  plot  program)  haa  baen  written 
to  plot  tha  data  atored  on  raeord  28  of  data  aat  6  by  tha  a/ A2  varaua  fraquancy 
program  of  aactlon  VII.  Tha  o/A2  varaua  fraquancy  plot  program  accapta  input 
data  on  data  aata  1  (punched  carda)  and  6  (direct  acceaa)  in  the  following  manner. 


R£AD(1,10)  NC 
10  FORMAT(20I3) 

DO  12  JC  -  1,  NC 
READ(1,13)  NF,  N6,  NE,  NEP,  B 
13  FORMAT(4l3,  E14.7) 

READ(1,10)(N2(I),  1-1,  NEP) 
READ(1,10)(N3(I),  1-1,  NEP) 
READ(1,17)(BK(I),  1-1,  NF) 

17  FORMAT (5E14. 7) 

REWIND  6 

SKIP  N6  RECORDS  ON  DATA  SET  6 
NEF  -  NE*NF 

READ (6) (SIG(J) ,  J-l,  NEF) 

12  CONTINUE 


Virtually  all  of  the  program  la  inaide  DO  loop  12.  The  quantity  o/A2  for 
the  aet  of  loada  and  the  Itl1  frequency  la  in  SIG((J-1)*NF+I)  where  J-1,2,...NE 
and  I»1,2,...NF.  Among  the  propagation  con8tanta  BK,  B  la  the  propagation  conatant 
correaponding  to  the  center  frequency.  The  Ith  curve  to  be  plotted  ia  that  of 
o/A2  for  the  N2(I)£^  load  aet.  If  N3(I)  >  0,  the  pen  draw8  the  horizontal  and 
vertical  axea  and  movea  to  the  next  frame  after  plotting  the  I1'*1  curve.  If 
N3(I)  <_  0,  the  pen  neither  plota  any  axl8  nor  move8  to  the  next  frame  after 
plotting  the  Ith  curve. 


Minimum  allocationa  are  given  by 


DIMENSION  N2 (NEP),  N3(NEP),  BK(NF) ,  SIG(NE*NF),  Y(NF) 

Do  loop  19  prepares  the  horizontal  coordinatea  BK  for  plotting.  The 
index  L  of  DO  loop  24  denote8  the  L curve  to  be  plotted.  DO  loop  25  put8 


I 
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the  vertical  coordinates  In  Y.  Statement  31  plots  o/X2  versus  frequency  by 
drawing  straight  Unas  between  data  points.  Statement  27  draws  the  vertical 
axis.  The  logic  between  and  Including  statements  28  and  29  puts  tha  scale  on 
the  vertical  axis.  Statement  30  draws  the  horlsontal  axis.  DO  loop  26  puts 
tha  scale  on  the  horizontal  axis. 


/ 
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A 


I 


t 

k. 

4- 


l 

1 

i 

> 


II  I  00  14 .  LL  ,  I  ,  I  .  ,  10  I  ,  'MAUI  7  ,  JOE  •  ,urr,  | IIN -  |  4  0* 

II  WM.  I.  MAUL  l  iNf^  r  I IM  V  |  It-  hill  NOW  AND  WILD  111  MLR  H  IP  H  L  A  f,  K  INDIA  IMP 
//  I  .If  MW  IC.r  l(.  ,  PARM.LIW  1  •  "MAP* 

//FI|M  1  .SYS  IN  DD  • 

D|  ML  NS  |  IN  A«L  A  I  4001  ,N2I  1001  .Ml  I001.RKI50I  .  X  |  (  4.  |  ,  VI  I  A  I  ,  S  |f.|M)«  I 

II I  PL  NS  |  IIN  Y(SOI 

CALL  MI.Hl  III 

CALL  ML01 S I  AH L  A , 4  00 1 

HLAIM  1,10)  NC 

10  FIIUMA1I20I3) 

WMIlLn.'ll  Nf. 

11  FtIHMAll  'ONO«'.I3> 

Dll  1 2  JC-I.NC 

M  F  AD  I  1,131  NL,NO,NF.NFP,R 
13  FORMA  1  ( 4 l 3  *  E 14 . 7 1 

WH  I  1 F I  3 , 1  a  |  NF ,NO,NE ,NEP  ,R 
1A  FOHMAIUO  NF  N  A  NE  NE  P  '  ,  A  X  ,  •  H  •  / 1  X  ,  3  I  3  ,  I  A  ,  E  1 A  .  7  1 
R  E  A  0  I  I • 10) IN2II 1 . 1=1  ,NFM  1 
WR I  IF ( 3, l S 1 <N2( 1 1 , 1 *1 ,NEM) 

15  FORMAT!  •  ON?  '  /  (  IX, 20131  1 
READ!  1  ,  101 (N3<  I) ,  I =1 ,NEM  I 
wR I Tfc ( 3, 16  1 1 N3 ( I  1,1  =  1 ,NEP ) 

10  FORMAT!  * 0N3 ' / ( 1 X, 201 3 1  I 
RFADI 1 , 17 1 (HK( I ) , I =1 ,NF 1 
17  FORMA  T ( 6F 14. 7 1 

WR  I  TE  (  3 *  1R  M  HK I  I  ) ,  I  =  1 , NF  I 
1H  FORMAT! • ORK 1 / ( IX,5E14.7 I ) 

Sl  =  l()./F 

DO  19  J= 1 , NF 

PK(J)=S1*RK(.I)-A. 

19  CONTINUE 
XI  I  1  1  =  1. 

XI  (21  =  1. 

Vlll  1  =  1. 

VI (21=7. 

X  l  (  3  1  =  1 . 

XI  I  41=0. 

Y  1 1  3  1  =  1 . 

VI <41=1. 

REWIND  0 

I F ( NO  1  20,20,21 

21  00  22  J= 1 , NO 
RFADI 01 

22  CONTINUE 

20  NEF=NF*NF 
RFAD!OHSir,(.l),.)  =  l,NEFl 
WR  I  TF ( 3, 23 1 <  SIG I J 1 , J= 1 ,0 1 

23  FORMAT! •OSIG'/I IX, OF  11 .A) ) 

0(1  24  L=  1  ,NEP 
J1=(N2(LI-1 1*NF 

1)0  25  J=  1  ,NF 
J  2  =  J  +  J 1 
Y ( J  1  *S I G ( J2 1 

I F  <  Y  < J) .LT.. 00001 1  Y(J)  =  . 00001 
Y(J)=6.+AL0G10(Y!J)  1 
25  CONTINUE 

31  CALL  LINEIHMl), Y(l), NF, 1,0,0! 

I F ( N3 ( L 1 1  24,24,27 


5  S»  ‘'4' 
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27  CALL  LINE!  XI  ,  Y  I  ,2,  1 ,0,0) 

?H  C  ALL  NUMBER  (  >14«  IO<  «0a  iO) 

CALL  SYMBOL  II., 7.,. 14,13,90.,-!  ) 

CALL  NUMBER!  .44,5.93, .14, 1 . ,0. ,0) 
call  symbol! i., 4.,. i4,i3,90.,-i » 

CALL  NUMBER  (  .52  ,  4  .93  ,  .  14  ,  .  I  ,0.  ,  1  ) 

CALL  SYMBOL (1., 5.,. 14, 13, 90., -1  ) 

CALL  NUMBER (  .40, 3. 93,. 14,. 01  ,0.,?) 

CALL  SYMBOL ( 1. ,4.,. 14, 13,90.,- 1 ) 

CALL  NUMBER  I  .?fl , ? .93 , . 14 , .  00 1 ,0 . , 3 ) 

CALL  SYMBOL! 1., 3.,. 14, 13,90., -1 ) 

CALL  NUMBER (  .14,  1  .93,. 14,. 0001 ,0., 4) 

?9  CALL  SYMBOL!  1  •  , ?•  ,  .  14, 13,90. ,-l ) 

30  CALL  L INF ( XI  (3) , Y1 (31 ,2,1 .0,0) 

XA  =  6. 

XH=1  .2 
00  24  J*  1 « 5 
XC*X A-.14 

CALL  NUMBER ( XC . . 7?, . 14, XB, 0. , 1 ) 

CALL  SYMBOL (X A, 1.,.  14,1 3,0., -1) 

X  A*  X  A- 1 . 

X  B  =  X  H— •  1 
24  CONTINUE 

CALL  BLOT ( 7 . ,0. , -3 ) 

24  CONTINUE 
12  CONTINUE 

CALL  PLOT ( 5  .  »  0. » -3  I 

STOP 

END 

/* 

//Oil  • F  T04F00 1  00  DSNAME  =  EE0O34.REVl,0lSP=OLD,UNIT*2314,  X 

//  VOLOME»SER=SU0004,nCB=<RECFM=VS,BLK$!ZE*2594,LRECL=2592,X 

//  RIIFNO=  1  ) 

//RO.SYSIN  DD  * 
l 

19  ?7  2  2  0.196 3495E+00 

l  2 
1  1 

0. 1570796E+00  0.  144B97 1E+00  0. 1767145E+00  0. 1R26050E+00  0. 1R65320E+00 
0. 1904590E+00  0. 1924225E+00  0. 1943B40E+00  0. 195347BE+00  0. 1943495E+00 
0. 1973312E+00  0. 19B3130F+00  0.2002766F+00  0. 2022 399E+00  0. 206 1670E+00 
0.21 00939E+00  0.2  159H45E+00  0.225R01RE+00  0 . 2356 193E +00 
/* 

// 

PRINTED  OUTPUT 


NC=  1 

NF  N6  NE  NEP  B 

19  27  2  2  0. 1963494F+00 

N2 

1  2 


N3 

1  1 


HK 

0. 1570796E  +  00  0. 166H97 1E  +  00  0. 1 747 14  5F+00  0. 1B26050F+00  0 . 1B65320E+00 
0.1904590E+00  0. 1924224E+00  0. 1943B59E+00  0. 195367BE+00  0. 1963494E+00 
0. 1973312E+00  0.19R3130E+00  0.2002766E+00  0.2022399E+00  0.2061670E+00 
0. 2 1009 39E  +00  0. 2159845E+00  0.225BC1BF+00  0. 2354192E+00 


/ 

i 


! 


SIC, 

0.85 14F-02  0.3194E-01  0.1049E+00  0.2552E-01 
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